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Abstract 


Objective.  The  goal  of  this  project  was  to  develop  a  practical  tool  for  optimizing  the  design  and 
operation  of  groundwater  remediation  systems  that  explicitly  considers  uncertainty  in  site  and 
remediation  system  characteristics,  performance  and  cost  model  limitations,  and  measurement 
uncertainties  that  affect  predictions  of  remediation  performance  and  cost.  The  project  was 
specifically  focused  on  chlorinated  solvent  contaminated  sites  with  dense  nonaqueous  phase 
liquid  (DNAPL)  sources. 

Technical  Approach.  The  method  is  based  on  a  semi-analytical  mathematical  model  to  simulate 
DNAPL  source  depletion  and  dissolved  phase  transport  in  response  to  natural  and  engineered 
conditions.  The  performance  model  is  coupled  with  cost  functions  for  thermal  source  zone 
treatment  and  enhanced  bioremediation.  Compliance  criteria  are  defined  by  statistical  rules.  The 
performance  model  is  also  coupled  with  an  inverse  solution  to  estimate  model  parameters, 
parameter  covariances,  and  residual  prediction  error.  A  stochastic  cost  optimization  (SCO) 
algorithm  is  used  to  determine  values  for  design  variables  that  minimize  expected  net  present 
value  cost  over  Monte  Carlo  realizations.  The  method  is  implemented  in  SCOToolkit  software. 
The  method  was  applied  to  two  well-characterized  sites  where  different  remedial  technologies 
were  used,  to  evaluate  its  ability  to  reduce  costs  and  improve  remedial  designs. 

Results.  Stochastic  cost  optimization  of  design  and  monitoring  variables  was  found  to  reduce 
expected  costs  by  50%  or  more  relative  to  conventional  design  methods,  and  to  substantially 
increase  the  probability  of  meeting  compliance  targets.  Although  additional  field  applications  to 
demonstrate  the  method  are  needed,  along  with  development  of  a  “user  friendly”  interface,  the 
method  was  shown  to  be  highly  effective  for  two  field  test  sites. 

At  the  first  site,  the  Fort  Lewis  East  Gate  Disposal  Yard  (EGDY)  site,  optimization  of  thermal 
source  treatment  indicated  a  need  for  a  much  larger  treatment  area  than  was  actually  employed, 
to  avoid  a  high  failure  probability  associated  with  source  delineation  uncertainty  based  on 
available  source  characterization  data.  Source  treatment  may  be  cost  effective  if  additional 
characterization  were  undertaken  to  reduce  source  zone  uncertainty.  The  method  was  also  used 
to  optimize  source  and  plume  bioremediation  at  the  EGDY  site,  using  whey  injection  without 
additional  source  reduction.  The  results  indicated  that  this  remedial  strategy  should  achieve  the 
maximum  contaminant  limit  by  2110,  with  a  93%  probability  of  success  when  using  relatively 
low  whey  injection  rates. 

The  methodology  was  applied  to  Dover  AFB  Area  5  to  optimize  enhanced  bioremediation.  The 
results  suggested  that  electron  donor  injection  may  be  needed  indefinitely  in  the  future  to  meet 
remediation  criteria,  as  a  result  of  the  uncertainty  associated  with  the  source  mass  estimates,  and 
the  difficulties  in  characterizing  and  treating  source  zones  without  disrupting  base  operations. 
Optimization  analyses  to  minimize  long  term  operating  costs  indicated  compliance  criteria  could 
be  met  using  only  five  of  the  current  ten  emulsified  vegetable  oil  injection  galleries,  with 
operating  costs  approximately  half  current  costs.  Recalibration  and  optimization  after  50  years, 
using  additional  data  from  this  period,  was  projected  to  further  reduce  operating  costs. 

Benefits.  The  results  indicate  that  SCOToolkit  has  a  potential  to  significantly  improve 
remediation  performance  and  reduce  costs.  Use  of  the  method  also  can  identify  the  most  critical 
data  gaps  and  the  uncertainties  that  will  most  affect  costs  and  performance  projections. 
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Chapter  1 
Introduction 


1.1  Background 

In  the  last  25  years,  there  have  been  a  number  of  studies  on  remediation  design  optimization, 
which  reflect  an  increasing  appreciation  for  the  importance  of  cost  optimization.  Teutsch  et  al. 
(2001)  emphasize  the  need  to  combine  physically  based  simulation  models  and  economic  models 
for  quantitative  decision-making.  Becker  et  al.  (2006)  reported  that  simulation-optimization 
methods  were  able  to  identify  solutions  that  cost  5%  to  50%  less  than  trial-and-error  results, 
translating  to  cost  savings  of  $600K  to  $10M  for  sites  studied. 

Much  of  the  early  optimization  work  focused  on  pump-and-treat  system  design  with  the 
objective  to  minimize  total  pumping  as  a  surrogate  for  operating  cost.  Gorelick  et  al.  (1984)  was 
the  first  to  combine  simulation  models  with  nonlinear  optimization  methods  to  identify  optimal 
design  variables.  McKinney  and  Lin  (1996)  argued  for  incorporating  fixed  as  well  as  operating 
costs,  leading  to  solutions  with  fewer  wells  and  higher  pumping  rates.  Regulatory  compliance 
(e.g.,  plume  containment)  has  been  variously  treated  as  an  optimization  constraint  (Wagner  and 
Gorelick,  1987;  McKinney  and  Lin,  1996),  as  a  “penalty  cost”  for  noncompliance  included  in  the 
objective  function  (Rizzo  and  Dougherty,  1996;  Chan-Hilton  and  Culver,  2005),  or  as  a 
component  in  multi-objective  programming  (e.g.,  Erickson  et  al.,  2002). 

Many  remediation  optimization  studies  disregard  uncertainty  in  simulated  performance  resulting 
in  overly  optimistic  estimates  of  cost  and  success  probability.  Wagner  and  Gorelick  (1987) 
considered  uncertainty  in  parameters  using  a  linearized  chance-constraint  formulation.  A  number 
of  studies  treated  uncertainty  explicitly  by  optimizing  the  expectation  of  system  performance 
over  an  ensemble  (Andricevic  and  Kitanidis,  1990;  Lee  and  Kitanidis,  1991;  Tucciarelli  and 
Pinder,  1991;  Wagner  et  al.,  1992;  Aly  and  Peralta,  1999;  Teutsch  and  Finkel,  2002;  Mugunthan 
and  Shoemaker,  2004;  Chan-Hilton  and  Culver,  2005;  Feyen  and  Gorelick,  2005). 

Although  pump-and  treat  methods  have  been  widely  used,  they  are  more  realistically  regarded  as 
containment  rather  than  a  remediation  measures.  This  is  particularly  true  for  sites  with  dense 
nonaqueous  phase  liquid  (DNAPL)  source  zones,  because  of  their  low  solubility  and  very  long 
persistence  (e.g.,  Cohen  and  Mercer,  1993;  National  Research  Council,  1994).  A  number  of 
remediation  strategies  have  been  developed  in  recent  years  to  deal  with  DNAPL  sites.  These 
include  source  mass  reduction  via  thennal  or  chemical  oxidation  or  surfactant  flushing  (e.g., 
Liang  and  Falta,  2008;  Heron  et  al.,  2009;  Thompson  et  al.,  2007)  and  introduction  of  electron 
donors  (EDs)  and  other  amendments  in  wells  or  trenches  to  enhance  in  situ  contaminant 
biodecay  within  the  source  and/or  dissolved  plume  (Wymore  et  al.,  2006). 

Effective  remediation  may  require  a  combination  of  technologies,  thus  complicating  the 
optimization  problem.  Furthermore,  a  more  complex  simulation  model  will  be  required  that  is 
capable  of  describing  the  individual  remediation  strategies  of  interest.  This  will  introduce 
additional  model  parameters  with  attendant  uncertainty,  which  will  propagate  to  performance 
predictions.  Mayer  and  Endres  (2007)  studied  cost  tradeoffs  associated  with  source  mass 
removal  and  dissolved  plume  remediation  measures  under  detenninistic  conditions.  Cardiff  et  al. 
(2010)  presented  a  semi-analytical  model  for  DNAPL  site  remediation  using  thermal  source 
reduction  (TSR)  and/or  ED  injection  for  enhanced  dissolved  plume  biodecay,  coupled  with  an 
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inverse  solution  to  estimate  model  parameters  and  their  uncertainty  from  measured  site  data;  a 
cost  module  to  compute  net  present  value  (NPV)  cost  to  meet  specific  objectives  and  constraints 
considering  fixed  and  operating  costs;  an  error  propagation  module  to  determine  expected  cost; 
and  an  optimization  module  to  detennine  design  variables  that  minimize  the  expected  cost.  A 
simple  test  problem  was  presented  to  illustrate  the  methodology. 

1.2  Objectives  and  Overview  of  Approach 

The  goal  of  this  project  was  to  develop  a  practical  tool  for  optimizing  the  design  and  operation  of 
groundwater  remediation  systems  that  explicitly  takes  into  consideration  uncertainty  in  site  and 
remediation  system  characteristics  and  inherent  model  limitations  that  affect  predictions  of 
remediation  performance  and  cost  -  with  a  specific  focus  on  chlorinated  solvent  contaminated 
sites  with  DNAPL  sources. 

The  general  approach  for  the  Stochastic  Cost  Optimization  Toolkit  (SCOToolkit)  described  in 
this  report  is  depicted  in  flowchart  form  in  Figure  1-1.  The  central  component  of  the  method  is  a 
semi-analytical  mathematical  model  to  simulate  DNAPL  source  depletion  and  dissolved  phase 
transport  of  a  target  chlorinated  hydrocarbon  over  time  in  response  to  natural  and  engineered 
conditions.  Currently,  the  performance  model  is  coupled  with  cost  functions  for  thermal  source 
zone  treatment  and  enhanced  bioremediation  using  ED  injection. 

In  the  inverse  modeling  mode,  historical  site  data  is  used  to  calibrate  the  simulation  model  and  to 
estimate  parameter  covariances  and  residual  prediction  error.  Forward  predictions  of  remediation 
performance  and  cost  are  performed  for  defined  remediation  strategies,  operating  rules  and 
remediation  criteria.  A  Monte  Carlo  (MC)  method  is  used  to  quantify  uncertainty  in  perfonnance 
and  cost  attributable  considering  uncertainty  in  model  parameters,  measurements  employed  for 
real-time  decisions,  and  cost  function  variables. 

Design  optimization  is  performed  to  determine  values  of  design  variables  that  minimize  the 
expected  value  (average  over  MC  realizations)  of  NPV,  which  may  include  “penalty  costs”  for 
failure  to  achieve  defined  remediation  objectives  within  a  specified  time  period. 

Chapter  2  of  this  report  describes  the  basic  transport  model.  Chapters  3  and  4  describe  the 
enhanced  bio  and  thennal  cost/performance  models,  respectively.  Chapter  5  outlines  calibration 
and  uncertainty  analysis  methods.  Chapter  6  presents  the  stochastic  design  optimization  approach 
and  Chapters  7-9  document  practical  applications  of  the  program.  A  program  user  guide  for 
SCOToolkit  is  available  in  a  separate  document. 
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Design  Optimization  Module - r  --  Model  Calibration  Module  - 


Figure  1-1.  Procedural  flowchart  for  SCOToolkit. 
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_ Chapter  2 _ 

Contaminant  Transport  Model  Formulation 

2.1  DNAPL  Source  Model 

Field-scale  DNAPL  source  dissolution  and  mass  depletion  over  time  is  described  by  the  model  of 
Parker  and  Park  (2004)  and  Park  and  Parker  (2005).  Considering  the  possibility  of  engineered 
manipulation  in  mass  transfer  kinetics,  we  describe  the  rate  of  contaminant  mass  dissolution  in  a 
source  zone,  J,  [M  T'1],  versus  time,  t,  by 


—  Fm,J cal 


f  M(t)V 

\Mcal  J 


(2.1) 


where  Jcai  =  J(t=tcai )  and  Mcai  =  M(t=tcai )  in  which  tcai  denotes  a  reference  time  selected  for 
model  calibration,  M  is  the  source  contaminant  mass  remaining,  (3  is  a  depletion  exponent  that 
reflects  the  DNAPL  source  “architecture,”  and  Fmt  is  a  time-dependent  dimensionless  mass 
transfer  enhancement  factor. 


Integration  of  the  source  mass  balance  equation  employing  (2.1)  as  described  by  Park  and  Parker 
(2005)  yields  source  mass  remaining  versus  time  after  the  release  date  t0  for  (3  1 

M(t)  =  -  F.,B(  1  -  m  -  (2.2) 

where  B  -  Jcal  /  Mpcal .  Considering  remedial  actions  at  dates  trem  /,  trem  2— trem  «  when  partial  source 

mass  removal  and/or  step  changes  in  Fmt  occur  and  stipulating  that  t0<tCaF  trem  1,  values  of  Mref, 
tref  and  Fmt  in  (2.2)  are  assumed  to  vary  with  time  as  follows 
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in  which  Mrefn  for  n  >  0  is  given  by 
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where  AMrem  „  is  the  mass  removed  from  the  source  at  time  trem  „  regarded  as  instantaneous. 
Sorenson  (2006)  reported  that  enhanced  source  zone  biodecay  caused  dissolution  rate 
coefficients  to  increase  by  factors  of  2  to  6  in  laboratory  studies  and  3  to  8  in  field  studies.  Parker 
and  Park  (2004)  also  have  shown  that  field-scale  dissolution  rate  coefficients  will  vary  inversely 
with  changes  in  source  zone  darcy  flux  (e.g.,  due  to  engineered  or  inadvertent  permeability 
decreases  due  to  amendment  injection).  Changes  in  Fmt  due  to  ED  injection  upgradient  of 
DNAPL  sources  are  described  in  Chapter  3. 
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2.2  Mapping  of  Groundwater  Flow  Field 

We  consider  multiple  chlorinated  hydrocarbon  (CH)  sources  that  may  occur  within  an  aquifer 
characterized  by  a  Cartesian  coordinate  system  in  field  mapping  units  (E,  N)  -  e.g.,  northing  and 
easting.  Flow  paths  will  generally  be  nonlinear.  To  apply  a  semi-analytical  solution  for 
contaminant  transport  to  mildly  nonlinear  flow  fields,  we  define  a  coordinate  transformation  for 
each  DNAPL  source  (and  later  for  ED  injection  galleries)  to  convert  from  field  coordinates  to 
linearized  local  coordinates  (x,  y)  and  back. 

For  each  source  j  we  define  the  origin  of  the  local  coordinate  system  to  be  at  the  center  of  the 
downgradient  plane  of  the  source  -  (E0,  N0)  in  field  coordinates.  A  streamline  may  be  drawn 
through  (Eo,  N0)  using  water  level  contours  and  dissolved  plume  data  and  selected  coordinates 
along  the  streamline  are  used  to  fit  a  cubic  polynomial  of  the  form  N  =  No  +  a(E-Eo)  +  b(E-Eo)2  + 
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c(E-Eo)  .  The  polynomial  function  is  used  to  define  the  (E,  N)— >(x,  y)  mapping  where  x  is  the 
distance  along  the  centerline  and  y  is  the  transverse  distance  orthogonal  to  the  centerline.  Given  a 
location  in  field  coordinates  (E*,  N*)  local  coordinates  can  be  found  as  follows  (Figure  2-1): 

1)  The  orthogonal  line  that  passes  through  (E*,  N*)  and  intersects  the  streamline  at  E  = 
Ecrossfliay  be  described  by  N  =  A+BE  where  A  =  N*  +BEcross  and  B  =  -(a+2b+3cX2)-1, 

2)  Solve  recursively  for  Ecross  using  Ecross  =  E*  initially  then  solve  the  cubic  equation  for  E 
where  the  orthogonal  line  intersects  the  streamline, 

3)  Compute  y  as  the  distance  from  (E*,  N*)  to  (Ecross,  Ncross),  and 

4)  Compute  x  as  distance  along  streamline  from  (E0,  N0)  to  (Ecross,  Ncross)  by  integrating  [1  + 
N’(E)2]1/2dE  from  Eo  to  Ecross  where  N’(E)=dN/dE.  Calculate  numerically  using  dE  =  10  m. 


Figure  2-1.  Curvilinear  streamline  in  (E,  N)  field  coordinates  and  mapping 
to  local  (x,  y)  coordinates  for  source  j. 
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2.3  Dissolved  Plume  Transport 

Dissolved  phase  transport  of  contaminants  emanating  from  DNAPL  source  zones  in  an 
unconfined  aquifer  is  described  by  a  2-D  vertically-averaged  semi-analytical  solution  that 
employs  the  source  zone  function  described  in  section  2-1.  The  solution  considers  linear  sorption 
and  first-order  decay.  Spatially-variable  decay  within  the  aquifer  can  be  described  with  up  to 
three  “zones”  at  different  distances  from  the  source  that  are  characterized  by  different  decay 
coefficients.  Zone  1  represents  the  region  x  <  Ln,  Zone  2  is  Ln<  x  <L23,  and  Zone  3  is  x  >  L23. 
Extension  of  the  model  to  simulate  ED-limited  biodecay  is  discussed  in  Chapter  3.  Leakage  of 
contamination  from  the  unconfined  aquifer  (unit  A)  to  an  underlying  semi-confined  aquifer  (unit 
B)  and  dissolved  transport  in  the  lower  aquifer  due  to  the  secondary  source  can  be  optionally 
considered.  The  solution  for  a  single  DNAPL  source  in  the  unconfined  aquifer  is 
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dr 


where  (x,  y )  are  local  coordinates  indexed  to  the  source,  t  is  Julian  date,  ta  is  the  source  release 
date,  x  is  a  dummy  integration  variable,  RA  =\+pkd/(, j>  is  the  retardation  factor  in  which  pkd  is  a 
dimensionless  sorption  coefficient,  vA  is  pore  velocity,  A  [  {  is  longitudinal  dispersivity,  A  is 

transverse  dispersivity,  LzA  is  aquifer  thickness,  LyA  is  width  of  the  source,  (f>A  is  porosity,  and  AAl 
is  a  first-order  decay  coefficient  for  Zone  i.  Integration  of  (2.4)  is  performed  numerically. 


For  Zone  1  A/U  =  Aai  and  computed  from  (2.1)  with  “normal”  calibrated  values  for  Jcai  and  Mcai. 
For  Zone  2  AAi  =  AA2  and  J ,,(/)  =  J  I2(J)  =  J  computed  from  (2.1)  with  values  for  Jcai  and  Mcai 
multiplied  by  a  scaling  factor  S2  defined  as 


g  _  CJx  Ln,y  0,t,AAI,JAl) 

Cm  (x  =  Ln,y-  0,  t;  AA2 ,  JAl ) 

For  Zone  3  AAi  =  An  and  ./  u(t)  =  J  n(t)  =  J  computed  from  (2.1)  with  values  for  Jcat  and  Mca/ 
multiplied  by  a  scaling  factor  S3  defined  as 


g  _  CAi(x  k23 , y  0 >t,AA2,JA2)  . 

CA\  (x  =  L23,y  =  0,t;  AA3 ,  JAl ) 

If  leakage  from  aquifer  A  to  B  is  considered,  Zone  2  can  be  treated  as  a  leakage  window  and  the 
decay  coefficient  is  computed  as  An  =  qv/(f>AkzA  where  qv  is  the  vertical  darcy  velocity  in  the 
leakage  “window”. 

The  solution  for  aquifer  B,  which  is  characterized  by  a  single  unifonn  decay  coefficient,  is 
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CB(x,y,t)  =  J 


Jatt-to-*) 


^^zB^yB^B  (^^B^LBVB^ 

erfc 


1/2 


exp 


f  V  (Rbx~vbt)2^ 


R , 


4RBAfBvBT  j 


<N 

-A 

+ 

1 

<N 

(g 

A' 

1 

1 

l  2(Atbvbt  /  RB)'n  ) 

erjc 

l  2(ATBvBr/RBy/2J 

(2.7) 


d  t 


where  Jg(t)  is  the  mass  leakage  rate  [M  T‘  ]  from  Aquifer  A  to  Aquifer  B  computed  as 


JB(t)  =  qvCAA  (2.8) 

where  A  is  the  area  [L“]  of  the  “window”  in  plan  view  (or  part  of  the  window  if  divided  into 
multiple  “panes”),  and  CA  is  the  average  concentration  in  the  window  (or  pane)  computed  in  the 

A  unit  aquifer  from  (2.4)  at  the  x  midpoint  of  the  window  or  pane  (i.e.,  compute  Ca  for  several 
points  on  a  transect  through  the  window  and  average  them).  The  computed  Jb  values  are  stored  at 
fixed  delta  t  increments  and  the  resulting  look-up  table  is  used  to  solve  (2.7). 


Contaminant  concentrations  resulting  from  multiple  sources  in  the  A  or  B  aquifers  are  computed 
by  superposition  of  the  individual  source  solutions  after  reverting  back  to  field  coordinates  as 


N 

1 ’  source 

C(E,N,t)  =  £  Cj(E, N,t)  (2.9) 

7=1 

Note  that  function  calls  on  the  RHS  of  (2.9)  require  (E,  N)— >(x,  y)  mapping  for  each  source  j  as 
illustrated  in  Figure  2-2. 


Figure  2-2.  Mapping  of  well  location  to  linearized  coordinates  for  two  adjacent  nonlinear 
streamlines  in  field  coordinates  for  solution  superposition. 
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_ Chapter  3 _ 

Enhanced  Bioremediation  Cost-Performance  Model 


3.1  Electron  Donor  Transport 

ED  injection  is  assumed  to  be  performed  in  one  or  more  injection  well  galleries  of  width  Led 
perpendicular  to  the  groundwater  flow  direction.  The  flow  field  for  each  ED  gallery  is 
characterized  by  a  polynomial  streamline.  Field  coordinates  (E,  N)  are  mapped  to  a  linearized 
local  coordinate  system  (x,  y)  for  each  gallery  with  the  local  coordinate  origin  at  the  center  of  the 
ED  gallery  as  described  in  Chapter  2. 

To  maximize  computational  efficiency,  ED  injection  for  each  gallery  i  is  approximated  as  a 
constant  mass  injection  rate,  MEDi  [M  T  ']  from  start  date  tEDo  to  stop  date  tEDf  with  MED  =  0 
otherwise.  Actual  aqueous  phase  ED  injection  will  generally  be  pulsed  to  reduce  well  fouling 
problems  and  injection  of  nonaqueous  phase  ED  will  be  performed  with  a  frequency  that 
depends  on  the  dissolution  rate  of  ED  material.  As  a  result,  temporal  variations  in  ED 
concentrations  will  occur  near  injection  galleries.  However,  since  these  variations  will  diminish 
markedly  with  distance  from  the  galleries,  modeling  ED  injection  with  a  time-averaged  rate  will 
generally  not  affect  ED  available  to  drive  biodecay  through  most  of  the  aquifer.  For  the  assumed 
step  function  injection,  the  ED  concentration  before  reactions  in  the  aquifer  attributable  to 
gallery  i  in  local  coordinates  is 


C™o-(x,y,t)  =  Q  for  t  <tEDo 

CEDi  ix’  y->  0  ~  C(x,  y,t  —  tEDo )  for  tEDo  <t<  tEDj-  (3 - 1 ) 

c'e°d >7  <X y> 0  =  c ( x ’ yJ-tuDo)  ~c yJ-  W )  for  t  >  tEDf 

where  C'EE.  is  the  ED  concentration  attributable  to  ED  gallery  i  prior  to  any  reactions  and 
C(x,y,f)  is  the  continuous  injection  solution  of  Domenico  (1987) 
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in  which  LED  is  the  ED  gallery  width  perpendicular  to  the  groundwater  flow  direction,  L:  is  the 
effective  aquifer  thickness,  q  is  the  darcy  velocity,  vR=q/(tf>RED)  is  the  mean  retarded  ED  plume 
velocity,  and  Ay  and  A  T  arc  longitudinal  and  transverse  dispersivity,  respectively. 


Kinetics  of  ED  reactions  is  approximated  assuming  the  fraction  of  injected  ED  concentration  that 
is  reactive  varies  exponentially  with  travel  time  as 


f  _ fy  yR  ^ 

C£?W,f)=  1-exp 

V  V  J 


CnEDi(x,y,t) 


(3-3) 
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where  (. x ,  v)  are  local  coordinates  for  ED  gallery  i,  cced  is  a  reaction  rate  coefficient  [T1],  v  is 
groundwater  pore  velocity  [L  T'1],  and  Red  is  the  ED  retardation  factor  [-]. 

For  multiple  ED  galleries,  the  total  ED  concentration  prior  to  reactions  C™"  at  a  given  location 
and  time  is  computed  by  superposition  in  field  coordinates  as 

C”ed(E,  N,t)  =  ^C^(E,N,t)  (3-4) 

i 

and  the  total  ED  concentration  available  for  reactions  is  computed  similarly  as 


C“”"  (E,  N,  1)  =  Xcs"  (E,  N,  t)  (3-5) 

i 

where  i  denotes  solutions  for  individual  ED  injection  galleries. 

The  foregoing  implicitly  treats  introduced  ED  as  an  aqueous  phase  material.  Nevertheless, 
nonaqueous  or  emulsified  ED  (e.g.,  various  vegetable  oil  formulations)  can  be  simulated  within 
this  mathematical  construct  by  suitable  adjustment  of  (Xed  and  Red  -  e.g.  by  calibration  to  pilot 
test  data. 


3.2  Electron  Donor  Reactions  with  Electron  Acceptors  and  Contaminants 

Biodecay  of  CH  species  is  assumed  to  be  limited  by  the  quantity  of  ED  species  relative  to 
electron  acceptor  (EA)  species  (Kamanth  et  ah,  2006).  To  estimate  the  attenuation  of  CH  due  to 
ED  addition,  a  superposition  method  is  used  that  is  analogous  to  that  described  by  Borden  and 
Bedient  (1986).  If  redox  reactions  occur  serially  in  order  of  decreasing  reaction  free  energy  (e.g., 
O2  >  NO3  >  SO4  >  Fe(III)  >  CH),  then  an  electron  balance  yields 


(-'CH  (serial)  =  maX 


0,  CCH  -  max 


0, 
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f'  avail  ,  n 
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H  nat  s-iH  \ 
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where  C’ch( serial)  is  the  aqueous  CH  concentration  after  serial  ED  reactions,  Cch  is  the  computed 
CH  concentration  before  ED  reactions  in  the  A  or  B  aquifer  as  needed  (as  described  in  Chapter 
2),  C‘"jf  is  the  aqueous  concentration  of  injected  ED  available  for  reactions,  Cf  "'"  is  the 
background  H-equivalent  ED  concentration  in  the  aquifer,  CfA  is  the  background  H-equivalent 

concentration  all  EA  species  in  the  aquifer,  f’cH  is  the  ratio  of  H-equivalent  to  actual  contaminant 
concentration,  f  Ed  is  the  ratio  of  H-equivalent  to  actual  injected  ED  concentration,  Rqh  is  the  CH 
retardation  factor,  and  Red  is  the  ED  retardation  factor  (no  retardation  is  assumed  for  EA 
species). 


H-equivalent  ratios  (Hstoch)  for  common  groundwater  ED,  EA,  and  solvent  species  are 
summarized  in  Table  3-1.  H-equivalent  ratios  for  EAs  are  estimated  as  f’cA  =  JeaI Eeff  wh ere  //.»  is 
the  stoichiometric  ratio  for  complete  EA  reduction  and  Eeff  is  fractional  energy  yield  for  the 
biologically-mediated  reaction  after  deducting  energy  consumed  for  cell  synthesis.  The  H- 
equivalent  ratio  for  CH  is  similarly  computed  as  f’cH  =  fed 'Eeff,  while  H-equivalent  ratios  for  ED 
species  are  computed  as  f  ed  =  fsoEeff,  since  ED  occurs  on  the  opposite  side  of  the  EA-ED 
balance  ledger. 
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Table  3-1.  H-equivalent  conversion  factors  for  selected  ED,  EA  and  CH. 


Species 

H stock 

_ EJ _ 

f* 

Chlorinated  solvents4 

PCE  (tetrachloroethene) 

0.058 

0.9 

0.064 

TCE  (trichloroethene) 

0.046 

0.9 

0.051 

DCE  (cis-l,2-dichloroethene) 

0.042 

0.9 

0.038 

VC  (vinyl  chloride) 

0.032 

0.9 

0.024 

Electron  donors 

Acetate 

0.13 

0.6 

0.078 

Butyrate 

0.21 

0.6 

0.126 

Ethanol 

0.26 

0.6 

0.156 

HRC 

0.12 

0.6 

0.072 

Lactate 

0.13 

0.6 

0.078 

Methanol 

0.19 

0.6 

0.114 

Molasses 

0.14 

0.6 

0.084 

Propionate 

0.18 

0.6 

0.108 

Vegetable  oil 

0.39 

0.6 

0.234 

Whey 

0.13 

0.6 

0.078 

Electron  accentors 

Oxygen 

0.125 

0.51 

0.245 

Sulfate 

0.081 

0.43 

0.188 

Nitrate 

0.083 

0.92 

0.091 

Iron  (II) 

0.018 

0.9 

0.020 

1  Hstoch  is  the  H-equivalent  ratio  computed  for  the  simple  chemical  redox  reactions  (Kamanth  et 


al.,  2006) 

2  Eeg  is  the  energy  conversion  efficiency  to  adjust  yield  for  cell  synthesis  for  relevant  microbial 
populations  (Rittman  and  McCarty,  2001) 

3  /’  is  the  net  H-equivalent  mass  ratio  =  HstochEefj  for  electron  donors  and  HstochIEeff  for  electron 

acceptors. 

4  Stoichiometry  based  on  conversion  to  ethene. 


If  reductive  dechlorination  of  CH  is  assumed  to  occur  under  anaerobic  conditions  with 
competition  among  microbial  populations  responsible  for  reduction  of  NO3,  SO4,  etc.,  an 
electron  balance  yields 
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1,  max 
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where  no  retardation  is  assumed  for  CH.  Assuming  that  actual  biodecay  can  be  approximated  as  a 
linear  combination  of  the  foregoing  pathways,  then 


C CH  ( mixed )  R serial^ CH  {serial)  +  ^  R serial)^ CH  {parallel)  (3  ^  ) 

where  Fseriai  is  the  fraction  of  reductive  dechlorination  that  follows  the  serial  pathway.  The  ED 
concentration  remaining  in  solution  after  reactions  with  EA  and  CH  may  be  computed  by 
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s~inet  _  s-inorx 

^  ED  ~  ^  ED 

which  is  used  to  determine  effects  of  ED  injection  on  mass  transfer  enhanced  in  (2.4). 

3.3  DNAPL  Source  Mass  Transfer  Enhancement  Due  to  ED  Injection 

The  current  version  of  SCOToolkit  considers  changes  in  DNAPL  mass  transfer  rate  by  a  factor 
Fmt  (see  Chapter  2)  due  to  ED  injection  upgradient  of  DNAPL  sources.  Whey  injection  studies  at 
the  Fort  Lewis,  Washington  East  Gate  Disposal  Yard  (EGDY)  DNAPL  site  by  Macbeth  and 
Sorenson  (2008)  indicate  an  approximately  linear  mass  transfer  enhancement  with  ED 
concentration  in  the  source  zone  (measured  as  COD)  that  may  be  described  by 

Fml=l  +  fmlC™‘  (3-10) 

where  is  the  net  aqueous  ED  concentration  in  the  source  zone  and  fmt  is  a  mass  transfer 

enhancement  coefficient.  If  increases  in  source  zone  dissolved  concentrations  are  observed 
during  an  ED  injection  pilot  test,  fnt  may  be  calibrated  by  adjusting  the  model  input  value  until 
the  computed  Fmt  value  agrees  with  observed  mass  transfer  enhancement. 

The  value  ol'C“j  is  computed  as  described  in  Chapter  3  at  coordinates  corresponding  to  the 
center  of  the  DNAPL  source  zone  at  a  date  equal  to  the  injection  start  date  plus  a  duration  t/ag  to 
allow  the  ED  concentration  to  reach  about  99%  of  the  steady-state  value  estimated  as  tulg  =  xR/v 
+  3(2AlxR)i/2/v  where  x  is  the  travel  distance  from  the  ED  injection  gallery  to  the  center  of  the 
DNAPL  source,  R  is  the  ED  retardation  factor,  v  is  the  mean  groundwater  pore  velocity,  and  Al 
is  the  aquifer  longitudinal  dispersivity. 

3.4  Electron  Donor  Injection  Cost  Function 

Capital  and  operating  costs  for  implementation  of  enhanced  bioremediation  by  injection  of  ED 
and  other  amendments  in  injection  well  galleries  are  described  as  follows: 


avail  .  o  H  nat 
ED  KED^ED  - 


CEA+RCHfcH 


(C< 


-C 

CH  ^CH(  mixed) 


)) 


RchIch 
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where  the  time  summation  is  over  integer  values  of  time  in  years, 

Cnpv'  is  ^e  total  NPV  cost  for  ED  injection  ($K), 

C,2pyP  is  the  total  NPV  fixed  ED  cost  ($K), 

ClmP  's  the  fixed  cost  per  ED  gallery  width  ($K/m), 

O'  is  the  construction  cost  per  operational  ED  monitoring  well  ($K/well), 
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CoiherP  is  any  other  fixed  ED  costs, 

C'^py  is  the  total  NPV  operating  cost  for  ED  injection  ($K), 

is  the  operating  cost  per  ED  gallery  width  for  maintenance  etc.  ($K/m), 
is  the  operating  cost  per  unit  ED  mass  injection  ($K/kg), 

C^amp  is  the  collection  and  analysis  cost  per  ED  monitoring  sample  ($K/sample), 

C other  is  other  ED  operating  costs  per  gallery  per  year  for  reporting  etc.  ($K/gallery/yr), 

Can°P  is  other  ED  operating  costs  regardless  of  the  number  of  galleries  ($K/yr), 
d  is  the  annual  discount  rate  (fraction), 

fin,,,  is  the  number  of  samples  per  well  per  year  for  ED  operational  monitoring, 
if  °  is  an  indicator  that  is  1  if  gallery  i  is  actually  implemented  else  0, 

LfD  is  the  width  of  ED  gallery  i  perpendicular  to  the  flow  direction  (m), 

MEDi  is  the  mass  injection  rate  of  ED  for  gallery  i  (kg/yr), 

iV™  is  the  number  of  operational  monitoring  wells  (not  injection  wells)  per  ED  gallery, 

Nf'f  is  the  number  of  potential  ED  galleries, 

OfD  is  an  indicator  that  is  1  if  gallery  i  is  operating  in  year  t  else  0, 

tff  is  the  start  date  for  ED  gallery  i  (yr), 

tref  is  the  reference  date  for  NPV  adjustment  (yr),  and 

tmax  is  the  maximum  simulation  date  (yr). 

3.5  Electron  Donor  Injection  Design  Variables 

The  current  version  of  SCOToolkit  allows  multiple  ED  injection  galleries  to  be  considered.  Each 
ED  gallery  can  have  different  gallery  coordinates,  width,  start  date,  and  a  constant  (average)  ED 
injection  rate.  Termination  of  ED  injection  may  be  specified  as  a  fixed  date  or  the  end  date  can 
be  based  on  operational  monitoring  of  contaminant  concentrations  ( Ced )  in  a  nearby  well.  For 
ED  galleries  downgradient  of  a  DNAPL  source  intended  to  enhance  dissolved  plume 
remediation,  ED  injection  is  terminated  when  the  annually-averaged  measured  contaminant 
concentration  in  a  designated  well  location  (or  locations)  upgradient  of  the  ED  gallery  is  less 
than  a  specified  value  Ced  (provided  any  other  upgradient  ED  galleries  have  terminated  injection 
earlier  than  the  travel  time  between  the  galleries  prior  to  present).  For  ED  galleries  placed 
immediately  upgradient  of  a  DNAPL  source  to  enhance  source  zone  remediation,  an  operational 
monitoring  well  is  designated  downgradient  of  the  source  to  monitor  total  solvents  and  decay 
products  as  an  indicator  of  source  zone  DNAPL  dissolution  rate.  When  the  total  concentration  is 
less  than  a  stipulated  Ced  value,  ED  injection  is  terminated.  Lognormal  noise  is  applied  to 
simulated  annually  averaged  “measurements”  with  a  ln-standard  error  equal  to  Si„ced  1/2 

where  Si„ced  is  the  ln-error  for  single  measurements  and  f  ™  is  the  number  of  measurements  per 
year  (number  of  ED  monitoring  wells  per  gallery  times  annual  sampling  frequency). 

Design  variables  that  may  warrant  optimization  for  ED  injection  galleries  include  the  following: 

Medi  is  the  ED  mass  injection  rate  for  gallery  i  (m), 
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LEDi  is  the  width  of  gallery  i  (m), 

tEDoi  is  the  start  date  for  ED  injection  in  gallery  i, 

Cedi  is  the  contaminant  concentration  at  designated  operational  monitoring  location(s)  below 
which  ED  injection  will  be  terminated,  and 

ffamp  's  the  number  of  measurements  per  gallery  per  year  for  ED  operational  monitoring. 

ED  injection  gallery  coordinates  may  be  optimized  directly  or  several  potential  locations  can  be 
pre-selected  and  optimization  used  to  detennine  which  are  most  effective. 

If  the  value  of  MEdi  is  below  a  specified  lower  cutoff,  the  gallery  is  treated  as  inoperative  and  no 
capital  or  operating  costs  are  applied  to  it. 
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_ Chapter  4 _ 

Thermal  Source  Remediation  Model 


4.1  Overview 

Source  mass  reduction  is  a  key  component  of  many  DNAPL  site  remediation  strategies.  Options 
include  excavation,  soil  vapor  extraction,  chemical  oxidation  and  thermal  treatment.  Our  focus 
here  will  be  on  electrical  resistance  heating  (ERH)  systems,  which  have  broad  applicability  for 
saturated  and  unsaturated  zone  DNAPL  sources.  Unlike  methods  that  depend  on  efficient 
advective  mass  transfer,  ERH  is  relatively  insensitive  to  heterogeneity  in  aquifer  permeability, 
and  is  capable  of  removing  a  high  percentage  of  mass  from  the  treatment  volume.  The  main 
limitation  of  ERH,  which  it  shares  with  most  source  reduction  methods,  is  that  mass  reduction  is 
largely  limited  to  the  treatment  zone  (some  increase  in  biotic  and  abiotic  decay  may  occur 
downgradient  due  to  elevated  temperatures  during  and  for  a  time  after  thermal  treatment).  Hence, 
its  overall  effectiveness  is  highly  dependent  on  the  reliability  of  DNAPL  source  identification 
and  delineation. 

ERH  utilizes  an  array  of  electrodes  installed  within  the  identified  treatment  zone,  which  may  be 
placed  at  discrete  depth  intervals  in  increments  of  the  electrode  length.  Soil  and/or  groundwater 
concentrations  within  the  source  zone  are  generally  determined  prior  to  commencing  treatment  to 
establish  baseline  conditions  and  monitoring  during  treatment  is  performed  to  decide  when  to 
shut  off  electrodes  to  achieve  a  target  reduction  in  soil  and/or  groundwater  concentrations  within 
the  treatment  zone.  In  this  chapter,  we  present  a  cost  estimation  formula,  alternative  operational 
monitoring  methods  to  manage  the  ERH  system,  and  the  modeling  approach  employed  to  assess 
mass  and  flux  reduction  effectiveness,  including  impacts  of  uncertainty  in  source  zone 
delineation. 


4.2  Thermal  Cost  Model 

A  cost  estimation  fonnula  for  ERH  as  a  function  of  treatment  zone  area  and  volume  was 
provided  by  TRS  Group,  Inc.  (Greg  Beyke,  personal  communication)  that  considers  capital  and 
operating  costs  based  on  energy  requirements  computed  to  reach  a  specified  relative  contaminant 
mass  reduction.  We  modified  the  formula  to  explicitly  consider  variations  in  operational 
sampling  costs  and  energy  costs  and  their  dependence  on  real-time  monitoring  for  aquifers  with 
multiple  NAPL  sources  as  follows 
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where 

2 

A  tsri  is  the  areal  extent  of  the  thermal  treatment  zone  (m  ), 

CTRpy  is  the  total  NPV  cost  for  thermal  treatment  for  all  sources  ($K), 

C™  is  a  fixed  cost  for  all  sources  at  a  site  ($K), 
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C™  is  a  cost  multiplier  per  unit  volume  of  the  treatment  zone  to  reach  design  energy  ($K/m3), 
C™„  is  a  cost  multiplier  per  unit  area  of  the  treatment  zone  to  reach  design  energy  ($K/m2), 
C™b  is  the  mobilization  cost  for  each  sampling  event  ($K/event), 

C™u  is  the  installation  cost  per  monitoring  well  ($K/well), 

C Gwsamp  's  the  sampling  and  analysis  cost  per  groundwater  sample  ($K/sample), 

ClRore  is  the  cost  per  soil  boring  ($K/boring), 

CsoiLsamp  's  the  cost  per  soil  sample  analyzed  ($K/sample), 

Efraci  is  the  ratio  of  actual  energy  consumed  when  TSR  tenninates  versus  the  design  estimate 
(model-computed  as  described  in  section  4.4), 
fE  is  the  fraction  of  non-monitoring  variable  costs  attributable  to  energy  use  (-0.22), 

Ft  is  the  total  number  of  soil  and/or  groundwater  samples  divided  by  the  number  in  the  pre¬ 
treatment  sampling  round  (model-computed  as  described  in  section  4.4), 

7™  is  1  if  thermal  treatment  is  performed  for  source  i,  else  0, 

NS0llrce  is  the  number  of  individual  source  zones, 

N™n  is  the  number  of  source  zone  groundwater  monitoring  wells, 

Nsamp/weii  's  the  number  of  sampling  depths  per  well, 

N™re  is  the  number  of  soil  boring  locations  for  each  sampling  time, 

^ HLp: bore  's  the  number  of  depth  intervals  sampled  per  boring, 

RTSR(soii)i  is  target  ratio  of  source  mass  remaining  following  thermal  treatment  to  the  mass 
immediately  prior  to  treatment, 
trsR  is  the  date  that  thennal  treatment  is  performed  (yr), 
tref  is  the  basis  date  for  present  value  (yr), 

Xj  is  a  binary  switch  (=  0  or  1)  to  select/deselect  thennal  treatment  for  specific  sources,  and 
ZrsRi  is  the  vertical  extent  of  the  thennal  treatment  zone  (m). 

Typical  values  forC™  ,  C™  and  C™fl  are  $320,000,  $25/m3  and  $  188/m2  (in  2009),  respectively, 

for  a  source  consisting  predominantly  of  TCE  excluding  design  and  oversight.  Unit  costs  will 
vary  with  the  spatial  configuration  of  the  source  region  (e.g.,  increasing  with  the  ratio  of  source 
perimeter  length  to  source  area),  source  composition  (e.g.,  solvent  species,  content  of  oil,  grease 
and  fuel  hydrocarbons  in  mixture),  local  labor  and  electricity  costs  (based  on  $0.08  per  kWh), 
groundwater  velocity,  and  containment  measures  required.  A  remediation  contractor  should  be 
consulted  to  develop  unit  costs  for  site-specific  conditions. 

In  certain  instances,  a  given  source  region  may  be  characterized  by  two  or  more  source  functions 
representing  different  DNAPL  “architectures”  -  for  example,  pools  and  residual  DNAPL  with 
different  mass,  depletion  exponents,  etc.  For  such  co-located  sources,  the  total  source  zone 
thickness  should  be  divided  among  the  various  for  the  co-located  sources  architecture  and  the 
source  area  similarly  divided  among  the  sources.  Well  installation  and  soil  boring  costs  for 
collocated  sources  must  be  computed  in  a  manner  that  does  not  double  count  drilling  through 
shallower  source  zones. 
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4.3  Operational  Methods  for  Real  Time  Monitoring 

The  cost  function  in  the  preceding  section  is  based  on  an  estimate  of  the  energy  requirement  to 
achieve  a  target  contaminant  mass  reduction  considering  the  heat  capacity  of  the  aquifer, 
advective  heat  loss,  latent  heat  of  volatilization  of  water  and  contaminants,  and  planned 
operating  temperature.  In  practice,  it  is  advantageous  to  base  decisions  on  when  to  terminate 
thermal  treatment  on  real-time  data  from  operational  monitoring  during  remediation.  The 
following  monitoring  and  system  operation  strategies  are  considered  here. 

Method  la  -  Prior  to  heating,  N™re  soil  borings  with  N™mp/bore  samples  per  boring  are  collected 

and  average  initial  soil  concentration  is  computed  for  each  source  treatment  zone  i.  After  a 
period  of  time  when  a  fraction  Efrac(init)i  of  the  theoretical  energy  requirement  has  been  applied, 
the  same  number  of  soil  samples  as  in  the  pre-heating  round  are  collected  and  analyzed.  If  the 
average  measured  soil  concentration  is  less  than  Rrsmtsoun  x  the  initial  average,  heating  is 
tenninated  for  all  heating  elements  within  the  source.  Otherwise,  the  sampling  procedure  is 
repeated  at  time  intervals  corresponding  to  a  fractional  energy  use  Efrac  until  the  criteria  are  met. 

Method  lb  -  Same  as  Method  la  except  N™ir  groundwater  monitoring  wells  with  N™mp/well 

sampling  depths  per  well  are  employed  rather  than  soil  samples  with  an  operational  objective  to 
reduce  groundwater  concentration  by  a  factor  RrsR(gw)i- 

Method  lc  -  Same  as  Method  la  and  lb  but  soil  and  groundwater  samples  are  taken  at  each 
sampling  event  and  heating  is  terminated  when  both  the  average  measured  soil  concentration  is 
less  than  RrsR(soii)i  x  the  initial  average  soil  concentration  and  the  average  measured  groundwater 
concentration  is  less  than  RrsR<gw)i  x  the  initial  average  groundwater  concentration. 

Method  2a  -  This  is  similar  to  Method  la  but  individual  electrodes  are  turned  off  when  the 
measured  soil  concentration  for  a  sample  within  the  zone  heated  by  the  electrode  is  less  than 
RrsR(soii)i  x  the  average  pre-heating  concentration.  A  one-to-one  relation  between  electrodes  and 
sampling  regions  is  assumed.  After  an  electrode  is  turned  off,  no  further  sampling  is  performed 
in  this  region. 

Method  2b  -  Same  as  Method  2a  using  groundwater  rather  than  soil  samples  with  a  target 
groundwater  concentration  reduction  factor  RrsR(gw)i. 

Method  2c  -  Same  as  Method  2a  and  2b  but  both  soil  and  groundwater  samples  are  taken  at  each 
sampling  event.  Electrodes  are  shut  off  incrementally  when  both  the  local  measured  soil 
concentration  is  less  than  RrsR(soii)i  x  the  initial  average  soil  concentration  and  the  average 
measured  groundwater  concentration  is  less  than  RrsR(gw)i  x  the  initial  average  groundwater 
concentration. 

4.4  Thermal  Source  Reduction  Model 

The  model  for  source  mass  reduction  versus  time  due  to  mass  transfer  to  the  dissolved  phase  was 
described  in  Chapter  2.  Here  we  focus  on  source  mass  reduction  during  thermal  treatment 
between  time  Esro  /  when  thermal  treatment  begins  for  source  i  to  trsRfi  when  thermal  treatment  is 
completed.  The  source  mass  at  the  commencement  of  thermal  treatment  ( Mtsr0  <)  is  computed 
from  (2.2)  at  t  =  trsRo  i- 


4-3 


To  consider  the  effects  of  imperfect  source  delineation  and  uncertainty  in  the  area  requiring 
thennal  treatment,  we  divide  Mtsroi  into  (1)  mass  within  the  treatment  zone  {Mtreatment  zone  o  d  and 
(2)  mass  outside  the  treatment  zone  (Mmiss ,),  such  that 


M 


treatment  zone  o  i 


M  TSR  o  i  Mn. 


(4-2) 


Assuming  Mmiss ,  does  not  change  substantially  during  thennal  treatment,  the  total  source  mass 
remaining  after  thermal  treatment  is 


^  TSR  fi  M treatment  zone  fi  ^  miss  i 
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where  Mtreatmentzonefi  is  the  treatment  zone  mass  following  source  heating.  The  “missing”  mass  is 
described  by 


AT 
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where  Fmiss  ,  is  the  fraction  of  the  pre-remediation  source  mass  that  is  outside  the  thermal 
treatment  zone  for  source  i,  i.e.,  Fmissi  >  0  indicates  that  a  fraction  Fmissi  of  the  estimated  source 
mass  is  not  within  the  treatment  volume  and  hence  will  not  be  treated  by  source  heating. 


We  treat  Fmiss ,  as  a  stochastic  variable  controlled  by  uncertainty  in  the  actual  source  volume 
approximated  as 
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where  Vtsr  i  is  the  soil  volume  treated  by  the  thermal  system  for  source  i  and  Vsource ,  is  the 
actual  volume  of  contaminated  soil  for  source  i.  The  latter  is  taken  as  a  stochastic  variable  by 
virtue  of  uncertainty  in  source  zone  dimensions  based  on  field  investigations  and/or  model 
calibration.  The  “actual”  source  volume  is  computed  in  MC  simulations  as  the  product  of  source 
dimensions  Lx,  Ly  and  Lz  generated  from  their  respective  best  estimates  and  covariances  (or  log- 
transforms  if  uncertainty  is  regarded  as  log-normally  distributed). 

Source  mass  within  the  ERH  treatment  zone  is  assumed  to  decrease  with  cumulative  applied 
energy  as 
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where  M treatment  zone  i(t)  is  the  mass  within  the  treatment  zone  at  time  t  during  treatment  at  which  a 
fraction  (or  multiple)  Efrac  j(t)  of  the  estimated  total  energy  has  been  applied.  The  treatment 
duration,  AAsr,  in  years  is  approximated  as 

/  X  -0.087 

^TSRi  =^-^^fracf\^TSR(soil)i)  (4"7) 


with  a  In  standard  deviation  of  -0.25.  Assuming  that  Efmci(i)  =  (t-t tsroi)/ ^trsRi  then 
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with  A trsRi  given  by  (4-7)  with  noise.  From  (2.1),  the  total  mass  discharge  rate  from  the  source 
during  thennal  treatment  is  described  by 
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while  the  discharge  rate  from  the  treatment  zone  alone  is  only 
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Average  treatment  zone  dissolved  concentrations  can  be  inferred  from  (4-10)  assuming  aqueous 
concentration  is  approximately  proportional  to  Jtreatement zone  i  and  the  darcy  velocity,  indicating 
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where  /3cai  i  is  the  best  estimate  of  the  source  depletion  exponent  for  source  i.  Note  that  for  MC 
analyses,  although  ft  is  generally  a  stochastic  variable,  RrsR(gw )  and  Rtsr(soU)  are  design  variables, 
which  must  be  deterministic.  Hence,  conversion  between  soil  and  groundwater  cleanup  levels 
must  be  made  using  the  best  estimate  of  beta  that  is  available  to  the  user. 


Note  that  average  soil  and/or  groundwater  concentrations  may  be  less  than  their  target  values  due 
to  overshoot  at  discrete  sampling  intervals  (including  use  of  a  value  for  Efrac(inu)i  that  overshoots 
the  target  on  the  first  operational  sampling).  Furthermore,  for  Methods  la-lc,  positive  or 
negative  deviations  can  occur  due  to  biases  in  average  measured  concentrations  before  and 
during  thennal  treatment  operation  and  to  biases  in  average  concentrations  before  treatment.  For 
Methods  2a-2c,  positive  or  negative  deviations  in  individual  concentrations  associated  with  each 
electrode  volume  may  occur.  Also,  although  Methods  2a-2c  will  generally  result  in  reductions  in 
energy  and  sampling  costs,  they  will  exhibit  greater  variability  in  actual  cleanup  levels  because 
individual  samples  will  exhibit  greater  variance  than  averaged  values. 


To  simulate  thermal  treatment  perfonnance  for  the  various  methods,  we  compute  average 
treatment  zone  soil  and/or  groundwater  concentrations  at  the  beginning  of  thermal  treatment  and 
at  specified  cumulative  energy  fractions  from  (4-2)  -(4-6)  and  (4-10)  and  apply  log-normal 
“noise”  to  represent  individual  measurement  variability  within  the  source  zone.  Uncertainty  in 
averaged  soil  or  dissolved  concentrations  are  taken  as  the  uncertainty  in  individual  values  (from 
calibration  results  or  statistical  analysis  of  field  data)  times  AT  where  N  is  the  number  of 
samples  averaged.  When  simulated  measurements  indicate  that  remediation  criteria  have  been 
met  for  Methods  la-lc,  the  “actual”  treatment  zone  mass  remaining  is  taken  as  the  noise-free 
treatment  zone  mass.  For  Methods  2a-2c  in  which  electrodes  are  turned  off  incrementally,  the 
mass  in  each  “cell”  within  the  treatment  zone  is  simulated  and  the  total  source  zone  remaining 
when  all  electrodes  have  been  shut  down  is  computed  by  summing  the  final  mass  in  each  cell. 
Note  that  reduced  energy  and  monitoring  costs  associated  with  incremental  electrode  shutdown 
may  be  offset  by  the  greater  uncertainty  in  individual  measurements  than  site  averages,  which 
may  necessitate  operating  longer  to  meet  more  stringent  nominal  cleanup  levels  that  provide  a 
safety  factor  against  performance  uncertainty. 
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4.5  Thermal  System  Design  Variables 

Design  variables  that  may  warrant  optimization  include  the  following  for  each  identified  DNAPL 
source: 

Atsr  is  the  total  soil  area  treated  by  the  thermal  system. 

Operational  monitoring  method  used  to  shut  off  TSR  (i.e.,  Method  1  or  2), 

Rtsr(soU)  is  the  ratio  of  final  soil  concentration  in  the  source  zone  to  the  initial  concentration, 
Efrac(init)  is  the  cumulative  fraction  of  estimated  total  energy  requirement  at  which  the  first 
samples  are  taken  (after  pre-remediation  baseline), 

A Efrac  is  the  incremental  fraction  of  estimated  total  energy  for  subsequent  sampling  events, 
N™u  is  the  number  of  source  zone  groundwater  monitoring  wells, 

Kip/well,  is  the  number  of  sampling  depths  per  well  (if  TV™,  >  0), 

N™.e  is  the  number  of  soil  boring  locations  for  each  sampling  time,  and 
N™mp/bore  's  lhe  number  of  depth  intervals  sampled  per  boring  (if  N™re  >  0). 

Note  that  if  groundwater  measurements  are  used  to  decide  when  to  terminate  thermal  treatment, 
then  the  ratio  of  final  groundwater  concentration  in  the  source  zone  to  the  initial  groundwater 
concentration  ( RrsRfgw >)  is  computed  internally  from  (4-11).  The  value  of  RrsR(gw)  corresponding 
to  a  specified  Rtsr(soU)  is  printed  by  the  code  for  field  use  in  making  system  termination  decisions 
from  source  zone  groundwater  monitoring  data. 
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_ Chapter  5 _ 

Model  Calibration  and  Uncertainty  Analysis 


5.1  Overview 

Before  a  model  can  be  used  to  evaluate  effects  of  various  design  variables  on  the  expected 
performance  and  cost  of  remediation  strategies  of  interest,  a  number  of  site-specific  parameters 
must  be  estimated  from  data  available  from  site  characterization  investigations.  Most  model 
parameters  will  be  subject  to  more  or  less  uncertainty,  depending  on  the  quality  and  quantity  of 
data  available.  In  addition  to  parametric  uncertainty,  models  are,  by  definition,  simplifications  of 
reality,  and  thus  are  subject  to  intrinsic  model  fonnulation  errors  associated  with  explicit  or 
implicit  simplifying  assumptions  that  have  been  invoked.  Measurements  used  to  estimate  model 
parameters  are  themselves  subject  to  sampling  and  analytical  errors,  which  also  contribute  to 
prediction  uncertainty. 

Therefore,  the  first  step  that  must  be  undertaken  prior  to  employing  a  model  to  evaluate  design 
options  is  to  perform  model  calibration  using  all  available  data  and  then  to  quantify  parametric, 
intrinsic  and  measurement  uncertainty.  The  model  can  then  be  used  to  make  forward  estimates  of 
performance  and  cost  and  to  quantify  the  prediction  uncertainty. 

5.2  Model  Calibration  and  Error  Analysis 

The  purpose  of  the  calibration/uncertainty  analysis  module  is  to  detennine  best  estimates  of  key 
model  parameters  and  to  quantify  parametric  and  residual  model  uncertainty.  The  model 
calibration  process  utilizes  field  measurements  (which  may  comprise  various  types  of  data,  e.g., 
contaminant  concentrations,  ED  concentrations,  and  contaminant  fluxes  at  various  locations  and 
times),  and  prior  information  about  parameter  values  and  their  uncertainty.  Assuming  Gaussian 
(or  log-Gaussian)  measurement  errors  and  prior  parameter  distributions,  we  seek  to  minimize  the 
negative  log  of  the  posterior  distribution,  L,  described  by 

E  =  I  wrw  (y  -  h(s))r  R  (9)  1  (y  -  h(s))  + 1  (s  -  s*  f  Q  ‘  (s  -  s* )  (5- 1 ) 

where  y  is  a  vector  of  field  measurements,  s  is  a  vector  of  parameter  values,  s*  is  a  vector  of 
prior  parameter  estimates,  h(s)  is  a  vector  of  model  predictions  corresponding  to  the  field 
measurements,  R  is  a  matrix  of  measurement  covariances  corresponding  to  the  vector  of  data 
types  0  (e.g.,  measured  contaminant  concentration,  source  mass,  source  mass  discharge  rate, 
etc.),  w  is  a  user-defined  weighting  matrix,  and  Q  is  the  covariance  matrix  of  prior  parameter 
estimates.  Off-diagonal  terms  are  disregarded  for  w,  R  and  Q. 

Each  model  parameter  and  each  data  type  may  be  log-transformed  prior  to  application  of  (5-1). 
For  parameters  or  data  types  that  are  physically  constrained  to  be  non-negative  and  that  are 
expected  to  exhibit  a  coefficient  of  uncertainty  greater  than  -20%,  log-transformation  is 
advisable.  Also,  calibration  data  types  that  exhibit  ranges  in  the  data  set  that  extend  over  several 
orders-of-magnitude,  log-transformation  may  be  desirable  if  comparable  relative  error  (as 
opposed  to  absolute  error)  is  desired  over  the  measurement  range.  Otherwise,  the  regression 
results  will  likely  be  controlled  by  absolute  errors  from  a  small  number  of  large  data  values. 
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The  magnitude  of  each  data  type’s  uncertainty  (i.e.,  diagonal  terms  in  R)  is  generally  not  known 
a  priori,  but  a  posterior  estimate  can  be  made  using  the  Restricted  Maximum  Likelihood  (RML) 
algorithm  (Kitanidis,  1987).  Note  that  the  final  estimate  of  residual  prediction  uncertainty  R  for 
each  data  type  represents  the  portion  of  data  variability  that  cannot  be  accounted  for  by  the 
model,  which  may  be  due  to  sampling  or  measurement  errors  and/or  to  intrinsic  limitations  of  the 
model  to  represent  all  processes  in  the  field.  For  simplicity,  we  will  refer  to  this  uncertainty  as 
the  “residual”  error.  A  gradient-based  nonlinear  optimization  algorithm  is  used  to  find  the 
solution  that  minimizes  (5-1). 

A  linearized  approximation  of  the  parameter  posterior  covariance  matrix  is  computed  from  the 
final  results  as 


cov(s)  « (HrR(0)_1  H  +  Q  1 )  ‘  (5-2) 

where  Hij=0hi/0sj  is  a  sensitivity  matrix.  Incorporating  prior  estimates  of  parameters  and  their 

uncertainty  into  the  regression  objective  function  greatly  reduces  nonuniqueness  problems  in  the 
inverse  solution  and  allows  many  more  parameters  to  be  calibrated  than  would  be  possible  with 
unconstrained  optimization.  This  not  only  allows  refinement  of  parameters  with  relatively  low 
uncertainty  that  may  otherwise  be  assumed  at  their  prior  estimates,  but  allows  interactions  among 
more  parameters,  through  the  covariance  matrix,  to  be  taken  into  consideration  in  the  error 
analysis.  In  addition  to  variable  constraints  on  parameters  associated  with  the  stipulation  of  prior 
parameter  uncertainty,  absolute  upper  and  lower  constraints  may  also  be  placed  on  any 
parameters. 

5.3  Monte  Carlo  Model 

Uncertainty  in  forward  simulations  of  remediation  performance  and  cost  is  characterized  using  a 
MC  modeling  approach.  Liu  et  al.  (2010)  have  shown  that  linearized  uncertainty  methods 
compare  well  with  more  rigorous  and  much  more  computationally  intensive  than  Markov  Chain 
Monte  Carlo  (MCMC)  methods  when  data  is  not  inordinately  noisy  and  reasonable  prior 
information  is  available  to  condition  parameter  estimates.  Thus,  we  utilize  linearized  uncertainty 
propagation  methods  to  generate  conditional  parameter  realizations. 

The  first  step  in  performing  MC  simulations  is  to  generate  Nmc  equiprobable  realizations  of 
calibrated  model  parameters  for  the  problem  under  consideration  using  standard  methods  for 
multivariate  Gaussian  distributions  (e.g.,  Press  et  al,  2007)  based  on  parameter  best  estimates  and 
covariances  determined  from  the  calibration.  If  log-transformations  of  parameters  were  used 
during  calibration,  the  same  transformations  are  used  for  parameter  generation.  In  addition, 
uncorrelated  normal  or  log-normal  “noise”  may  be  optionally  applied  to  additional  non- 
calibrated  model  inputs,  including  cost  model  coefficients. 

The  performance  model  is  used  to  simulate  soil  and/or  groundwater  concentrations  in  source 
zones  during  thermal  treatment,  groundwater  concentrations  at  ED  injection  gallery  monitoring 
locations,  and  groundwater  concentrations  at  compliance  wells,  which  are  used  for  making  “real¬ 
time”  operational  decisions  during  the  simulations  (e.g.,  to  turn  remediation  systems  off, 
terminate  monitoring,  or  incur  “penalty  costs”).  Log-normal  noise  is  added  to  these  simulated 
concentrations  for  each  parameter  set  realization  using  a  user-specified  “measurement  error.”  For 
data  types  employed  for  model  calibration,  these  estimates  are  obtained  from  the  residual 
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calibration  error.  Otherwise,  they  are  based  on  independent  statistical  analyses  or  the  user’s 
experience.  NPV  cost  is  computed  for  each  MC  realization  and  the  average  NPV  cost  for  all 
realizations  is  computed  as  an  estimate  of  the  expected  NPV  cost  value. 
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_ Chapter  6 _ 

Remediation  Design  Optimization 


6.1  Overview 

Predictions  of  remediation  perfonnance,  no  matter  how  sophisticated  the  performance  model  that 
is  employed,  are  always  uncertain  to  a  significant  degree.  Quantification  of  performance 
uncertainty  is  critical  for  remediation  design,  since  system  “failure”  nearly  always  occurs  as  a 
result  of  deviations  from  best  estimates  of  performance.  Therefore,  remediation  design  (indeed, 
engineering  design  in  general)  is  an  exercise  in  detennining  the  optimal  degree  of  overdesign  to 
achieve  a  “factor  of  safety”  that  balances  the  cost  of  overdesigning  against  the  cost  of  failing  to 
meet  design  objectives.  Our  approach  to  this  tradeoff  problem  is  to  seek  a  design  that  minimizes 
the  “expected”  cost  across  the  spectrum  of  equi-probable  outcomes  considering  all  sources  of 
prediction  uncertainty.  More  specifically,  we  seek  to  minimize  the  average  NPV  cost  to  reach 
compliance  criteria  over  a  finite  number  of  MC  realizations. 

6.2  Compliance  Rules 

The  overall  goal  of  groundwater  remediation  is  to  reduce  contaminant  concentrations  in 
compliance  wells  below  a  risk-based  or  regulatory-mandated  level  within  a  certain  timeframe  at 
minimum  cost.  For  a  given  remediation  design,  one  of  three  possible  outcomes  is  assumed  to 
apply  for  a  given  remedial  action  at  the  actual  site  and  hence  for  each  MC  model  realization: 

No  Further  Action  (NFA).  If  compliance  well  monitoring  data  meet  specified  NFA  criteria  at  a 
date  after  tn/a  and  prior  to  the  maximum  simulation  date,  tmax,  then  all  remediation  and  monitoring 
activities  will  be  tenninated  (“no  further  action”). 

Non-Compliance  (NC).  If  contaminant  concentrations  for  one  or  more  compliance  wells  exceed 
certain  noncompliance  criteria  after  a  specified  “penalty  date”,  tpenaity,  or  if  NFA  has  not  been 
achieved  prior  to  tmax,  then  a  specified  fixed  present  value  “penalty  cost”  will  be  added  to  the  cost 
function  and  the  simulation  is  tenninated.  No  penalty  cost  will  be  considered  if  tpenaity  >  tmax- 

Implement  Conditional  Containment  (PT).  If  contaminant  concentrations  at  a  designated  PT 
trigger  well  exceed  certain  criteria  after  a  specified  trigger  date  tPT,  then  a  pump-and-treat  (PT)  or 
other  containment  system  is  implemented  upgradient  of  the  trigger  location.  Discounted  capital 
and  annual  operating  costs  for  the  containment  system  are  accrued  until  NFA  or  NC  conditions 
occur  or  tmax  is  reached.  Other  remedial  actions  and  their  costs  continue  to  be  simulated  until  they 
meet  their  termination  criteria  or  tmax  is  reached.  The  PT  system  is  not  explicitly  simulated,  but  is 
assumed  to  contain  the  plume.  Model-predicted  concentrations  at  compliance  wells  are  computed 
as  though  unaffected  by  the  containment  system  -  i.e.,  attainment  of  NFA  is  conservatively 
assumed  to  be  unaffected  by  the  containment  system.  No  PT  costs  are  considered  if  tpr  >  tmax 

We  assume  site-wide  monitoring  of  contaminant  concentrations  at  a  frequency  of  f^np  per  year 

at  each  designated  compliance  well  location.  Since  concentration  measurements  at  a  given  well 
can  exhibit  considerable  temporal  variability,  compliance  rules  must  be  carefully  defined  to  filter 
out  measurement  “noise”  so  that  the  likelihood  of  a  false  NFA  detennination  is  small. 
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Two  compliance  rule  options  to  filter  measurement  noise  are  considered.  The  first  option 
employs  annual  averaging  of  f^np  compliance  well  concentrations  per  year  and  requires  all 

annual  averages  within  a  lookback  period  of  Ni00kback  years  to  be  below  a  specified  target 
concentration  Cmax  to  attain  NFA.  This  approach  is  very  conservative,  but  noisy  data  may 
inordinately  extend  the  remediation  duration. 

The  second  option  which  has  been  recommended  by  the  US  EPA  (Levine,  2010)  involves  a 
linear  regression  of  concentration  versus  time  for  all  measurements  within  a  lookback  period  and 
computes  confidence  limits  of  the  mean  regression  value  at  the  end  of  the  regression  period.  If 
the  upper  one-sided  a-probability  confidence  limit  of  the  regression  is  less  than  Cmax  for  each 
compliance  well,  then  NFA  is  attained.  This  approach  effectively  handles  data  noise  in  the 
decision  process,  but  is  less  sensitive  to  outliers.  A  complete  summary  of  the  compliance  rules  is 
given  in  Table  6-1. 


Table  6-1.  Compliance  rule  protocol  options. 


Extreme  Value  (EXV)  Rule 


Regression  Confidence  Limit  (RCL)  Rule 


NFA  If  annually  averaged  contaminant 

concentrations  for  all  compliance  wells 
are  less  than  Cmax  for  each  of  the  last 
Ni00kback  years  ending  on  or  after  t„fa,  then 
all  remediation  and  monitoring  activities 
are  terminated. 

NC  If  annually  averaged  contaminant 

concentrations  for  any  compliance  wells 
exceed  Cmax  in  any  of  the  last  Ni00kback 
years  ending  on  or  after  tpenaity  or  if  NFA 
has  not  been  achieved  prior  to  tmax,  then  a 
specified  present  value  penalty  cost  Dpenaity 
is  added  to  the  cost  function  and  the 
simulation  is  terminated. 

PT  If  annually  averaged  contaminant 

concentrations  in  a  designated  PT  trigger 
well  exceed  Cpr  for  all  of  the  last  N/00kback 
years  ending  on  or  after  tpp,  then 
discounted  capital  and  operating  costs  for 
pump-and-treat  or  other  plume 
containment  system  are  accrued.  The 
simulation  continues  until  NFA  is 
achieved  or  tmax  is  reached. 


If  the  upper  confidence  limit  of  the  current 
value  of  a  regression  of  contaminant 
concentration  vs.  time  over  the  last  Ni00kback 
years  ending  on  or  after  t„fa  is  less  than  Cmax 
for  all  compliance  wells,  then  all 
remediation  and  monitoring  activities  are 
terminated. 

If  the  upper  confidence  limit  of  the  current 
value  of  a  regression  of  contaminant 
concentration  vs.  time  over  the  last  Ni00kback 
years  ending  on  or  after  tpenaity  exceeds  Cmax 
for  any  compliance  wells  or  if  NFA  has  not 
been  achieved  prior  to  tmax,  then  a  present 
value  penalty  cost  Dpenaity  is  added  to  the 
cost  function  and  the  simulation  is 
terminated. 

If  the  lower  confidence  limit  of  the  current 
value  of  a  regression  of  contaminant 
concentration  vs.  time  over  the  last  Ni00kback 
years  ending  on  or  after  tPT  exceeds  Cpr  for 
a  designated  PT  trigger  well,  then 
discounted  capital  and  operating  costs  for 
pump-and-treat  or  other  plume  containment 
system  are  accrued.  The  simulation 
continues  until  NFA  is  achieved  or  tmax  is 
reached. 


1  Based  on  Levine  (2010) 
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The  earliest  possible  compliance  date,  tnfa,  in  Table  6-1  is  computed  as 

Kfa  =  max,y  (*577}  +  Atij  )  +  ^ lookback  (6‘  1 ) 

where  tsrpi  is  the  date  that  ED  injection  stops  in  gallery  i,  Aty  is  the  travel  time  from  ED  injection 
gallery  i  to  the  compliance  well  j,  and  max,y  is  the  maximum  over  all  ij  pairs  with  gallery  i 
upgradient  of  compliance  well  j.  The  compliance  date  is  used  to  ensure  that  low  compliance  well 
concentrations  due  to  upgradient  ED  injection  are  not  misinterpreted  as  permanently  clean. 

Annual  average  concentrations  in  the  EXV  rule  are  treated  as  geometric  averages.  Log-normal 
measurement  “noise”  is  applied  to  annual  averages  in  MC  realizations  with  a  standard  error 
equal  to  the  In  C  standard  deviation  from  model  calibration  divided  by  f™pm  .  Note  that  the 

effect  of  increasing  sampling  frequency  is  to  decrease  measurement  noise,  which  will  allow 
earlier  NFA  attainment,  and  hence  lower  NPV  cost.  The  optimal  sampling  frequency  will  depend 
on  the  tradeoff  between  NPV  cost  reduction  and  the  additional  sampling  and  analytical  costs. 
The  same  tradeoff  occurs  with  the  RCL  rule,  since  increasing  sampling  frequency  decreases  the 
regression  uncertainty  by/™ip12  . 

The  “penalty  cost”  may  be  a  real  cost  (e.g.,  for  last  ditch  containment  measures)  or  a  fictitious 
cost  applied  to  constrain  the  failure  probability.  Penalty  cost  is  specified  in  NPV  dollars  (i.e.,  no 
discount  is  applied  internally).  Care  should  be  taken  not  to  specify  a  penalty  date  ( tpenaity )  that 
cannot  be  realistically  achieved  with  current  site  conditions  and  proposed  remedial  technologies 
lest  compliance  will  be  nearly  impossible  to  meet  and  the  optimization  problem  will  be  ill- 
defined.  Likewise,  if  tmax  is  too  early  to  achieve  a  high  NFA  probability,  low  remediation  costs 
may  be  misleading.  An  exception  may  be  if  long-term  containment  with  institutional  controls  is 
under  consideration,  tmax  may  be  set  to  reach  a  pseudo-steady  state  condition,  fixed  and  penalty 
costs  set  to  zero,  and  design  variables  optimized  to  minimize  operating  costs. 

Note  that  pump-and-treat,  if  considered,  is  not  explicitly  simulated.  Pump-and-treat  is  assumed 
to  contain  the  plume,  but  model-predicted  concentrations  at  compliance  wells  are  treated  as 
though  unaffected  by  the  PT  system. 

Real-time  implementation  of  Option  2  would  involve  perfonning  regressions  of  the  last  N  = 
fZp  Niookback  ln-concentration  values  after  each  sampling  round  versus  Julian  date,  T,  to  fit  a 
linear  trend  model 


In  C  -  a  +  bT 

where  a  and  b  are  regression  coefficients.  The  regression  error  is 

2  _  SSE 
*  N-2 


(6.2) 


(6.3) 


where  SSE  is  the  sum  of  squared  In  C  regression  residuals.  Uncertainty  in  the  mean  regression  is 

S„ 


S,  =  ■ 


r  NU2 

and  one-sided  confidence  limits  at  probability  a  are  computed  by 


(6.4) 
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LCL„  =exp (a  +  bT-t(a,fZfNk^t  -1)S,) 

,  ;  (6.5) 

UCL„  =  exp  (a  +  bT  +t(a,fZ„  -1)S,) 

where  t(a,  df)  is  the  one-sided  Student-/  value  for  probability  a  with  df  degrees  of  freedom. 

To  reduce  computational  effort  in  the  MC  model,  we  generate  a  single  (geometric)  average 
concentration  value  for  each  year  in  each  realization  with  noise  S/nc/  fssa„p 1 2  rather  than  f  ™np 

sample  values  per  year  with  noise  Si„c,  where  Si„c  is  the  posterior  In  C  calibration  error.  The 
regressions  are  performed  on  time  series  with  N  =  Ni00kback  concentration  values,  which  yield  a 
lower  SSE  due  to  annual  averaging.  Equations  (6.2)  through  (6.5)  are  employed  to  compute 
confidence  limits  using  N  =  Ni00kback ■  Note  that  in  actual  field  practice,  regressions  would  be 
performed  on  N  =  f  ™np  Ni00kback  concentration  values  to  detennine  when  NFA  is  met. 

It  should  be  self-evident  that  the  duration  of  time  allowed  before  cleanup  criteria  are  required  to 
be  met  will  have  a  significant  effect  on  remediation  design  and  cost.  Specifically,  if  tpenaity  or  tpj 
are  decreased  or  the  magnitude  of  contingent  containment  or  penalty  costs  are  increased,  then 
earlier,  more  aggressive  and  more  costly  remediation  will  be  favored.  Likewise,  if  the  time 
discount  factor,  d,  is  decreased,  future  costs  are  less  sharply  discounted,  thus  favoring  earlier  and 
more  aggressive  action.  In  the  case  of  no  contingent  penalty  or  containment  costs  (or  /pe„a/^>  tmax 
and  tpT  >  tmax)  with  a  positive  discount  rate,  the  cost  optimal  solution  will  be  to  simply  monitor 
until  tmajc  or  NFA  is  reached,  regardless  of  the  probability  of  achieving  NFA.  Some  cost 
consequence  of  “failure”  must  be  stipulated  to  induce  active  remediation. 

6.3  Site-Wide  Costs  and  Total  Cost 

The  NPV  cost  for  site-wide  monitoring,  reporting  and  maintenance  are  computed  as  follows 


s~iSWtot  _ 

^NPV  ~ 

s~lSWcap  _  (  S~iSWcap  -KjSW  s^SWcap  \  /I  _  start  ref 

^ NPV  —y-'well  well  ^ other  )  VA  a) 


t  t  start 


and  the  total  NPV  cost  for  all  site  remediation  activities  and  penalty  costs  is 


s-iSWcap  ,  s~iSWop 
^ NPV  ^ NPV 


(6.6) 


Call  _  s~t 

MPT/  ^  . 


■'NPV 


SWtot 

NPV 


Z  0(1 


,  J  /~i  pen 
”r  1  pen^  NPV 


where 


CmpV  is  the  total  NPV  remediation  and  penalty  cost  ($K), 
Cfjp'y  is  the  total  NPV  site -wide  cost  ($K), 

C^p'y  is  the  total  NPV  ED  injection  system  cost  ($K), 
C™y  is  the  total  NPV  thermal  remediation  cost  ($K), 
CZ7  is  the  total  PT  fixed  cost  ($K), 


(6.7) 
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CZai  is  the  total  PT  operating  cost  per  year  ($K/yr), 

C^pV  is  the  NPV  penalty  cost  ($K), 

Cs^yP  is  the  total  NPV  site-wide  fixed  cost  ($K), 

C^ap  is  the  fixed  cost  for  monitoring  well  construction  ($K/well), 

CZ7  is  any  other  site- wide  fixed  costs  ($K), 

is  the  total  NPV  operating  cost  for  site-wide  monitoring  and  reporting  ($K), 

Csl°p  is  the  cost  per  sample  for  site- wide  monitoring  ($K/sample), 

CoI°r  is  other  annual  site-wide  operating  costs  ($K/yr), 

f^up  is  the  number  of  samples  per  well  per  year  taken  for  site-wide  monitoring, 

/  is  1  if  a  penalty  cost  is  triggered,  else  0, 

7pris  1  if  PT  implementation  is  triggered,  else  0, 

N'^i  is  the  number  of  site-wide  monitoring  wells  (including  compliance  wells), 
tPT  the  year  that  PT  is  triggered  (yr), 

tnfa  is  the  year  compliance  is  achieved  or  the  max  simulation  date  for  NC  (yr), 
tstart  is  the  first  year  capital  costs  are  incurred  (yr), 
tref  is  the  basis  date  for  NPV  adjustment  (yr),  and 
d  is  the  annual  discount  rate  (fraction). 

Total  NPV  remediation  cost  with  penalty,  Ca"PV ,  is  the  objective  function  to  be  minimized. 

6.4  Design  Optimization 

The  design  optimization  problem  is  formulated  by  defining  the  specific  NFA,  NC  and 
(optionally)  PT  rules  and  remediation  options  to  be  considered  involving  possible  thermal  source 
remediation,  potential  locations  for  ED  injection  galleries  upgradient  and/or  downgradient  of 
DNAPL  source  zones,  and  options  for  PT  or  alternative  Plan  B  systems  to  be  employed  in  the 
event  of  failure.  Design  variables  to  be  optimized  will  need  to  be  identified  and  unit  costs 
determined  considering  site  specific  factors. 

The  mathematical  objective  for  the  design  optimization  problem  is  to  find  an  optimal  set  of 
design  variables  that  minimizes  total  expected  NPV  (ENPV)  cost.  NFA,  NC  and  PT  rules  serve 
as  implicit  constraints  on  the  optimization  conditioned  through  the  cost  and  performance  models. 
In  addition  to  any  remediation  technology  specific  design  variables,  site -wide  variables  (e.g., 
sampling  frequency,  fsamP,  and  lookback  period,  Ni00kback )  may  also  be  optimized. 

An  efficient  and  robust  genetic  optimization  algorithm,  which  is  resistant  to  local  minima  in  the 
response  function,  is  employed  to  perform  the  optimization. 

Optimization  output  includes  optimal  ENPV  costs  (with  and  without  penalty  costs),  cost 
breakdowns,  and  cost  probability  distributions,  in  addition  to  optimized  design  variables.  Cost 
and  perfonnance  for  individual  MC  realization  as  well  as  for  expected  values  can  be  extracted. 
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_ Chapter  7 _ 

Application  to  a  Hypothetical  Problem 


7.1  Problem  Description 

A  hypothetical  problem  is  considered  involving  a  TCE  plume  in  an  unconfined  aquifer.  A 
DNAPL  release  was  assumed  to  commence  in  1965  resulting  in  a  DNAPL  pool  at  the  bottom  of 
the  aquifer  and  residual  DNAPL  in  the  upper  portion  of  the  aquifer  designated  as  Source  1  and 
Source  2,  respectively. 

Vertically  averaged  dissolved  plume  concentrations  were  computed  using  the  “true”  parameters 
at  26  monitoring  well  locations  between  1985  and  2009.  Simulated  concentrations  (Q  were 
subjected  to  ln-normally  distributed  “noise”  to  represent  deviations  due  to  sampling  and 
measurement  error,  spatial  and  temporal  variability,  and  model  simplifications  with  a  ln-standard 
deviation  (Si„c)  of  0.5.  The  resulting  noisy  data  were  assumed  to  represent  field  observations 
available  for  model  calibration.  The  hypothetical  site  characterized  by  different  parameter  sets 
are  summarized  in  Table  7-1.  Well  locations  and  “measured”  concentrations  in  2009  are 
illustrated  in  Figure  7-1. 


1300 

1200 

E  1100 

|  1000 
z 

900 
800 
700 

900  1000  1100  1200  1300  1400  1500  1600  1700  1800  1900  2000  2100 

East,  m 

Figure  7-1.  TCE  concentrations  in  monitoring  wells  in  2009  and  locations  of  DNAPL  source, 
compliance  well,  ED  injection  galleries,  and  wells  for  monitoring  upgradient  of  ED  galleries. 
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Model  calibration  was  performed  as  described  in  Chapter  5  using  the  “data”  described  above  and 
prior  estimates  of  model  parameters  summarized  in  Table  7-1  based  on  information  assumed  to 
be  obtained  from  site  characterization  and/or  literature  studies.  Final  parameter  estimates  and 
their  estimation  standard  errors  are  summarized  in  Table  7-1  (error  covariances  between  all 
parameters  were  also  computed  but  are  not  shown).  The  final  a  posteriori  estimate  of  residual 
error  in  TCE  concentrations  not  accounted  for  by  the  model  was  S/nc  =  0.56.  Observed  versus 
calibrated  concentrations  are  compared  in  Figure  7-2. 
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Table  7-1.  True  model  parameters,  prior  information,  and  final  estimates. 


Parameter 

PDF1 

True  Value 

Prior  Estimates 
Value  St  Dev3 

Calibrated  Values 
Value  St  Dev3 

Source  1  mass  on  calib.  date  (kg) 

LN 

3646 

2715 

0.60 

2948 

0.52 

Source  1  rate  on  calib.  date2  (kg/d) 

LN 

0.088 

0.120 

0.60 

0.067 

0.21 

Source  1  depletion  exponent  (-) 

N 

0.60 

0.75 

0.25 

0.73 

0.24 

Source  2  mass  on  calib.  date2  (kg) 

LN 

56.4 

45.0 

0.60 

44 

0.48 

Source  2  rate  on  calib.  date2  (kg/d) 

LN 

0.012 

0.010 

0.60 

0.008 

0.37 

Source  2  depletion  exponent  (-) 

N 

1.30 

1.26 

0.25 

1.32 

0.21 

Source  length,  Lx  (m) 

LN 

20.0 

20.0 

0.00 

20.00 

0.00 

Source  width,  Ly  (m) 

LN 

20.0 

23.0 

0.50 

18.3 

0.39 

Release  date  (yr) 

N 

1965.0 

1963.0 

3.00 

1963.9 

2.53 

Plume  darcy  velocity  (m/d) 

LN 

0.070 

0.065 

0.60 

0.049 

0.15 

Aquifer  porosity  (-) 

N 

0.30 

0.35 

0.04 

0.36 

0.04 

Longitudinal  dispersivity  (m) 

LN 

20.0 

17.0 

0.50 

19.8 

0.08 

Transverse  dispersivity  (m) 

LN 

2.00 

4.00 

0.50 

2.1 

0.03 

Coordinate  rotation,  a  (degrees) 

N 

5.00 

5.99 

0.00 

5.99 

0.00 

Aquifer  thickness  (m) 

LN 

30.0 

34.0 

0.10 

34.7 

0.10 

Fl-equiv  EA  concentration  (g/m3) 

LN 

2.75 

2.39 

0.15 

- 

- 

Fl-equiv  O2  concentration  (g/m3) 

LN 

0.57 

0.49 

0.15 

- 

- 

H-equiv  ED  ratio  (-) 

LN 

0.31 

0.26 

0.15 

- 

- 

H-equiv  ratio  CH  (-) 

LN 

0.061 

0.051 

0.15 

- 

- 

Serial  decay  fraction  (-) 

LN 

0.32 

0.50 

0.15 

- 

- 

1  Assumed  probability  distributions:  N  =  normal,  or  LN  =  lognormal. 

2  Source  mass  and  discharge  rate  on  tcaj  =  1990. 

3  Standard  deviations  of  LN  variables  are  log-transformed  (dimensionless);  all  other  values  are  in  specified  units. 


Observation,  ppb 


Figure  7-2.  Observed  versus  calibrated  TCE  concentrations. 
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7.2  Remediation  Options  and  Unit  Costs 

Remediation  options  to  be  considered  in  this  hypothetical  problem  are  source  zone  mass 
reduction  using  thennal  source  removal  by  ERH,  and  dissolved  plume  bioremediation  with 
emulsified  vegetable  oil  (EVO)  injection  as  ED  in  a  single  injection  well  gallery.  ERH  is  to  be 
initiated  in  January  2010.  The  duration  of  ERH  and  timing  and  operation  of  ED  injection 
galleries  will  be  detennined  through  optimization  of  relevant  design  variables.  The  remediation 
goal  is  to  achieve  contaminant  concentrations  in  the  compliance  well  below  Cmax  =  5  pg/L  by 
tpenaity  =  2050  with  NFA  being  met  by  tmax  =  21 10.  If  the  latter  is  not  accomplished,  a  penalty  cost 
will  be  incurred.  Compliance  rules  described  in  Chapter  6,  ED  operation  variables  given  in 
Chapter  3,  and  TSR  variables  in  Chapter  4  will  be  investigated  with  cost  variables  summarized  in 
Table  7-2.  Uncertainty  of  cost  variables  was  not  considered  in  this  hypothetical  problem. 


Table  7-2.  Cost  variables  used  in  hypothetical  problem. 


Definition 

Variable 

Value 

Unit 

Thermal  Treatment  Costs 

Fixed  costs  for  design,  operation  and  reporting 

s~iTR 
^  site 

320 

$K 

Treatment  volume  cost  coefficient 

s~iTR 
^ vo^ 

0.025 

$K/m3 

Treatment  area  cost  coefficient1 

s~iTR 

area • 

0.188 

$K/m2 

Cost  per  source  soil  boring1 

siTR 

K~'borei 

3.5 

$K 

Analysis  cost  per  soil  sample 

s~iTR 

^  SOILs  amp 

0.1 

$K/sample 

Source  monitoring  well  installation  cost  per  well1 

s-iTR 
^ welli 

7.0 

$K 

Analysis  cost  per  groundwater  sample 

ED  Iniection  Costs 

s~iTR 

^ GWsamp 

0.5 

$K/sample 

Fixed  cost  per  meter  of  ED  gallery 

fiEDcap 
^  width 

0.096 

$K/m 

Fixed  cost  per  ED  performance  monitoring  well 

s^iEDcap 

mw 

10 

$K/well 

Other  fixed  cost  for  ED  injection 

s^iEDcap 

other 

65 

$K 

Operating  cost  per  ED  mass  injected 

s~iEDop 

mass 

0.00632 

$K/kg 

Other  operating  cost  per  year  for  ED  injection 

s~iEDop 

other 

240 

$K/year 

Cost  per  ED  monitoring  event 

/~iEDop  J\J  ED 
^ samp  v  mw 

5 

$K/event 

Other  Costs 

2 

Penalty  cost  incurred  for  non-compliance 

/n  pen 
^ NPV 

25,000 

$K 

Cost  per  site -wide  monitoring  event 

s~iSWop  jySW 
samp  v  well 

20 

$K 

Reference  year  for  present  value  cost 

tref 

2010 

year 

Discount  rate  used  for  incurred  costs 

d 

0.03 

- 

1  Since  sources  1  and  2  are  collocated  (i.e.,  pool  and  residual  within  the  same  area),  the  treatment  thickness  and 
indicated  costs  are  allocated  15%  to  Source  1  and  85%  to  Source  2. 

2  Penalty  cost  is  a  fictitious  value  employed  to  limit  noncompliance  probability. 
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7.3  Case  Study  Results 

7.3.1  Optimization  for  different  TSR  operational  monitoring  strategies  (Cases  1  and  2) 

A  set  of  optimization  simulations  was  performed  to  evaluate  effects  of  the  six  TSR  operational 
monitoring  strategies  described  in  Chapter  4  on  performance.  The  simulations  are  designated  as 
follows: 

•  Case  la  -  Optimize  design  using  TSR  Method  1  with  soil  concentration  data, 

•  Case  lb  -  Optimize  design  using  TSR  Method  1  with  dissolved  concentration  data, 

•  Case  lc  -  Optimize  design  using  TSR  Method  1  with  soil  &  dissolved  concentration  data, 

•  Case  2a  -  Optimize  design  using  TSR  Method  2  with  soil  concentration  data, 

•  Case  2b  -  Optimize  design  using  TSR  Method  2  with  dissolved  concentration  data,  and 

•  Case  2c  -  Optimize  design  using  TSR  Method  2  with  soil  &  dissolved  concentration  data. 

A  95%  RCL  criteria  with  a  5-yr  lookback  period  ( Ni00kback )  was  employed  to  define  NFA  and  NC 
for  all  simulations.  The  following  design  variables  were  optimized: 

•  TSR  treatment  area,  A  TSr  (same  for  Source  1  and  2), 

•  TSR  removal  fractions,  Rtsr(soU)i  for  source  1  and  Rtsr(soU)2  for  source  2, 

•  Energy  fractions  at  first  sampling  during  TSR,  Efrac(init)i  and  Efrac(init)2  for  sources  1  and  2, 

•  ED  gallery  injection  start  date,  tEDo, 

•  ED  injection  rate,  Med, 

•  TCE  concentration  in  monitoring  well(s)  upgradient  of  ED  gallery  below  which 
injection  is  stopped,  Ced,  and 

•  Compliance  and  ED  monitoring  frequency,/^  =  f™p  =  f™'p . 

The  following  variables  were  fixed  during  optimization: 

•  TSR  start  date  =  2010, 

•  Energy  increment  for  sampling  after  Efmc(inu )  during  TSR,  AEfmc  =  0.2, 

•  TSR  treatment  thickness,  Ztsr  =  3.6  m  for  Source  1  and  26.4  m  for  Source  2, 

•  Number  of  source  zone  soil  borings  or  groundwater  wells  for  TSR,  N™n  =  N™.(,  =  5, 

•  Number  of  source  zone  soil  samples  per  boring  or  groundwater  samples  per  well  for 
TSR>  Klpi.eu,  =  KLpibore,  =  2  for  Source  1  and  15  for  Source  2,  and 

•  Width  of  ED  injection  gallery,  LEd  =  350  m. 

An  initial  estimate  of  the  ED  injection  rate  was  obtained  based  on  an  electron  balance  to  reduce 
observed  levels  of  EAs  and  contaminant.  Assuming  no  aqueous-solid  partitioning,  the  required 
ED  injection  rate  is  estimated  as 


SF  q  W Lz 


M 


IpAi^FM  +  fcH^i 


max 

CH 


V  i 


ED(est) 


f 


ED 


(7-1) 
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where  q  is  darcy  velocity,  W  is  the  injection  well  gallery  width,  Lz  is  saturated  aquifer  thickness, 
C"))s  is  the  maximum  contaminant  concentration  entering  the  gallery,  SF  is  a  safety  factor,  and 
other  variables  were  defined  previously.  Employing  EA  concentrations  and  stoichiometry  from 
Table  3-1  with  SF  =  2  yields  Med  ~  16  kg/d  for  a  350  m  width  gallery.  During  optimization,  if 
Med  drops  below  MEDmm  =  8  kg/d  (SF=1),  no  ED  injection  is  assumed  to  be  performed  and  all 
ED  costs  are  set  to  zero. 

Measurement  “noise”  (Si„c)  for  thermal  system  monitoring  in  the  source  zone  was  assumed  to  be 
0.25  for  soil  concentrations  and  0.15  for  groundwater  concentrations  for  Source  1  and  0.50  and 
0.30,  respectively,  for  Source  2. 

Results  of  Case  la-2c  optimization  simulations  are  tabulated  in  Table  7-3.  Total  EPNV  costs, 
excluding  penalty  costs,  ranged  from  $1852K  to  $1939K  and  maximum  NPV  costs  for  100  MC 
realizations  ranged  from  $1922K  to  $2409K.  Simulated  median  concentration  at  the  compliance 
location  versus  time  with  95  and  99%  confidence  limits  are  shown  in  Figures  7-3a  to  7-8a  and 
cost  probability  histograms  are  given  in  Figures  7-3b  to  7-8b  for  all  cases.  Note  that 
concentrations  are  truncated  at  an  assumed  detection  limit  corresponding  to  10%  of  the  cleanup 
target. 

Probability  of  noncompliance  was  3%  for  all  cases  except  lb,  which  was  2%.  For  all  cases,  the 
optimization  algorithm  determined  that  ED  injection  was  not  cost  effective.  The  optimum 
sampling  frequency  for  compliance  monitoring  was  determined  to  be  once  per  year  ( fsamp=  1)  for 
all  cases  except  lb,  which  was  semi-annual  (fsamp=  2). 

All  six  optimizations  yielded  design  thermal  treatment  areas  ( Ajsr )  about  2.35  times  the  best 
estimate  of  source  area  (Aca/)  to  reduce  the  probability  of  leaving  untreated  DNAPF  beyond  the 
treatment  zone.  The  large  multiplier  reflects  relatively  large  uncertainty  in  the  calibration  source 
area  and  high  sensitivity  of  total  cost  to  untreated  mass.  Since  the  cost  of  thermal  treatment  could 
be  reduced  on  the  order  of  $500K  by  decreasing  the  thermal  treatment  area  by  half,  additional 
source  characterization  effort  is  probably  warranted,  especially  considering  that  we  really  do  not 
know  exactly  where  the  unaccounted  mass  occurs. 

Optimized  values  of  the  soil  mass  reduction  ratio  (Rtsr(soU))  for  Source  1  (DNAPF  pools)  ranged 
from  0.026  to  0.082  for  Method  1  and  from  0.009  to  0.031  Method  2.  However,  this  variable  is 
effectively  merely  a  minimum  mass  reduction  factor  that  is  overridden  a  high  initial  energy 
fraction  ( Efmc(mit f)  for  the  first  sampling  event  after  system  startup.  The  fact  that  the  average 
number  of  sampling  events  (including  the  baseline  sampling  before  heating)  is  exactly  2  for  all 
cases  for  Source  1  (Table  7-3)  indicates  that  optimization  sets  Efrac(init)  so  large  (~1 .8-3.5  times 
the  estimated  energy  requirement)  that  the  Rtsr(soU)  criteria  is  met  for  all  MC  realizations  at  the 
first  post-startup  sampling  event.  In  other  words,  Efrac(init),  which  is  a  surrogate  for  thermal 
treatment  duration,  becomes  the  primary  design  variable  rather  than  the  stipulated  mass 
reduction  ratio  for  Source  1  (DNAPF  pool). 

This  is  not  as  much  the  case  for  Source  2  (residual  DNAPF),  which  has  a  lower  Efrac(init)  (0.82  to 
2.99)  and  an  average  number  of  source  sampling  events  between  2.01  and  3.5.  The  optimized 
soil  mass  reduction  ratio  ranges  from  0.020  to  0.096.  The  computed  groundwater  reduction  ratio 
(RrsRfgw/)  for  Source  2  is  several  times  lower  than  Rtsr(soU)  reflecting  a  source  depletion  exponent 
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(/?)  greater  than  1,  while  it  is  several  times  higher  than  Rtsr(soU)  for  Source  1  reflecting  a  source 
depletion  exponent  (/?)  less  than  1  (see  (4-11)). 

The  further  the  source  depletion  exponent  deviates  from  1,  the  greater  the  difference  between 
soil  and  groundwater  reduction  ratios.  This  is  important  to  recognize  when  setting  and 
monitoring  source  cleanup  targets.  It  is  notable  that  a  significantly  smaller  (more  aggressive) 
minimum  soil  reduction  ratio  is  needed  for  the  low  /3  source  (Source  1,  pools)  than  for  the  high  f) 
source  (Source  2,  residuals)  to  obtain  a  comparable  reduction  in  dissolved  concentrations.  This  is 
somewhat  offset  by  the  observation  that  the  optimal  minimum  groundwater  reduction  ratio 
( RrsR(gw >)  for  Source  1  (0.032  to  0.103)  is  less  aggressive  (higher)  than  that  for  Source  2  (0.012  to 
0.046),  which  likely  reflects  the  fact  that  the  rate  of  mass  flux  reduction  with  time  ( dJ/dt ) 
increases  with  time  if  ft  <  1  while  it  decreases  with  time  if  J3>  1. 

Comparison  of  Cases  la  through  2c  indicates  that  monitoring  of  soil  concentration  alone  during 
thennal  treatment  yields  the  lowest  ENPV  for  Method  1  (turning  off  all  heating  units  at  the  same 
time  when  the  average  soil  and/or  groundwater  concentration  in  the  source  drops  below  a 
specified  level),  while  soil  and  groundwater  data  yield  the  lowest  cost  for  Method  2  (turning  off 
the  electrodes  incrementally  based  on  local  concentration  reductions).  Method  1  with  soil  data 
yields  the  lowest  cost  of  all  six  cases  studied.  Including  the  probability-weighted  penalty  cost  in 
the  ENPV  leads  to  the  same  conclusion.  Thus,  for  this  particular  case,  Method  1  with  soil  data 
appears  most  cost  effective.  Although  measurement  error  for  soil  concentration  data  is  assumed 
to  be  higher  than  for  groundwater  data,  the  average  cost  per  sample  for  groundwater  data  is 
higher  due  to  well  construction  costs  if  wells  are  only  sampled  2-3  times. 

The  average  number  of  sampling  events  for  monitoring  Source  2  thermal  treatment,  hence  the 
average  thermal  monitoring  cost,  is  significantly  higher  using  Method  2  than  Method  1.  This 
reflects  greater  variability  in  individual  measurements  used  to  incrementally  turn  off  heating 
units  compared  to  averaged  measurements  used  to  turn  off  the  entire  system.  For  the  example 
problem  considered  here,  higher  thermal  monitoring  costs  more  than  offset  cost  savings 
associated  with  incremental  shutdown  of  heating  units.  Due  to  the  large  number  of  factors  and 
tradeoffs  that  affect  the  total  expected  cost  NPV  cost,  it  is  no  more  likely  that  given  sampling 
strategy  will  be  universally  optimal  than  that  any  other  design  variable  will  be  optimal  for  all 
conditions. 

The  difference  between  the  highest  and  lowest  optimized  ENPV  cost  for  the  six  different 
operational  methodologies  for  monitoring  thermal  treatment  is  $426K  without  including  penalty 
costs,  which  corresponds  to  potential  cost  reduction  of  19%  relative  to  the  most  expensive 
thermal  monitoring  strategy.  Thus,  in  addition  to  optimizing  various  treatment  design 
parameters,  evaluating  various  options  for  monitoring  thermal  performance  can  result  in 
significant  additional  cost  savings. 
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Table  7-3.  Summary  of  design  simulations 
(bold  =  optimized  values;  italic  =  fixed  values,  normal  type  =  performance  results). 


Case 

la 

lb 

lc 

2a 

2b 

2c 

3 

4 

5a 

5b 

Results 

ENPVtotal  ($K)' 

1852 

2278 

1932 

1907 

1939 

2025 

11525 

1867 

1305 

1336 

-  TSRtotal  ($K) 

1435 

1450 

1514 

1490 

1521 

1607 

826 

1433 

887 

888 

-  EDtotai  ($K) 

0 

0 

0 

0 

0 

0 

8377 

0 

0 

0 

-  TSRmonitor($K) 

40 

97 

135 

51 

133 

208 

35 

40 

34 

34 

-  Site  wide  ($K) 

417 

828 

418 

417 

418 

418 

2322 

434 

418 

448 

Max  NPV  ($K) 

1922 

2409 

2037 

1968 

2008 

2100 

13036 

1908 

1379 

1379 

Non-compliance  (%) 

3 

2 

3 

3 

3 

3 

71 

8 

3 

53 

Comnliance  and  ED  Monitoring 

Compliance 

RCL 

RCL 

RCL 

RCL 

RCL 

RCL 

RCL 

EXV 

RCL 

RCL 

fsamp  (yr  ) 

1 

2 

1 

1 

1 

1 

4 

1 

1 

1 

Thermal  Treatment  Parameters  (Both  Sources) 

TSR  Method 

1 

1 

1 

2 

2 

2 

1 

1 

1 

1 

TSR  Data 

S 

C 

s+c 

S 

C 

S+C 

S 

S 

S 

S 

Arsn/Aeal 

2.37 

2.34 

2.35 

2.35 

2.34 

2.35 

1.00 

2.12 

1.00 

1.00 

Thermal  Treatment  Parameters  (Source  1  - 

Pool) 

RtSR(soU) 

0.0444 

0.0260 

0.0824 

0.0090 

0.0308 

0.0217 

0.0100 

0.0266 

0.0320 

0.0320 

R-TSRfgw) 

0.1027 

0.0693 

0.1612 

0.0320 

0.0786 

0.0608 

0.0345 

0.0706 

0.0807 

0.0807 

Efracfinit) 

3.23 

1.82 

3.52 

2.18 

2.27 

3.23 

0.80 

3.87 

2.31 

2.31 

Avg  sample  events 

2.00 

2.00 

2.00 

2.00 

2.00 

2.00 

3.48 

2.00 

2.00 

2.00 

Thermal  Treatment  Parameters  (Source  2  - 

Residual) 

RtSR(soU) 

0.0712 

0.0414 

0.0963 

0.0199 

0.0351 

0.0899 

0.0100 

0.0327 

0.0860 

0.0860 

RrsR(gw) 

0.0306 

0.0150 

0.0456 

0.0057 

0.0120 

0.0416 

0.0023 

0.0110 

0.0392 

0.0392 

Efracfinit ) 

1.08 

1.03 

1.09 

2.99 

0.84 

0.82 

0.80 

1.34 

1.07 

1.07 

Avg  sample  events 

2.01 

2.08 

2.01 

2.64 

3.3 

3.54 

3.55 

2 

2.11 

2.11 

ED  Injection  for  Plume  Bioremediation 

tEDo  (yr) 

- 

- 

- 

- 

- 

- 

2010.0 

- 

- 

- 

Med  (kg/d) 

- 

- 

- 

- 

- 

- 

16 

- 

- 

- 

CED  (gg/L) 

- 

- 

- 

- 

- 

- 

5 

- 

- 

- 

Med  avg  (yr) 

- 

- 

- 

- 

- 

- 

71.1 

- 

- 

- 

1  NPV  costs  do  not  include  penalty  cost 
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Figure  7-3.  Results  for  Case  la:  (a)  TCE  concentrations  at  compliance  well 
and  (b)  NPV  cost  distribution  (without  penalty  cost). 


Figure  7-4.  Results  for  Case  lb:  (a)  TCE  concentrations  at  compliance  well 
and  (b)  NPV  cost  distribution  (without  penalty  cost). 


Year  NPV,  $M 

Figure  7-5.  Results  for  Case  lc:  (a)  TCE  concentrations  at  compliance  well 
and  (b)  NPV  cost  distribution  (without  penalty  cost). 
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Concentration,  ppb  Concentration,  ppb  Concentration,  ppb 


Year  NPV,  $M 

Figure  7-6.  Results  for  Case  2a:  (a)  TCE  concentrations  at  compliance  well 
and  (b)  NPV  cost  distribution  (without  penalty  cost). 


Figure  7-7.  Results  for  Case  2b:  (a)  TCE  concentrations  at  compliance  well 
and  (b)  NPV  cost  distribution  (without  penalty  cost). 


Year  NPV,  $M 

Figure  7-8.  Results  for  Case  2c:  (a)  TCE  concentrations  at  compliance  well 
and  (b)  NPV  cost  distribution  (without  penalty  cost). 
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7.3.2  Improvement  attributable  to  optimization  (Case  3) 

A  non-optimized  MC  analysis  of  remediation  performance  and  cost  was  performed 
corresponding  to  Case  la  using  design  variables  that  reflect  common  field  engineering  practices 
(Case  3).  Thermal  source  zone  remediation  was  assumed  to  start  in  2010.  Thermal  treatment  area 
and  volume  were  taken  as  their  best  estimates  from  model  calibration  (i.e.,  Atsr=Acui).  The  first 
round  of  soil  sampling  after  starting  thermal  treatment  is  assumed  to  occur  at  a  fraction  Efrac(init)  = 
0.8  of  the  theoretical  energy  requirement  for  each  source  with  additional  sampling  rounds  at 
intervals  A Efrac  =  0.2  until  shutdown  according  to  TSR  Method  1.  Thermal  remediation  is 
terminated  when  the  average  measured  soil  concentration  in  the  source  zone  is  below  Rtsr(soU)  = 
0.01. 

ED  injection  was  also  assumed  to  start  in  2010  ( tEDo )  with  an  annual  average  injection  rate  based 
on  (7-1)  with  SF  =  2  yielding  MEd  =16  kg/d.  ED  injection  was  terminated  when  the  annually 
averaged  measured  (simulated  +  noise)  TCE  concentration  upgradient  of  the  ED  injection  gallery 
dropped  below  5  pg/L. 

The  Case  3  thermal  treatment  design  is  considerably  less  effective  than  the  optimized  designs  for 
Cases  1  and  2  even  though  the  source  mass  reduction  target  is  more  aggressive  than  nearly  all  of 
the  Case  1  and  2  targets  and  ED  injection  is  utilized  in  Case  3.  The  poor  performance  is 
attributable  to  a  high  probability  of  contaminant  mass  beyond  the  limited  thermal  treatment  area, 
which  is  60%  smaller  than  the  optimized  area.  Lower  thermal  costs  result,  but  there  is  a  high 
probability  of  source  mass  remaining  above  target  levels  after  thermal  treatment  that  cannot  be 
removed  by  natural  source  attenuation  by  2110.  This  results  in  a  71%  probability  of 
noncompliance  that  leads  to  long  durations  for  ED  injection  and  site- wide  monitoring  resulted 
with  a  total  ENPV  cost  of  $1 1.5M  -  more  than  six  times  higher  than  for  optimized  Case  la.  The 
cost  reduction  for  Case  la  relative  to  Case  3  is  about  84%. 

Long  duration  ED  injection  can  maintain  concentrations  at  the  compliance  location  below  the 
target  level  with  a  high  probability  (Figure  7-9a).  However,  source  mass  flux  is  unlikely  to  drop 
low  enough  to  shut  down  the  ED  injection  gallery  prior  to  the  maximum  simulation  date  of  2 1 10, 
leading  to  a  noncompliance  condition.  The  cost  probability  distribution  has  a  positive  bias 
(Figure  7-9b)  reflecting  the  large  number  of  cases  with  operating  costs  continuing  up  to  21 10. 


Year  NPV, 


Figure  7-9.  Results  for  Case  3:  (a)  TCE  concentrations  at  compliance  well 
and  (b)  NPV  cost  distribution  (without  penalty  cost). 
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7.3.3  Effect  of  compliance  rules  (Case  4) 

To  evaluate  the  effect  of  compliance  rules  on  remediation  design,  cost  and  performance,  we 
performed  an  optimization  identical  to  Case  la  except  that  the  EXV  rule  described  in  Chapter  6 
was  used  in  lieu  of  the  RCL  rule.  The  remedial  design  was  similar  to  that  for  Case  la  with  no  ED 
injection.  The  optimized  thermal  treatment  area  was  somewhat  less  aggressive,  while  mass 
reduction  targets  and  Efrac(init)  were  slightly  more  aggressive.  The  net  result  was  a  very  slightly 
lower  average  thermal  treatment  cost  and  a  slightly  higher  average  site-wide  monitoring  duration 
and  cost  (Table  7-3).  The  95%  upper  confidence  limit  for  concentration  versus  time  shift  a  few 
years  to  the  right,  while  the  99%  confidence  limit  becomes  nearly  flat  above  the  target 
concentration  all  the  way  to  2110  (Figure  7-10).  The  more  stringent  and  “noisy”  EXV  rule 
results  in  a  NC  probability  that  is  5  percentage  points  higher  than  for  Case  la. 
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Figure  7-10.  Results  for  Case  4:  (a)  TCE  concentrations  at  compliance  well 
and  (b)  NPV  cost  distribution  (without  penalty  cost). 


7.3.4  Effect  of  source  delineation  uncertainty  on  thermal  performance  (Case  5) 

To  evaluate  the  effect  of  uncertainty  in  source  zone  delineation  on  remediation  design,  cost  and 
performance,  we  first  performed  an  optimization  identical  to  Case  la  except  that  the  thermal 
treatment  area  was  fixed  at  the  best  estimate  of  the  source  area  from  calibration  and  not 
optimized.  However,  the  treatment  area  was  assumed  to  fully  encompass  the  region  with 
DNAPL  and  no  “missing  mass”  outside  the  treatment  area  was  simulated.  This  optimization  is 
designated  as  Case  5  a. 

The  optimized  design  variables  determined  in  Case  5  a  were  then  employed  in  a  MC  simulation 
without  further  optimization  with  “missing  mass”  outside  the  assumed  treatment  zone  treated  as 
a  stochastic  variable.  This  simulation  is  designated  as  Case  5b. 

Case  5a  yields  predicted  performance  comparable  to  Case  la.  The  probability  of  NC  at  3%  is  the 
same  as  Case  la  and  concentration  versus  time  probabilities  are  similar  (Figure  7-11).  The 
ENPV  cost  for  thermal  treatment  is  lower  due  to  the  smaller  system  footprint,  leading  to  a  lower 
simulated  total  ENPV  cost. 
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However,  Case  5a  is  unrealistically  optimistic  because  it  disregards  uncertainty  in  source  zone 
delineation  and  hence  in  the  thermal  treatment  area  necessary  to  reliably  reduce  the  source  mass. 
The  effects  of  this  Pollyannaish  treatment  of  source  delineation  reliability  is  seen  in  Case  5b 
results,  which  reveal  a  high  probability  of  contaminant  concentrations  above  target  levels  (Figure 
7-12)  and  a  53%  probability  of  NC  by  21 10  (Table  7-3). 

While  thermal  source  treatment  is  generally  very  effective  at  reducing  contaminant  mass  within 
the  treated  soil  volume,  it  is  sensitive  to  error  in  source  delineation  and  disregarding  this  fact  can 
lead  to  much  poorer  performance  than  expected. 


Year  NPV, 


Figure  7-11.  Results  for  Case  5a:  (a)  TCE  concentrations  at  compliance  well 
and  (b)  NPV  cost  distribution  (without  penalty  cost). 
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Figure  7-12.  Results  for  case  5b:  (a)  TCE  concentrations  at  compliance  well 
and  (b)  NPV  cost  distribution  (without  penalty  cost). 
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7.4  Conclusions 

Stochastic  cost  optimization  was  found  to  decrease  the  expected  NPV  cost  for  site  remediation  as 
much  as  85%  compared  to  conventional  non-optimized  approaches,  while  also  increasing  the 
probability  of  achieving  “no  further  action”  status  in  a  specified  timeframe  by  more  than  60%. 

Optimizing  monitoring  frequency  for  compliance  wells  used  to  make  no  further  action 
determinations  as  well  as  operational  monitoring  used  to  make  decisions  on  individual 
remediation  system  components  involve  tradeoffs  between  increased  direct  costs  for  sampling 
and  analysis  versus  decreased  construction  and  operating  costs  that  arise  because  more  data 
increases  decision  reliability,  which  ameliorates  the  extent  to  which  systems  must  be 
overdesigned  to  overcome  data  uncertainty. 

Operational  monitoring  and  heating  unit  shutdown  protocols  for  thermal  source  treatment 
(incremental  vs.  all-or-none  shutdown,  soil  vs.  groundwater  sampling,  number  and  frequency  of 
samples)  can  significantly  affect  system  effectiveness  and  cost.  Optimizing  these  protocols  alone 
was  found  to  effect  cost  savings  of  more  than  20%. 

The  statistical  formulation  of  compliance  rules  can  also  significantly  affect  remediation  design, 
performance  and  cost.  We  found  that  defining  compliance  based  on  95%  confidence  limits  of  a 
regression  on  a  moving  time  window  was  more  reliable  than  using  measured  extreme  on  a 
similar  lookback  period,  resulting  in  decreased  expected  cost  and  lower  probability  of  failure. 

Disregarding  uncertainty  in  DNAPL  source  delineation  can  lead  to  marked  undersign  of  thermal 
treatment  systems,  which  can  result  in  a  low  probability  of  achieving  remediation  objectives  and 
ultimately  to  higher  costs  to  later  rectify  the  misjudgments. 

Application  of  the  SCOToolkit  methodology  on  a  given  site  should  increase  the  probability  of 
meeting  cleanup  objectives  with  the  least  cost.  Employed  over  a  large  number  of  sites  should 
yield  significant  total  cost  reductions. 
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_ Chapter  8 _ 

Application  to  Fort  Lewis  EGDY  site 


8.1  Site  Description 

The  East  Gate  Disposal  Yard  (EGDY)  is  a  source  of  groundwater  contamination  at  the  Logistics 
Center  National  Priority  List  Site  located  at  Fort  Lewis,  Washington  (Figure  8-1).  EGDY  was 
used  between  1946  and  the  mid  1970s  as  a  waste  disposal  site  for  solvents  from  cleaning  and 
degreasing  operations.  Material  was  transported  to  the  disposal  yard  in  barrels  and  vats  from 
various  areas.  About  seven  barrels  of  liquid  waste  per  month  were  disposed  during  peak 
operation.  A  TCE  plume  in  the  shallow  aquifer  evolved  from  the  disposal  site  with 
concentrations  in  the  range  of  hundreds  pg/L  in  the  source  area  and  concentrations  exceeding  5 
pg/L  over  4  km  downgradient  (Dinicola,  2005;  USAGE,  2008). 

The  climate  of  Fort  Lewis  is  characterized  by  warm  dry  summers  and  cool  wet  winters  with  a 
mean  annual  temperature  of  about  13  degrees  Celsius  and  mean  annual  precipitation  of  about 
1000  mm.  Fort  Lewis  is  underlain  by  a  complex  and  heterogeneous  sequence  of  glacial  and  non¬ 
glacial  deposits  including  a  shallow  aquifer  (Vashon)  and  a  deep  aquifer  (Sea  Level  Aquifer, 
SLA).  The  Vashon  aquifer  is  unconfined  and  continuous  throughout  the  Fort  Lewis  area.  It 
ranges  in  thickness  between  about  30  to  60  meters.  The  Vashon  and  SLA  aquifers  are  separated 
by  a  mostly  continuous  low  permeability  aquiclude.  However,  a  “window”  occurs  about  2  km 
downgradient  of  the  disposal  area  that  allows  water  and  contaminants  from  the  shallow  Vashon 
aquifer  to  migrate  to  the  deep  SLA  aquifer. 


Figure  8-1.  Location  of  EDGY  site  and  TCE  plumes  as  of  2004  (Dinicola,  2005). 
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Groundwater  at  Fort  Lewis  generally  flows  to  northwest  in  the  Vashon  aquifer  and  west- 
southwest  in  the  SLA  aquifer.  A  simplified  geologic  cross  section  of  the  Fort  Lewis  site  is  shown 
in  Figure  8.2.  More  details  on  the  site  geology  are  found  in  Dinicola  (2005),  Truex  et  al.  (2006), 
and  USACE  (2008). 
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Figure  8-2.  Hydrogeologic  section  of  Fort  Lewis  (Dinicola,  2005). 


Several  remediation  actions  have  been  performed  at  the  EGDY  site  to  contain  the  existing 
contaminant  plume  or  reduce  DNAPL  mass  in  the  source  zone.  Disposal  trenches  were 
excavated  in  2000  to  remove  contaminated  waste  buried  above  the  water  table.  About  1260 
drums  of  contaminant  were  removed.  To  reduce  DNAPL  mass  below  the  water  table,  TSR  was 
implemented  for  three  source  zones  between  late  2003  and  early  2007.  Tables  8-1  and  8-2 
summarize  the  activities. 
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Table  8-1.  EGDY  site  remediation  history  (USAGE,  2008) 


Date 

Activity 

Location 

1995  -2005 

Pump-and-treat  systems  installed  in  Vashon 
Aquifer 

One  near  EGDY 
second  near  highway  1-5 

2003  -  2005 

Integrated  pump  test  in  Areas  1  and  3  in  Nov 
2003  and  Sep  2005,  respectively 

EGDY 

2003  -  2005 

Source  flux  measurements  in  Areas  1  and  3  in 
Nov  2003  and  Sep  2005,  respectively 

EGDY 

2003  -  2006 

TSR  and  monitoring  in  Areas  1,  2  and  3  in  Dec 
2003  -  Aug  2004,  Feb  2005  -  Aug  2005,  and  Oct 
2006  -  Jan  2007,  respectively 

EGDY 

2005-2006 

Whey  injection  pilot  tests 

EGDY 

2005  -  2007 

Post-TSR  monitoring  in  Areas  1 ,2  and  3  in  May 
2005,  Sep.  2005,  and  Feb  2007,  respectively 

EGDY 

2006  -  2008 

Post-treatment  soil  coring  in  Areas  1,2  and  3  in 
Apr  2006,  Apr  2006,  and  Mar  2008, 
respectively 

EGDY 

2009 

Pump-and-treat  system  installed  in  Sea  Level 
Aquifer 

Near  hospital 

Table  8-2.  Summary  of  TSR  operations  at  EGDY  site  (USAGE,  2008) 


Variable 

Area  1 

Area  2 

Area  3 

TSR  treatment  area  (m  ) 

2360 

2080 

1691 

TSR  max  depth  below  ground  surface  (m) 

10 

16 

9 

TSR  treatment  volume  (m3) 

23625 

135953 

15368 

Energy  on  date 

12/17/2003 

02/14/2005 

10/11/2006 

Energy  off  date 

08/04/2004 

08/05/2005 

01/26/2007 

Duration  (days) 

231 

172 

107 

Mass  removal,  TCE  +  DCE  (kg) 

2990 

1340 

1120 

8.2  Model  Calibration 

8.2.1  Characterization  of  groundwater  flow  field 

Groundwater  flow  at  the  Fort  Lewis  site  was  characterized  in  USAGE  (2008).  The  model  used  in 
this  project  simulates  groundwater  flow  and  transport  with  curvilinear  streamlines  as  described  in 
Chapter  2.  Streamlines  commencing  from  each  DNAPL  source,  actual  or  planned  ED  injection 
gallery  and  from  the  “window”  between  the  upper  and  lower  aquifer  units  were  digitized  and 
fitted  to  third  order  polynomial  equations  of  the  form  y  =  ax+bx2 +cx3 .  The  model  computes 
travel  distances  from  sources  to  the  wells  of  interest  along  their  streamlines.  Coefficients  of 
individual  streamline  equations  are  presented  in  Figure  8-3.  Since  the  locations  of  ED  galleries  1 
to  3  are  immediately  upgradient  of  Areas  1  to  3,  their  streamlines  are  similar  to  those  of  Areas  1 
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to  3.  The  streamline  of  ED  gallery  4  starts  upgradient  of  the  “window”  and  follows  the 
streamline  of  Area  2.  Polynomial  coefficients  can  be  estimated  using  the  provided  function 
‘trajectory.exe’  or  using  Excel  to  solve  the  nonlinear  least  squares  problem  using  the  SOLVER 
function. 


Origin 

a 

b 

c 

Area  1  (Source  1) 

4.275E-01 

5.61 5E-04 

2.756E-08 

Area  2  (Source  2) 

2.605E-01 

7.257E-04 

6.437E-08 

Area  3  (Source  3) 

1.129E-01 

1.244E-03 

2.758E-07 

Win  1 

-9.570E-01 

-8.070E-04 

-1.510E-07 

Win  2 

-7.770E-01 

-6.800E-04 

-1.230E-07 

Win  3 

-1.171E+00 

-9.984E-04 

-1.780E-07 

ED  1 

4.952E-01 

5.81  IE-04 
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Figure  8-3.  Streamlines  used  to  model  groundwater  flow  at  the  Fort  Lewis  site. 


8.2.2  Calibration  using  pre-  and  post-ERH  data 

Chlorinated  solvent  concentrations  in  groundwater  reported  by  Truex  et  al.  (2006)  were  utilized 
to  construct  time-series  for  each  monitoring  well.  Chlorinated  solvent  species  were  converted  to 
“TCE-equivalent”  concentrations  such  that  their  H-demand  for  complete  reduction  is  equal  to 
that  of  TCE.  The  sum  of  H-equivalent  concentrations  of  all  chlorinated  ethenes  was  taken  as  the 
total  TCE-equivalent  solvent  concentration.  Locations  of  monitoring  wells  with  data  used  for 
model  calibration  are  shown  in  Figure  8-4.  Pre-ERH  data  included  dissolved  concentration 
measurements  from  26  wells  from  1995  to  late  2001.  Post-ERH  data  included  longer  time  series 
for  pre-ERH  wells  through  2009,  data  from  14  additional  wells  closer  to  the  source  areas  (Figure 
8-4  inset),  and  measurements  of  source  mass  flux  and  mass  removed  by  ERH. 
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Figure  8-4.  Well  locations  used  in  calibration. 


To  model  source  zone  mass  dissolution  and  transport  downstream,  initial  estimates  of  source  and 
aquifer  parameters  were  estimated  from  information  in  various  reports  (Dinicola,  2005;  Truex  et 
al.,  2006;  USACE,  2008)  summarized  in  Table  8.3.  We  calibrated  model  parameters  to  site  data 
using  2000  as  the  reference  year  ( tcat )  for  source  mass  and  source  flux  using  both  pre-  and  post¬ 
remediation  data.  Data  include  dissolved  concentration  data  (total  TCE-equivalents),  source  flux 
measurements,  and  DNAPL  mass  removed  by  TSR.  Details  of  calibration  data  are  tabulated  in 
Appendix  A. 

We  assume  that  most  of  the  contamination  present  in  the  SLA  was  transported  from  the  Vashon 
aquifer  through  the  “window”  in  the  aquiclude.  The  Vashon  aquifer  is  divided  into  three  zones 
with  different  decay  coefficients.  Zone  1  extends  from  the  DNAPL  sources  to  the  window;  Zone 
2  encompasses  the  window  itself;  and  Zone  3  is  the  region  downgradient  of  the  window.  The 
model  internally  computes  a  first-order  “decay”  coefficient  for  the  window  zone  as  described  in 
Chapter  2  to  account  for  the  advective  flux  to  the  SLA  unit.  The  window  zone  is  divided  laterally 
into  three  sections  to  approximate  the  average  TCE  flux  from  the  upper  to  lower  aquifer.  Fluxes 
are  tabulated  at  discrete  times  in  a  lookup  table  and  interpolated  to  define  a  smooth  source 
function  for  the  Vashon  aquifer.  The  model  uses  0.0001  d'1  as  a  prior  estimate  of  biodecay 
coefficients  of  for  Zones  1  and  3. 
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Table  8-3.  Fort  Lewis  site  characterization  data. 


Area 

Parameters 

Prior  Value1 

STD2 

Reference 

Area  1 

Mass  at  2000  (kg) 

10330 

0.63 

USAGE,  2008 

Flux  at  2000  (kg/d) 

0.75 

1.00 

USAGE,  2008 

Release  date 

1970 

3.00 

USAGE,  2008 

Width  (m) 

47 

0.20 

USAGE,  2008 

Area  2 

Mass  at  2000  (kg) 

6495 

0.56 

USAGE,  2008 

Flux  at  2000  (kg/d) 

0.32 

1.00 

USAGE,  2008 

Release  date 

1970 

3.00 

USAGE,  2008 

Width  (m) 

42 

0.20 

USAGE,  2008 

Area  3 

Mass  at  2000  (kg) 

7987 

0.60 

USAGE,  2008 

Flux  at  2000  (kg/d) 

0.42 

1.00 

USAGE,  2008 

Release  date 

1970 

3.00 

USAGE,  2008 

Width  (m) 

34 

0.20 

USAGE,  2008 

Vashon 

q  (m/d) 

0.4 

0.15 

Truex  et  al.,  2006;  Dinicola,  2005 

Porosity 

0.29 

- 

Truex  et  al.,  2006;  Dinicola,  2005 

Retardation 

1.21 

0.15 

Typical  from  literature 

Longitudinal  dispersivity  (m) 

20 

0.20 

Estimated  from  plume  length 

Transverse  dispersivity  (m) 

2 

0.20 

Estimated  fcrom  plume  length 

Saturated  depth  (m) 

30 

0.20 

Truex  et  al.,  2006 

ED  average  (H-eq  ppb) 

48 

0.15 

Dinicola  2005 

EA  average  (H-eq  ppb) 

1977 

0.15 

Dinicola  2005 

SLA 

Darcy  velocity  (m/d) 

1 

0.15 

Truex  et  al.,  2006;  Dinicola,  2005 

Porosity 

0.22 

- 

Truex  et  al.,  2006;  Dinicola,  2005 

Retardation 

1.27 

0.15 

Typical  from  literature 

Longitudinal  dispersivity  (m) 

20 

0.20 

Estimated  from  plume  length 

Transverse  dispersivity  (in) 

2 

0.20 

Estimated  from  plume  length 

Saturated  depth  (m) 

30 

0.20 

Truex  et  al.,  2006 

Avg  natural  ED  (H-eq  ppb) 

768 

0.15 

Dinicola,  2005 

Avg  natural  EA  (H-eq  ppb) 

1977 

0.15 

Dinicola,  2005 

Window 

qz  (m/d) 

0.05 

0.50 

Truex  et  al.,  2006 

1  Prior  estimates  represent  arithmetic  mean  for  release  date,  geometric  mean  for  other  parameters. 

2  Standard  deviations  of  prior  estimates  are  dimensionless  statistics  for  ln-transformed  values  for  all  parameters 


except  release  dates,  which  are  in  actual  units. 
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Model  calibration  yields  refined  estimates  of  model  parameters  as  described  in  Chapter  5.  Final 
estimates  of  calibrated  parameters  are  tabulated  in  Table  8-4.  Residual  uncertainty  in  dissolved 
concentration  predictions  (Si„c)  was  estimated  as  0.79  for  the  pre-ERH  calibration  and  0.77  with 
post-ERH  data  indicating  that  additional  data  during  and  after  ERH  improved  model  calibration. 
Comparisons  of  observed  and  predicted  concentrations  for  pre-  and  post-ERH  calibration  are 
shown  in  Figure  8-5.  The  most  notable  change  in  the  post-ERH  calibration  is  a  decrease  in  the 
best  estimate  of  total  source  mass  in  2000  and  a  decrease  in  source  mass  uncertainty  using  post- 
ERH  data.  Interestingly,  the  estimated  source  flux  decreased  while  its  uncertainty  increased. 

Table  8-4.  Calibration  summary  for  Fort  Lewis. 

Pre-ERH  Post-ERH 


Parameter 

Best 

STD 

Best 

STD 

Post-ERH  Data 

McaU  (kg) 

16580 

0.58 

3831 

0.05 

2990  kg  removed  by  Aug  2004 

Mcal2(  kg) 

9851 

0.55 

5339 

0.54 

1340  kg  removed  by  Aug  2005 

Mcai  3  (kg) 

3403 

0.20 

5133 

0.25 

1 120  kg  removed  by  Jan  2007 

MCal  sum  (kg) 

29834 

0.37 

14302 

0.22 

J cal  1  (kg/d) 

0.81 

0.05 

0.43 

0.23 

0.75  kg/d  in  Nov  2003 

Jcal  2  (kg/d) 

0.34 

0.03 

0.06 

0.36 

0.32  kg/d  in  Nov  2003 

Jcal3  (kg/d) 

0.72 

0.10 

0.50 

0.22 

0.42  kg/d  in  Apr  2006 

J cal  sum  (kg/ d) 

1.87 

0.04 

0.99 

0.15 

Additional  monitoring  well  data 

Pi 

0.924 

0.25 

0.982 

0.24 

P2 

1.001 

0.25 

1.030 

0.25 

Ps 

1.150 

0.23 

1.090 

0.24 

Ai  (m2) 

2748 

0.20 

2736 

0.20 

A 2  (m2) 

2754 

0.20 

2822 

0.20 

A3  (m2) 

2579 

0.19 

6774 

0.08 

Notes:  Calibration  date  tcai  is  2000  for  both  pre-  and  post-ERH  calibrations.  Best  values  denote  parameter  best 
estimates.  STD  values  are  posterior  estimates  of  standard  deviations  of  ln-transformed  values  for  all  parameters 
except  P  which  is  not  transformed.  A,  denotes  the  area  of  source  zone  i. 


Observation,  ppb  Observation,  ppb 


Figure  8-5.  Observed  vs.  simulated  concentrations  for  pre-  and  post-ERH  calibrations.  Blue  (o)  = 

Vashon  aquifer;  Red  (+)  =  SLA  aquifer. 
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8.3  Remedial  Design  Evaluation  and  Optimization 

8.3.1  Long-term  Monte  Carlo  simulations  with  no  further  action 

Non-optimized  MC  simulations  were  performed  to  assess  the  long-term  effectiveness  of  ERH 
source  remediation  already  performed  at  the  site  based  on  pre-  and  post-remediation  parameter 
estimates.  Simulations  of  ERH  were  made  such  that  the  contaminant  mass  removed  from  each 
source  was  equal  to  the  measured  removal.  The  following  simulations  were  performed: 

•  NoOpt  la  -  No  ERH  or  any  subsequent  remedial  actions  based  on  pre -ERH  calibration, 

•  NoOpt  lb  -  Actual  ERH  with  no  further  remedial  actions  using  pre -ERH  calibration, 

•  NoOpt  2a  -  No  ERH  or  any  subsequent  remedial  actions  based  on  post-ERH 
calibration,  and 

•  NoOpt  2b  -  Actual  ERH  with  no  further  remedial  actions  usingpost-ERH  calibration. 

Long  term  performance  for  the  NoOpt  lb  case  differs  very  little  from  NoOpt  la  with  the 
predicted  time  for  the  median  concentration  to  decrease  to  5  pg/L  only  10  years  earlier  in 
response  to  ERH  (Figure  8-6). 


Year 


Year 

Figure  8-6.  TCE-equivalent  concentration  at  the  compliance  well  before  (NoOpt  la,  upper  panel) 
and  after  actual  source  removal  (NoOpt  lb,  lower  panel)  based  on  pre-ERH  calibration 
parameters  (red  arrow  indicates  median  time  to  reach  MCL). 
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The  small  predicted  impact  of  ERH  can  be  attributed  to  the  small  fraction  of  the  calibrated 
source  volume  that  is  recovered  by  the  actual  ERH  system,  which  is  at  least  partially  due  to 
overestimation  of  source  mass  in  the  pre-ERH  calibration.  The  confidence  bands  for  the  pre- 
ERH  simulations  are  rather  wide  and  increase  over  time. 

Confidence  bands  for  MC  simulations  based  on  post-ERH  calibration  parameters  are  narrower, 
reflecting  the  higher  reliability  of  parameter  estimates  conditioned  on  post-ERH  data  (Figure  8- 
7).  Long-term  performance  for  NoOpt  2b  shows  the  predicted  time  for  the  median  concentration 
to  reach  the  MCL  is  about  80  years  shorter  than  for  the  no  action  NoOpt  2a  case. 

Nevertheless,  with  a  greater  than  50%  probability  of  taking  more  than  100  years  to  reach  the 
MCL  following  ERH  with  no  additional  remedial  action,  the  performance  of  thermal  treatment 
conducted  at  the  site  taken  cannot  be  judged  an  unqualified  success. 


Year 


Year 


Figure  8-7.  TCE-equivalent  concentration  at  the  compliance  well  before  (NoOpt  2a,  upper  panel) 
nd  after  actual  source  removal  (NoOpt  2b,  lower  panel)  based  on  post-ERH  calibration 
parameters  (red  arrow  indicates  time  to  reach  MCL). 
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8.3.2  ERH  design  optimization 

The  small  mass  recovery  by  ERH  at  the  EGDY  site  relative  to  the  post-ERH  calibrated  source 
mass  suggests  that  a  significant  quantity  of  contaminant  mass  occurs  beyond  the  limits  of  the 
various  ERH  source  treatment  volumes.  Poor  perfonnance  may  also  be  due  to  a  certain  extent  to 
less  than  optimal  mass  recovery  within  the  treatment  zones.  To  evaluate  the  potential  for 
improvement  in  ERH  performance,  we  turn  our  attention  in  this  section  to  ERH  design 
optimization.  We  will  base  the  optimization  analysis  on  the  post-ERH  calibration,  which  is 
expected  to  be  significantly  more  reliable  than  the  pre-ERH  calibration. 


Table  8-5.  Cost  variables  used  in  design  optimization. 


Description 

Variable 

Value 

Unit 

Thennal  treatment  costs 

Fixed  cost  for  site 

s~*TR 
^  site 

3201 

$K 

Thermal  system  volume  cost  multiplier 

siTR 

volt 

0.09751 

$K/m3 

Thermal  system  area  cost  multiplier 

s~iTR 
^  area  i 

0.73321 

$K/m2 

Cost  per  soil  boring 

s~iTR 
^ boref 

0.1 

$K 

Analysis  cost  per  soil  sample 

s~iTR 

^  SOILs  amp 

0.1 

$K/sample 

Source  zone  monitoring  well  installation  cost 

siTR 

wellf 

0.2 

$K 

Sampling  and  analysis  costs  per  groundwater  sample 

ED  iniection  system  costs 

s-iTR 

GWsamp 

0.5 

$K/sample 

Fixed  cost  per  ED  injection  gallery  width 

s~iEDcap 
^ width 

3.32 

$K/m 

Fixed  cost  per  ED  performance  monitoring  well 

siEDcap 

mw 

102 

$K/well 

Other  fixed  costs  for  ED  injection 

siEDcap 

other 

552 

$K 

Operating  cost  per  mass  of  ED  injected 

s~*EDop 

mass 

0.0074272 

$K/kg 

Operating  cost  per  year  per  meter  ED  gallery  width 

s~iEDop 
^  width 

0.7332 

$K/yr/m 

Cost  per  sample  for  ED  performance  monitoring 

s~ lEDop 
^  samp 

2.52 

$K/sample 

Other  operating  cost  per  year  for  ED  injection 

s~iEDop 
^ all 

61. 12 

$K/year 

Other  costs 

Cost  per  compliance  sampling  event 

fiSWop  -ktSW 
^  samp  v  well 

2.5 

$K/event 

PT  operating  costs  and  other  site-wide  costs 

s-iPTop 
^  total 

87.5 

$K/year 

Penalty  cost  for  non-compliance 

/o  pen 
^ NPV 

50000 

$K 

Reference  year  for  cost  discounting 

tref 

2010 

year 

year'1 

Discount  rate 

d 

0.03 

Based  on  G.  Bayke,  TRS  Inc.,  personal  communication  (2009)  Note:  values  do  not  include  source 
characterization  costs,  post-remediation  flux  measurements  or  groundwater  pumping  at  source 

2  Sorenson  and  Macbeth  (2009) 

3  J.  Gillie,  CTR  US  USA  IMCOM,  personal  communication  (2010) 
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The  optimization  simulation  assumes  ERH  operational  monitoring  following  Method  la  (all 
heating  units  within  a  given  source  treatment  zone  are  turned  off  at  the  same  time  when  the 
average  soil  concentration  drops  below  a  target  fraction  of  the  initial  concentration  -  see  Chapter 
4)  and  compliance  is  based  on  the  EXV  rule  (see  Chapter  6)  for  a  5 -year  lookback  period 
(N lookback)  with  quarterly  groundwater  sampling  ( f™mp  =  4).  A  cost  penalty  (  CF™V  )  is  applied  for 
NC  after  tpenaity  =  2050  or  if  cleanup  criteria  are  not  met  by  tmax  =  2110. 

Cost  variables  used  for  design  optimization  are  summarized  in  Table  8-5.  Uncertainty  in  all  unit 
costs  were  assumed  to  follow  a  uniform  probability  distribution  with  a  range  from  0.9  to  1.1 
times  the  best  estimates. 

Optimized  ERH  design  variables  are  summarized  in  Table  8-6.  Optimal  mass  reduction  factors 
(Rtsr(soU))  range  from  1.35  to  4.59%  of  the  post-ERH  estimates  of  source  mass.  Optimized 
thermal  treatment  areas  range  from  1.35  times  the  actual  treatment  areas  used  at  the  site  during 
actual  operations  in  Area  1  to  3.60  times  in  Area  3.  The  high  ratio  for  Area  3  suggests  a  large 
amount  of  mass  outside  the  actual  treatment  zone,  which  is  consistent  with  the  high  observed 
TCE  concentrations  in  wells  near  Area  3  after  thermal  treatment. 

The  model  assumes  that  uncertainty  in  untreated  mass  is  directly  proportional  to  uncertainty  in 
source  area  and  hence  optimizes  treatment  area  to  reduce  the  quantity  of  untreated  mass  beyond 
the  treatment  perimeter  and  thus  decrease  the  duration  and  cost  of  dissolved  plume  remediation 
and  the  probability  of  NC.  While  this  assumption  is  simplistic  (mass  may  be  distributed 
nonuniformly  with  small  “hot  spots”),  large  uncertainty  in  source  area  and  large  spreads  between 
optimized  thermal  treatment  area  and  the  best  estimate  of  source  area  should  raise  a  warning  that 
additional  characterization  may  be  needed  to  delineate  the  source  more  accurately  unless 
characterization  costs  exceed  the  cost  of  overdesigning  the  thermal  system. 

The  NPV  cost  ranged  from  $20M  to  26M  with  an  expected  value  of  $22. 9M  (Figure  8-8),  not 
including  costs  incurred  prior  to  2010  or  penalty  costs.  The  ENPV  for  ERH  only  is  about  $22M, 
which  is  75%  higher  than  the  cost  of  the  actual  ERH  system  (~$  1 2.6M  excluding  groundwater 
pumping  and  source  flux  measurements),  primarily  reflecting  the  larger  areal  extent  of  the 
optimized  system.  NC  probability  was  estimated  to  be  on  the  order  of  1%  and  the  99%  upper 
confidence  limit  of  the  time  to  reach  NFA  was  about  30  years  after  commencing  ERH  (Figure  8- 
9)..  The  results  indicate  that  more  refined  source  characterization,  especially  in  the  vicinity  of 
Area  3,  followed  by  additional  thermal  treatment  of  “missed”  areas  could  significantly  increase 
the  odds  of  success  compared  with  the  actual  ERH  system  that  was  implemented  at  the  site. 


Table  8-6.  Optimization  of  ERH  design  for  EGDY  DNAPL  sources. 


Source 

2 

Thermal  Treatment  Area  (Ajsr,  m  ) 

RTSR(soil)  (%) 

Optimized  Value 

Actual  System 

Area  1 

3175 

2360 

4.59 

Area  2 

3108 

2080 

1.36 

Area  3 

6099 

1691 

1.35 
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Figure  8-8.  Cost  distribution  of  optimized  ERH  system  excluding  penalty  costs. 


Year 

Figure  8-9.  TCE-equivalent  concentration  at  the  compliance  point  based  on  optimized  ERH 
variables  in  Fort  Lewis  (red  arrow  indicates  99%  confidence  limit  of  time  to  reach  MCL). 

8.3.3  Optimization  of  enhanced  bioremediation  with  whey  injection 

Model  calibration  and  optimization  analyses  indicate  that  significant  amounts  of  DNAPL  were 
not  recovered  by  past  ERH  operations  at  the  site,  especially  near  Area  3.  Considering  the  high 
cost  of  a  second  round  of  source  characterization  and  thermal  treatment,  we  will  investigate  the 
effectiveness  of  employing  whey  injection  to  accelerate  source  zone  DNAPL  dissolution  and  to 
enhance  dissolved  phase  anaerobic  dechlorination.  Pilot  studies  in  Area  3  indicate  that  whey 
injection  increases  DNAPL  dissolution  kinetics  (Macbeth  and  Sorenson,  2008). 

We  consider  four  potential  whey  injection  gallery  locations:  three  upgradient  of  each  of  the  three 
source  areas,  designated  EDI  through  ED3,  and  a  fourth  upgradient  of  the  “window”  from  the 
shallow  to  deep  aquifer,  designated  as  ED4.  The  purpose  of  EDI  to  ED3  is  two-fold.  First,  we 
wish  to  accelerate  DNAPL  dissolution,  and  second,  we  wish  to  enhance  downgradient  aqueous 
plume  biodecay.  The  purpose  of  ED4  is  solely  to  enhance  aqueous  plume  biodecay  and  thus 
reduce  the  flux  of  contaminant  downgradient  of  the  gallery  in  the  Vashon  aquifer  and  to  reduce 
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the  contaminant  flux  through  the  window  to  the  SLA.  We  propose  to  monitor  total  TCE  and 
daughter  products  in  performance  monitoring  wells  downgradient  of  the  respective  ED  gallery 
and  source  area  and  terminate  injection  when  the  source  is  depleted  to  the  point  that  the  annual 
average  concentration  is  less  than  a  specified  concentration  Cedi  for  gallery  i.  For  the  ED4 
gallery,  we  propose  to  monitor  a  location  upgradient  of  ED4  and  to  shut  off  injection  after  the 
annual  average  concentration  drops  below  Ced4  and  all  upgradient  ED  galleries  have  been  off 
longer  than  the  max  travel  time  to  ED4.  A  single  compliance  well  is  employed  for  optimization 
simulations,  which  is  on  the  plume  centerline  downgradient  of  ED4  and  upgradient  of  the 
“window”  (Figure  8-10). 


DNAPL  sources 


Figure  8-10.  Conceptual  model  of  ED  injection  gallery  configuration,  ED  performance 
monitoring  and  compliance  monitoring  locations. 


Mass  transfer  enhancement  is  characterized  by  a  factor  Fmt  that  represents  the  ratio  of  mass 
transfer  under  “natural”  conditions  to  that  with  ED  injection  as  described  in  Chapter  2.  The  mass 
transfer  enhancement  factor  that  is  assumed  to  increase  proportional  the  ED  concentration  as 
Fmt  =  1  +  fjC£ .  An  estimate  of  fmt  was  made  from  a  pilot  test  performed  at  the  site  (Area  3 

treatment  cell  2)  in  which  a  10%  whey  solution  was  injected  in  well  IW-2,  screened  over  a  5  m 
interval,  at  a  rate  of  1800  gal/month  (22.7  kg  whey/day)  from  Nov  2005  to  Feb  2006  (Macbeth 
and  Sorenson,  2008).  Fmt  was  determined  for  wells  2A1  and  2D1  which  were  3  and  10  m, 
respectively  downgradient  of  IW-2,  as  the  ratio  of  measured  TCE  concentration  after  whey 
injection  to  that  prior  to  injection.  A  simulation  of  whey  injection  and  transport  for  the  pilot  test 
was  performed  to  estimate  whey  concentrations  C%  at  each  well  and  fmt  was  computed 
accordingly.  Results  summarized  in  Table  8-7  yield  an  average  value  of fmt  =  10°  L/pg. 


Table  8-7.  COD  measured  during  whey  injection  pilot  test. 


Observation  well 

Distance  from  IW-2 

Fmt 

fmt  (L/pg) 

2A1 

3  m 

5.6 

0.8xl0'5 

2D1 

10  m 

4.5 

1.3xl0'5 
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The  first  date  of  whey  injection  for  optimization  simulations  was  assumed  to  be  2012.  Since 
ERH  was  completed  prior  to  2010,  ERH  costs  were  not  considered  in  the  ED  cost  optimization. 
ERH  was  simulated  by  simply  removing  the  measured  ERH  mass  recovery  from  Area  1-3 
sources.  ED  injection  gallery  widths  were  fixed  at  70,  90,  and  90  m  for  galleries  1-3  (slightly 
wider  than  the  calibrated  source  widths)  and  400  m  for  gallery  4  (about  half  the  5  pg/L  plume 
width).  One  kilogram  of  whey  is  assumed  to  have  a  reducing  potential  equivalent  to  0.078  kg  of 
hydrogen  (Table  3-1). 

Whey  injection  rate  (Med),  TCE-equivalent  concentrations  to  terminate  ED  injection  ( Ced ),  and 
ED  performance  monitoring  frequency  ( f  ™  )  were  optimized  as  decision  variables.  A  penalty 

date  (t penalty)  was  set  as  2100  and  the  maximum  simulation  date  (tmax)  as  21 10  since  ED  injection 
is  considered  a  long  term  plume  control  strategy,  but  simulation  was  extended  to  2300  to  assess 
longer-term  perfonnance.  The  compliance  criteria  to  terminate  monitoring  is  for  the  95%  RCL  of 
total  TCE-equivalent  concentration  for  a  5-year  lookback  period  ( Ni00kback )  to  be  below  5  pg/L. 

Optimization  results  (Table  8-8)  indicate  only  galleries  1  and  2  should  be  operated  with  whey 
injection  rates  of  25.0  and  29.1  kg/d,  respectively,  which  correspond  to  1983  and  2339  gallons 
per  month  of  10%  whey  solution.  At  these  injection  rates  over  the  gallery  widths,  the  time- 
average  increase  in  source  mass  transfer  rates  are  less  than  about  20%.  Quarterly  monitoring  of 
ED  performance  wells  is  determined  to  be  optimal  with  stop  criteria  (Ced)  of  26  and  16  pg/L  for 
galleries  1  and  2,  respectively. 

The  ENPV  cost  for  the  optimized  whey  injection  design  including  costs  of  site-wide  monitoring 
and  existing  pump-and-treat  system  operation  but  excluding  penalty  costs  is  $4667K  with  a  93% 
probability  of  meeting  compliance  criteria  by  2100  and  NFA  by  2110  and  with  7%  NC 
probability  (penalty  costs  applied  in  optimization).  The  median  date  at  which  compliance  is  met 
is  about  2040  with  99%  probability  confidence  limits  ranging  from  2030  to  slightly  greater  than 
2300  (Figure  8-11).  The  range  in  NPV  cost  spans  from  about  $1.5M  to  $10M  (Figure  8-12). 


Table  8-8.  Optimization  of  ED  injection  after  actual  ERH  at  Fort  Lewis. 


ENPV* 

($K) 

ED  operation 

Average 

monitoring 

period 

(yr) 

Non-compliance  Gallery 
(%) 

Injection 

(kg/d) 

Ced 

(Pg/L) 

Avg  ED 
duration 

(yr) 

rED 

J  samp 

(yr-1) 

4667 

7  ED  1 

25.0 

26 

1.0 

5 

46 

ED  2 

29.1 

16 

19.7 

5 

ED  3 

- 

- 

- 

- 

ED  4 

- 

- 

- 

- 

*  ENPV  does  not  include  penalty  costs 
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2010  2042  2074  2107  2139  2171  2203  2236  2268  2300 

Year 


Figure  8-11.  Confidence  limits  of  TCE-equivalent  concentration 
at  the  compliance  point  versus  time. 


Figure  8-12.  Cost  probability  distribution  for  whey  injection  optimization. 
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8.4  Conclusions 


Optimization  of  thermal  treatment  for  the  Fort  Lewis  EGDY  site  using  pre-thennal  treatment 
calibration  parameters  indicated  that  the  thermal  treatment  areas  needed  to  be  much  larger  than 
the  calibrated  source  areas  to  successfully  capture  the  source  due  to  uncertainty  in  source 
delineation.  MC  simulations  of  the  thermal  system  actually  implemented  at  the  site  indicated  that 
its  probability  of  meeting  remediation  objectives  was  low.  The  actual  Implementation  of  a 
thermal  system  Recalibration  after  thermal  treatment  provided  refined  estimates  of  source  mass 
and  other  parameters  and  confirmed  that  the  thermal  treatment  system  implemented  at  the  site 
failed  to  capture  much  of  the  source  mass. 

Although  not  reported  here,  we  also  performed  model  calibration  for  the  three  known  sources 
plus  a  fourth  source  with  an  unknown  location.  Coordinates  of  the  fourth  source  as  well  as  other 
source  parameters  were  calibrated.  Forward  simulations  of  thermal  treatment  for  the  three  known 
sources,  even  assuming  100%  capture  by  implemented  thennal  systems,  predicted  similar  poor 
performance  due  to  mass  remaining  in  the  untreated  fourth  source.  Results  were  not  presented 
because  uncertainty  in  calibrated  coordinates  of  the  fourth  source  were  large  and  the  results  led 
to  the  same  conclusions  -  namely,  that  further  site  characterization  is  needed  to  better  delineate 
DNAPL  sources  within  the  EGDY  area. 

Prediction  uncertainty  was  reduced  by  recalibrating  the  model  using  data  obtained  during  and 
after  thennal  treatment,  suggesting  that  an  iterative  calibration-optimization  approach  over  time 
as  additional  data  is  collected  may  allow  remediation  system  design  and  operation  to  be  refined, 
increasing  the  probability  of  success  and  reducing  cost. 

Optimization  of  whey  injection  systems  indicated  that  source  zone  injection  galleries  could  be 
useful  in  conjunction  with  other  source  zone  mass  removal  technologies,  such  as  thermal 
treatment,  although  injection  at  rates  high  enough  to  effect  substantial  mass  transfer 
enhancement  did  not  appear  to  be  cost  effective.  Improved  source  delineation  and  targeted 
source  treatment  appear  to  be  the  most  pressing  needs  at  the  site. 

Refinement  of  the  site  model  calibration  with  additional  source  characterization  data  and 
stochastic  optimization  of  additional  source  remedial  actions,  enhanced  bioremediation, 
monitoring  strategies,  and  their  interactions  should  enable  implementation  of  a  remedial  action 
plan  that  has  a  high  probability  of  success  with  the  least  cost. 
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_ Chapter  9 _ 

Application  to  Dover  AFB  Area  5 


9.1  Site  Description 

Dover  Air  Force  Base  (AFB)  is  located  in  Kent  County,  Delaware,  and  has  been  in  operation 
since  1942.  Base  operations  have  generated  numerous  wastes,  including  solvents  and 
hydrocarbons,  which  were  historically  buried  in  drums  or  disposed  in  the  storm  drainage  system. 
Wastes  were  disposed  in  various  on-base  locations,  which  have  resulted  in  several  NAPL  sources 
and  contaminant  plumes.  Our  focus  in  this  study  is  a  merged  groundwater  plume  associated  with 
several  sources  within  the  West  Management  Unit  (WMU)  designated  as  Area  5  (ORNL,  2008). 

Dover  AFB  is  underlain  by  unconsolidated  sediments  of  the  Atlantic  Coastal  Plain  (USGS, 
2000).  The  units  of  interest  are  the  Columbia  Formation  and  the  underlying  Calvert  Formation. 
The  Columbia  Formation  is  mainly  composed  of  sands,  silts  and  gravels.  A  clayey  silt  unit 
separates  it  from  the  Frederica  aquifer  in  the  upper  Calvert  Formation.  The  water  table 
configuration  at  Dover  AFB  is  generally  controlled  by  surface  water  bodies  including  Little 
River,  St.  Jones  River,  and  Delaware  Bay  (USGS,  2000)  and  recharge  from  precipitation.  The 
hydraulic  conductivity  is  variable  ranging  from  0.03  m/day  to  76  m/day  and  the  hydraulic 
gradient  varies  over  three  orders  of  magnitude  with  groundwater  velocity  ranging  from  about 
3x1 0'5  m/day  to  0.08  m/day. 

The  five  major  known  contaminant  sources  in  Area  5  are  designated  as  OT51,  OT50,  SS20, 
OT44,  and  OT41  (Figure  9-1).  They  were  fonnerly  oil/water  (O/W)  separator  sites  (OT51, 
OT50,  OT44),  a  vehicle  refueling  site  with  underground  storage  tanks  (UST)  (SS20),  and  the 
Industrial  Waste  Collection  Drain  (IWCD)  site  (OT41). 


Figure  9-1.  Contaminant 
plume,  streamlines,  and 
locations  of  sources  and 
EVO  injection  galleries 
at  Dover  AFB  Area  5. 
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Chlorinated  ethenes,  ranging  from  VC  to  PCE  and  soluble  fuel  hydrocarbons  (BTEX)  occur 
throughout  the  area  (Table  9-1).  A  number  of  remedial  actions  have  been  undertaken,  the  most 
important  of  which  are  summarized  in  Table  9-2.  Currently,  injection  of  EVO  in  10  injection 
well  galleries  (Figure  9-1)  is  being  employed  to  enhance  microbial  reductive  dechlorination  of 
chlorinated  ethenes.  The  system  is  referred  to  as  Accelerated  Anaerobic  Bioremediation  with 
Pennanent  Injection  Circulation  Transects  (AAB  PICT). 


Table  9-1.  Detected  contaminants  (ORNL,  2008). 


Contaminant 

OT51 

OT50 

SS20 

OT44 

OT41 

Benzne 

X 

X 

X 

X 

cis-l,2-DCE 

X 

X 

X 

• 

X 

PCE 

X 

X 

X 

X 

TCE 

X 

X 

X 

X 

X 

Toluene 

X 

X 

X 

VC 

X 

X 

• 

• 

• 

X  Contaminant  present  at  this  source. 

•  Potential  contaminants  of  concern  (COC)  due  to  the  breakdown  of  other  COCs. 


Table  9-2.  Chronology  of  Events  (ORNL,  2008). 


Date 

Event 

April  2006 

Record  of  decision  (ROD)  prepared  for  OT51,  OT50,  SS20,  OT44, 
OT41/Motor  Pool,  and  Area  5  (WMU) 

May  2006 

Finalized  Area  5  remedial  action  work  plan  (RAWP) 

May  8  -  24,  2006 

Collected  membrane  interface  probe  data  delineate  Source  Areas  D 
and  E  prior  to  treatment  injections 

May  1 6  -  June  9,  2006 

AAB  direct  injection  performed  at  Source  Areas  E,  F,  G,  and  H 

May  -  July  2006 

Installation  of  AAB  and  monitoring  wells 

May  -  October  2006 

First  round  of  AAB  PICT  injections 

Week  of  June  12,  2006 

Baseline  groundwater  sampling  conducted  in  upgradient  portion  of 
Area  5  (OT51/OT50) 

June  2006  -  on  going 

AAB  groundwater  monitoring 

July  2006 

Natural  Attenuation  monitoring  conducted  for  Area  5 

January/February  2007 

Natural  Attenuation  monitoring  conducted  for  Area  5 

July  2007  -  November  2007 

Second  round  of  AAB  PICT  injections;  stopped  before  completion  b 
ecause  of  cold  weather 

9.2  Model  Calibration 

Groundwater  flow  and  contaminant  transport  is  approximated  by  the  two-dimensional  depth- 
averaged  model  described  in  Chapter  2.  Groundwater  streamlines  originating  at  each 
contaminant  source  were  estimated  from  water  table  contours  and  the  contaminant  plume  and 
described  by  simple  polynomial  regressions  (Figure  9-1).  Streamline  functions  starting  from  each 
EVO  injection  gallery  were  also  determined  (not  shown  in  figure). 
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The  forward  model  at  Area  5  simulates  DNAPL  source  dissolution  and  transport  for  each  source 
as  well  as  transport  of  injected  ED  from  each  EVO  injection  gallery.  Since  the  model  does  not 
explicitly  simulate  incomplete  dechlorination,  we  model  total  chlorinated  ethenes  as  a  pseudo¬ 
species.  Since  TCE  is  the  primary  contaminant,  we  sum  the  concentration  of  all  chlorinated 
ethene  species  after  converting  to  their  stochiometric  equivalent  quantity  of  TCE  as  described  in 
Chapter  3.  Model  calibration  was  performed  initially  using  chlorinated  solvent  data  collected 
from  1994  to  2006,  before  ED  injection  started. 

For  each  contaminant  source,  we  calibrated  five  parameters:  contaminant  mass  remaining  in  the 
source  in  2006,  source  discharge  rate  in  2006,  source  width,  initial  release  date,  and  mass 
depletion  exponent.  To  reduce  the  number  of  parameters,  we  assume  that  all  sources  have  the 
same  depletion  exponent.  Besides  parameters  related  to  sources,  there  are  also  eight  aquifer 
parameters  of  which  porosity  (0.3),  natural  ED  concentration  (1200  H-eq  ppb)  and  natural  EA 
concentration  (2400  H-eq  ppb)  levels  were  assumed  known.  Prior  distributions  of  all  parameters 
are  represented  as  lognormal  except  for  the  source  release  date,  which  is  represented  as  nonnal. 
Prior  best  estimates  of  all  parameters  and  their  standard  deviations  are  listed  in  Table  9-3. 


Table  9-3.  Parameter  prior  estimates  and  their  uncertainty. 


Parameters 

Prior  value1 

STD2 

Reference 

Site  OT51 

Mass  in  2006  (kg) 

3000 

1.0 

ORNL,  2008 

(Source  1) 

Source  flux  in  2006  (kg/d) 

0.01 

1.0 

ORNL,  2008 

Release  date  (y) 

1975 

3.0 

ORNL,  2008 

Width  (m) 

50 

0.20 

ORNL,  2008 

Depletion  exponent  (-) 

1.3 

0.10 

Site  OT50 

Mass  in  2006  (kg) 

3000 

1.0 

ORNL,  2008 

(Source  2) 

Source  flux  in  2006  (kg/d) 

0.01 

1.0 

ORNL,  2008 

Release  date  (y) 

1975 

3.0 

ORNL,  2008 

Width  (m) 

50 

0.20 

ORNL,  2008 

Depletion  exponent  (-) 

1.3 

0.10 

Site  SS20 

Mass  in  2006  (kg) 

3000 

1.0 

ORNL,  2008 

(Source  3) 

Source  flux  in  2006  (kg/d) 

0.01 

1.0 

ORNL,  2008 

Release  date  (y) 

1975 

3.0 

ORNL,  2008 

Width  (m) 

50 

0.20 

ORNL,  2008 

Depletion  exponent  (-) 

1.3 

0.10 

Site  OT41 

Mass  in  2006  (kg) 

3000 

1.0 

ORNL,  2008 

(Source  4) 

Source  flux  in  2006  (kg/d) 

0.01 

1.0 

ORNL,  2008 

Release  Date  (y) 

1975 

3.0 

ORNL,  2008 

Width  (m) 

50 

0.20 

ORNL,  2008 

Depletion  exponent  (-) 

1.3 

0.10 

Site  OT44 

Mass  in  2006  (kg) 

3000 

1.0 

ORNL,  2008 

(Source  5) 

Source  flux  in  2006  (kg/d) 

0.01 

1.0 

ORNL,  2008 

Release  date  (y) 

1975 

3.0 

ORNL,  2008 

Width  (m) 

50 

0.20 

ORNL,  2008 

Depletion  exponent  (-) 

1.3 

0.10 

Aquifer 

Darcy  velocity  (m/d) 

0.03 

0.10 

USGS,  2000 

Retardation  factor  (-) 

1.10 

0.56 

Typical  from  literature 

Longitudinal  dispersivity  (m) 

20 

0.20 

Estimated  from  plume  length 

Transverse  dispersivity  (m) 

2 

0.20 

Estimated  from  plume  length 

Saturated  thickness  (m) 

12 

0.20 

USGS,  2000 

1  Prior  values  are  initial  best  estimates  representing  the  geometric  mean  for  all  variables  except  release  date  and  source-depletion 


exponent,  which  are  represented  by  the  arithmetic  mean. 

2  Standard  deviations  are  for  natural  logarithms  of  all  variables  except  release  date  and  source-depletion  exponent,  which  are  the 
standard  error  of  actual  values. 
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For  lognormally  distributed  parameters,  prior  best  estimates  represent  the  geometric  mean  of  the 
variable  and  the  standard  deviation  is  the  statistic  for  the  natural  logarithm  (S/„).  Note  that  the 
99%  confidence  limit  for  a  log-normally  distributed  variable  ranges  from  approximately  exp(- 
3 Si„)  to  exp(3 Si„)  times  the  geometric  mean.  For  example,  the  prior  estimate  of  Source  1 
contaminant  mass  in  2006  has  a  geometric  mean  of  3000  kg  with  a  standard  deviation  of 
ln(mass)  equal  to  1  (dimensionless)  indicating  a  99%  confidence  range  of  -0.05  to  20  times  3000 
kg  or  150  to  60,000  kg. 

The  calibration  approach  described  in  Chapter  5  implemented  in  SCOToolkit  was  used  to 
estimate  model  parameters  and  their  uncertainties  as  well  as  the  measurement  error  in 
concentration  data.  Best  estimates  of  parameters  as  well  as  their  posterior  standard  deviations  are 
summarized  in  Table  9-4  and  the  data  fit  is  illustrated  in  Figure  9-2. 


Table  9-4.  Summary  of  calibration  results  for  Dover  AFB  Area  5. 


Parameters 

Best  Estimates 

STD1 

Mass  at  source  1  in  2006  (kg) 

3065 

1.0 

Mass  at  source  2  in  2006  (kg) 

3224 

1.0 

Mass  at  source  3  in  2006  (kg) 

3069 

1.0 

Mass  at  source  4  in  2006  (kg) 

3017 

1.0 

Mass  at  source  5  in  2006  (kg) 

3010 

1.0 

Total  Mass 

15385 

0.45 

Flux  at  source  1  in  2006  (kg/d) 

0.006 

0.58 

Flux  at  source  2  in  2006  (kg/d) 

0.010 

0.53 

Flux  at  source  3  in  2006  (kg/d) 

0.006 

0.61 

Flux  at  source  4  in  2006  (kg/d) 

0.002 

0.90 

Flux  at  source  5  in  2006  (kg/d) 

0.001 

0.73 

Total  Flux 

0.025 

0.70 

Release  date  of  source  1  (y) 

1974 

2.9 

Release  date  of  source  2  (y) 

1985 

2 

Release  date  of  source  3  (y) 

1977 

2.5 

Release  date  of  source  4  (y) 

1975 

3.0 

Release  date  of  source  5  (y) 

1975 

3.0 

Width  of  source  1  (m) 

45 

0.19 

Width  of  source  2  (m) 

51 

0.20 

Width  of  source  3  (m) 

81 

0.15 

Width  of  source  4  (m) 

48 

0.20 

Width  of  source  5  (m) 

51 

0.20 

Source  depletion  exponent  (-) 

1.3 

0.1 

Darcy  velocity  (m/d) 

0.015 

0.056 

Retardation  factor 

1.10 

0.56 

Longitudinal  dispersivity  (m) 

11.3 

0.18 

Transverse  dispersivity  (m) 

1.6 

0.15 

Saturated  thickness  (in) 

15.8 

0.20 

1  Standard  deviations  are  for  natural  logarithms  of  all  variables  except  release  date  and  source  depletion 
exponent,  which  are  the  standard  error  of  the  actual  values. 

"  The  estimated  value  is  the  upper  limit  constraint  set  in  optimization,  so  uncertainty  estimate  from  linear 
analysis  is  not  available. 
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Darcy  velocity  and  release  dates  are  reasonably  well  characterized  by  the  calibration.  Estimates 
of  source  mass  and  its  uncertainty  barely  changed  from  the  prior  estimates,  indicating  that  the 
available  data  is  insufficient  to  reliably  estimate  source  mass.  Note  that  99%  confidence  limits 
for  the  mass  of  each  source  in  2006  range  from  about  150  to  60,000  kg.  Other  parameters  exhibit 
intermediate  reliability. 
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Figure  9-2.  Observed  versus  calibrated  contaminant  concentrations. 

9.3  Design  Optimization 

Given  the  large  uncertainty  in  source  mass  estimates,  reliable  prediction  of  the  duration  of 
DNAPL  dissolution  that  will  feed  the  dissolved  plume  is  not  expected.  However,  due  to  the  low 
estimated  DNAPL  dissolution  rates  relative  to  source  mass,  forward  simulations  indicate  a  high 
probability  that  without  source  remediation,  ED  injection  or  some  other  means  of  plume  control 
will  be  required  for  a  very  long  time  (decades  to  centuries).  Furthermore,  due  to  the  number  of 
known  DNAPL  sources,  a  possibility  of  additional  smaller  source  zones,  and  difficulty  to 
perform  detailed  source  zone  characterization  or  remediation  due  to  buildings,  infrastructure  and 
base  operation  disruptions,  we  will  focus  on  investigating  remedial  options  capable  of 
minimizing  the  cost  of  long-temi  plume  control  using  EVO  injection  with  ten  of  the  ED  injection 
galleries  that  have  already  been  constructed  within  Area  5. 

Specifically,  our  goal  is  to  minimize  total  costs  over  a  100  year  period  subject  to  the  constraint 
that  after  a  grace  period  the  TCE-equivalent  concentration  at  each  of  the  five  compliance  wells, 
signified  by  the  “@”  sign  on  Figure  9-1,  should  be  less  than  MCL  for  each  MC  realization. 
Formally,  we  minimize  the  expected  present  value  of  the  total  cost  over  100  years.  Note  that  this 
criterion  has  a  conservative  bias  since  other  chlorinated  ethane  species  present  at  levels 
comparable  to  TCE  have  lower  MCLs.  One  hundred  realizations,  conditioned  on  the  actual  data 
used  for  calibration,  are  used  in  the  optimization.  By  enforcing  that  all  100  realizations  need  to 
meet  MCL  constraints  at  compliance  wells,  we  implicitly  ensure  that  the  probability  of  failure  of 
the  optimal  remediation  design  is  on  the  order  of  1%.  The  grace  period  is  computed  as  the  travel 
time  of  ED  from  the  nearest  injection  gallery  to  a  compliance  well  plus  1  year.  The  MCL  value  is 
5  TCE-equivalent  pg/L.  Cost  variables  for  Dover  site  are  summarized  in  Table  9-5. 
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Table  9-5.  Summary  of  cost  variables  used  in  optimization. 


Description 

Value 

Capital  Costs 

Gallery  construction 

Injection  equipment 

$333/  gallery  width  (in) 
$65,000/  unit1 

Operating  Costs 

Fork  lift  rental  and  operation 
Generator  rental  and  operation 
Personnel 

$21,600/  year 
$21,600/  year/unit 
$65,000/year/gallery 

Material  Cost 

ED  material  (in  solution) 

$1.3  /  kg  solution1 

Annual  Discount  Rate 

Risk-free  return  minus  inflation  rate 

3% 

1  The  H-equiv  mass  per  kg  of  solution  is  approximately  0.063kg 


Optimization  under  uncertainty  requires  a  significant  amount  of  computational  resources.  In  this 
application,  where  we  only  consider  plume  control  with  constant  ED  injection  over  a  fixed  time 
period  and  without  source  zone  remediation,  the  ED  concentration  is  a  linear  function  of 

injection  rates  at  each  gallery;  thus,  the  contaminant  concentration  can  be  calculated  easily  by 
taking  advantage  of  the  linear  relationship  between  injection  rates  and  ED  concentrations.  Once 
we  calculate  the  DNAPL  concentration  without  ED  injection,  and  the  ED  concentration  due  to 
unit  ED  injection  rate  (i.e.,  1  kg-solution/day)  from  each  injection  gallery  over  the  simulation 
period,  MC  simulations  can  be  replaced  by  a  simple  linear  equation.  Furthermore,  with  the  linear 
relationship  between  DNAPL  concentration  and  ED  injection  rate,  we  can  find  the  optimal  ED 
injection  strategy  for  Dover  site  using  linear  programming.  More  precisely,  mixed  binary  integer 
programming  was  used  in  this  study  to  determine  whether  one  gallery  is  used  or  not  and  the 
amount  of  ED  injected  in  case  it  is  used. 

We  present  optimization  results  for  the  following  optimization  scenarios: 

•  Optimize  ED  injection  rates  in  the  ten  existing  ED  injection  galleries  at  the  site  over  a 
100  year  time  period  to  minimize  overall  remediation  costs,  and 

•  Optimize  ED  injection  rates  in  2011  and  implement  the  optimal  strategy  for  50  years. 
In  2061,  model  parameters  are  recalibrated  using  additional  data  since  2011.  Then,  with 
the  updated  calibration,  reoptimize  ED  injection  rates  for  the  next  50  years. 

9.3.1  ED  injection  optimization  for  100  years 

ED  is  injected  at  a  constant  rate  in  each  gallery  for  100  years.  Rates  can  differ  from  gallery  to 
gallery.  The  optimal  strategy  for  the  operating  galleries  are  shown  in  Table  9-6.  Expected  values 
for  capital  cost,  cumulative  NPV  operating  cost,  total  NPV  cost,  and  annual  operating  cost  are 
summarized  in  Table  9-7.  Confidence  limits  for  dissolved  TCE-equivalent  concentrations  at 
compliance  well  DM329S  (Figure  9-3)  indicate  that  by  2030,  there  is  a  high  likelihood  of 
meeting  the  cleanup  level.  The  predicted  spatial  distribution  of  groundwater  contamination  in 
2111  using  best  estimates  of  model  parameters  and  optimized  design  variables  (Figure  9-4) 
shows  that  EVO  injection  is  reasonably  successful  at  containing  the  dissolved  plume,  although 
improvements  may  be  possible  if  gallery  locations  were  optimized. 
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Table  9-6.  Single-stage  optimization  of  ED  injection  rates  for  Dover  AFB  Area  5. 


ED  Injection  Gallery 

Injection  Rate  (kg-sol/yr) 

1 

0 

2 

0 

3 

620 

4 

0 

5 

365 

6 

219 

7 

0 

8 

219 

9 

0 

10 

365 

All 

1788 

Table  9-7.  Expected  cost  details  for  single-stage  ED  optimization. 


Capital  Cost 

Operating  Cost 

Annual  Operating  Cost 

Total  Cost 

(ENPV  $K) 

(ENPV  $K) 

($K/yr) 

(ENPV  $K) 

431 

17,415 

549 

17,846 

Year 


Figure  9-3.  TCE-equivalent  concentration  at  a  compliance  point  DM329S. 
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Figure  9-4.  Contour  plot  of  TCE-equivalent  concentration  in  21 1 1  with  best  parameter  estimates 
(red  dots  and  numbers  represent  operating  gallery  locations 
and  corresponding  ED  injection  rates). 


9.3.2  ED  injection  multi-stage  optimization  with  parameter  recalibration 

In  this  section,  we  consider  two-stage  optimization  for  ED  injection  rates.  Injection  rates 
determined  from  the  initial  calibration/optimization  are  assumed  to  be  implemented  in  2011  for 
50  years.  In  2061,  model  parameters  are  re-calibrated  using  the  additional  observed 
concentrations  in  monitoring  wells  since  initial  implementation.  Then,  with  new  model 
parameters,  ED  injection  rates  are  re-optimized  and  applied  for  the  next  50  years.  In  this  way, 
the  remediation  strategy  can  be  adjusted  based  on  field  observations.  This  observation-feedback 
approach  should  reduce  the  probability  of  failure  and  expected  costs.  To  obtain  concentration 
“observations”  for  the  first  50  years,  "true"  model  parameters  are  selected  from  a  parameter 
realization  generated  from  the  original  calibration  which  is  not  used  in  optimization. 

The  results  show  that  after  recalibration  and  optimization,  ED  galleries  6,  8,  and  10  are  shut 
down,  and  the  total  injection  rate  is  reduced  from  1788  kg-sol/yr  to  876  kg-sol/yr  (Table  9-8), 
resulting  in  a  reduction  in  annual  (undiscounted)  operating  costs  from  $549  K/yr  to  $207  K/yr 
(Tables  9-7  and  9-9).  It  is  also  observed  that  after  reoptimization,  injection  rates  in  injection 
galleries  3  and  5  were  adjusted. 
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Table  9-8.  Two-stage  optimization  of  ED  injection  rates  for  Dover  AFB  Area  5. 


ED  Injection 

Injection  Rate  (kg-sol/yr) 

Gallery 

2011-2060 

2061-2111 

1 

0 

0 

2 

0 

0 

3 

620 

803 

4 

0 

0 

5 

365 

73 

6 

219 

0 

7 

0 

0 

8 

219 

0 

9 

0 

0 

10 

365 

0 

All 

1788 

876 

Table  9-9.  Cost  details  for  two-stage  ED  optimization. 


Capital  Cost 

Operating  Cost 

Annual  Operating  Cost1 

Total  Cost 

(NPV  $K) 

(NPV  $K) 

($K/yr) 

(NPV  $K) 

431 

15,484 

207 

15,915 

1  Annual  operating  costs  for  years  2011-2060  /  2061-2110 


The  spatial  distribution  of  groundwater  contamination  in  2061  indicates  that  some  contamination 
moves  beyond  the  most  downgradient  monitoring  wells  before  being  attenuated  by  EVO 
injection  (Figure  9-5).  Given  the  paucity  of  data,  this  is  difficult  to  confirm  or  refute.  With 
additional  downgradient  data,  the  model  accuracy  could  be  evaluated  and  the  calibration  refined. 

Comparison  of  the  two-stage  optimization  plume  simulation  in  2 1 1 1  (Figure  9-6)  with  the  one- 
stage  optimization  (Figure  9-4)  reveals  that  although  both  cases  yield  high  probabilities  of 
meeting  the  MCE  at  the  selected  compliance  locations,  fingers  of  contamination  may  slip  by 
compliance  wells.  This  risk  can  be  ameliorated  by  visual  inspection  of  the  plume  distribution  and 
by  using  the  optimization  to  select  compliance  well  locations. 
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Figure  9-5.  Contour  plot  of  TCE-equivalent  concentration  in  2061  using  the  best  parameter 
estimate  after  recalibration  in  2061  (red  dots  and  numbers  represent  operating 
gallery  locations  and  corresponding  ED  injection  rates). 
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Figure  9-6.  Contour  plot  of  TCE-equivalent  concentration  in  21 1 1  using  the  best  parameter 
estimate  after  recabbration  in  2061  (red  dots  and  numbers  represent  operating 
gallery  locations  and  corresponding  ED  injection  rates  assigned  in  2061). 


9-10 


9.4  Conclusions 


Estimates  of  source  mass  could  not  be  calibrated  reliably  from  available  site  data.  Ninety-nine 
percent  confidence  limits  for  the  mass  of  the  five  identified  DNAPL  sources  in  Dover  AFB  Area 
5  ranged  from  about  150  to  60,000  kg.  Forward  simulations  indicate  a  high  probability  that 
without  source  remediation,  long-term  measures  are  likely  to  be  required  to  control  off-base 
plume  migration.  Due  to  the  number  of  known  DNAPF  sources,  the  possibility  of  additional 
smaller  source  zones,  and  the  difficulty  to  perform  detailed  source  zone  characterization  or 
remediation  due  to  infrastructure  and  base  operation  disruptions,  plume  control  using  EVO 
injection  may  be  a  practical  long-term  solution. 

Optimization  of  EVO  injection  rates  in  ten  existing  injection  galleries  indicated  a  high 
probability  that  contaminant  concentrations  in  specified  compliance  wells  could  meet  cleanup 
targets  with  continued  injection  in  only  5  of  the  current  10  injection  galleries,  which  would 
reduce  operating  costs  by  about  50%. 

A  multi-stage  calibration-optimization  strategy  was  investigated  in  which  indicated  that  periodic 
re-calibration  and  re-optimization  using  additional  monitoring  data  can  further  reduce  expected 
costs  and  improve  system  performance  and  the  reliability.  Re-optimizing  after  50  years  was 
predicted  to  allow  a  reduction  in  ED  injection  rates  that  reduced  (non-discounted)  annual 
operating  costs  by  an  additional  60%. 

A  possible  strategy  that  may  be  investigated  in  future  studies  to  reduce  source  mass  and  model 
uncertainty  would  be  to  inject  EVO  immediately  up  gradient  of  each  of  source  zones  for  several 
injection  cycles.  High  EVO  concentrations  in  the  source  zones  is  expected  to  increase  DNAPL 
dissolution  rates  and  thus  accelerate  source  mass  depletion.  Recalibration  using  monitoring  data 
obtained  during  accelerated  source  depletion  should  diminish  model  uncertainty  and  enable  more 
reliable  system  design  and  perfmroance. 
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_ Chapter  10 _ 

Conclusions,  Future  Research  and  Implementation 


A  new  approach  for  optimizing  the  design  and  operation  of  DNAPL  site  remediation  systems 
was  developed  and  implemented  in  the  computer  program  SCOToolkit,  which  explicitly 
considers  uncertainty  in  site  and  remediation  system  characteristics,  measurement  data, 
performance  models  and  unit  costs,  and  specific  compliance  criteria.  By  seeking  to  minimize 
NPV  cost  averaged  over  equiprobable  Monte  Carlo  realizations  of  future  performance, 
conditioned  on  all  available  site  data,  with  penalties  applied  for  nonattainment  of  remediation 
objectives,  the  approach  constrains  the  probability  of  failure.  Unlike  conventional  design  and 
operational  optimization  protocols,  the  proposed  approach  takes  into  account  the  many  complex 
interactions  and  tradeoffs  among  design,  monitoring  and  characterization  variables  and  cleanup 
criteria.  For  example,  remediation  optimization  often  attempts  to  minimize  sampling  frequency, 
while  SCOToolkit  can  determine  an  optimum  sampling  frequency  at  which  the  increased  cost  of 
additional  samples  is  just  offset  by  cost  savings  associated  with  improved  decision  reliability 
(i.e.,  the  magnitude  of  overdesign  is  adjusted  to  minimize  expected  cost). 

The  method  was  applied  to  two  well-characterized  sites  where  different  remedial  technologies 
were  used  to  evaluate  its  ability  to  reduce  costs  and  improve  remedial  designs.  At  the  first  site, 
the  Fort  Lewis  East  Gate  Disposal  Yard  (EGDY)  site,  optimization  of  thermal  source  treatment 
indicated  a  need  for  a  much  larger  treatment  area  than  was  actually  employed  to  avoid  a  high 
failure  probability  associated  with  source  delineation  uncertainty  based  on  available  source 
characterization  data.  The  method  was  also  used  to  optimize  source  and  plume  bioremediation  at 
the  EGDY  site,  using  whey  injection  without  additional  source  reduction.  The  results  indicated 
that  this  remedial  strategy  is  expected  to  achieve  Maximum  Contaminant  Limits  by  2100  with  a 
93%  probability  of  success  when  using  relatively  low  whey  injection  rates.  Additional  source 
treatment  efforts  may  be  cost  effective  if  additional  characterization  were  undertaken  to  reduce 
source  zone  uncertainty.  Further  studies  may  be  undertaken  using  SCOToolkit  to  evaluate  the 
potential  monetary  value  of  additional  source  zone  investigations  in  conjunction  with  a  second 
round  of  thermal  or  other  source  zone  treatment. 

SCOToolkit  was  also  applied  to  Dover  AFB  Area  5  to  optimize  enhanced  bioremediation.  The 
results  suggested  that  electron  donor  injection  may  be  needed  indefinitely  in  the  future  to  meet 
remediation  criteria,  as  a  result  of  the  large  uncertainty  in  estimated  source  mass  and  difficulties 
in  characterizing  and  treating  source  zones  without  disrupting  base  operations.  Optimization 
analyses  to  minimize  long  term  operating  costs  indicated  a  high  probability  of  meeting 
compliance  criteria  using  only  five  of  the  current  ten  emulsified  vegetable  oil  injection  galleries, 
with  approximately  half  the  current  annual  operating  cost.  Recalibration  and  optimization  after 
50  years,  using  additional  data  from  this  period,  was  projected  to  reduce  prediction  uncertainty 
allowing  further  reductions  in  operating  costs.  Recommendations  for  further  optimization 
analyses  include  investigating  EVO  injection  immediately  upgradient  instead  of  downgradient  of 
each  source  zone  to  increase  DNAPL  dissolution  rates,  reduce  source  mass  more  quickly,  and 
recalibrate  at  shorter  time  intervals  to  accelerate  operating  cost  reductions. 

Stochastic  cost  optimization  of  design  and  monitoring  variables  was  shown  to  reduce  expected 
costs  by  50%  or  more  compared  to  conventional  design  methods,  and  to  substantially  increase 
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the  probability  of  meeting  compliance  targets.  Expected  cost  reductions  were  also  shown  to  be 
possible  by  judicious  specification  of  statistical  criteria  for  no  further  action  (NFA)  and  by 
stipulating  procedures  for  technology  deployment,  operation  and  shut-down.  For  example, 
source  zone  soil  and/or  groundwater  monitoring  may  be  used  to  operationally  decide  when  to 
tenninate  electrode  heating  for  thermal  treatment  all  at  once,  within  zones,  or  individually.  The 
optimal  protocol  depends  on  tradeoffs  between  the  magnitude  of  measurement  noise,  unit  costs 
for  well  installation  and  monitoring,  and  other  variables.  Over  many  sites,  minimization  of 
expected  cost  should  yield  significant  cost  savings.  For  any  single  site,  the  probability  of  failure 
is  reduced  to  a  level  that  the  probability- weighted  cost  of  failure  is  just  offset  by  the  cost  of 
overdesigning  the  system. 

Various  enhancements  to  SCOToolkit  could  be  implemented  to  extend  applicability  to  a  broader 
range  of  remedial  technologies  and  to  make  the  program  easier  to  use.  With  limited  effort,  the 
DNAPF  source  remediation  module  could  be  extended  to  add  the  capability  to  simulate  in  situ 
chemical  oxidation  (ISCO),  thermal  conduction  and  steam  injection  perfonnance  and  cost,  in 
addition  to  the  current  electrical  resistance  heating  capability.  Pre-  and  post-processors  would 
simplify  problem  setup,  execution  and  interpretation  of  output. 

In  addition  to  more  extensive  technology  demonstrations  on  the  Dover  AFB  and  Fort  Fewis 
sites,  as  discussed  above,  further  testing  would  be  beneficial  to  demonstrate  the  applicability  and 
benefits  of  SCOToolkit  on  DoD  sites  that  span  ranges  of  site  complexity,  data  availability,  and 
remediation  technologies,  including  enhanced  source  zone  and  plume  bioremediation,  various 
thermal  technologies,  and  ISCO,  used  individually  or  in  various  combinations.  In  addition  to 
single-pass  calibration-optimization  applications,  demonstrations  of  iterative  calibration- 
optimization  protocols  are  desirable  to  assess  potential  cost  reductions  due  to  additional  site 
characterization  efforts. 

Finally,  technology  transfer  will  be  critical  to  fully  realize  cost  savings  promised  by  SCOToolkit. 
Responsible  parties  and  regulators  must  gain  an  understanding  of  the  technology  and  its 
capabilities,  and  training  must  be  available  for  consultants  and  other  users  to  learn  how  to 
effectively  and  efficiently  use  the  software. 
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Appendix  A.  Concentration  and  Source  Data  for  Fort  Lewis  Calibration 


Table  A-l.  TCE-equivalent  concentrations  used  for  Fort  Lewis  site  calibration. 


No. 

Well 

East  (m) 

North  (m) 

Date 

Concentration  (ppb) 

Aquifer 

1 

LC-03 

455341.94 

200345.95 

06/01/1995 

0.80 

1 

2 

LC-03 

455341.94 

200345.95 

06/01/1996 

0.50 

1 

3 

LC-03 

455341.94 

200345.95 

06/01/1997 

0.60 

1 

4 

LC-03 

455341.94 

200345.95 

06/01/1998 

1.10 

1 

5 

LC-03 

455341.94 

200345.95 

06/01/1999 

1.00 

1 

6 

LC-03 

455341.94 

200345.95 

06/01/2000 

5.40 

1 

7 

LC-03 

455341.94 

200345.95 

06/01/2001 

1.60 

1 

8 

LC-05 

454413.21 

200342.91 

06/01/1995 

46.50 

1 

9 

LC-05 

454413.21 

200342.91 

06/01/1996 

22.40 

1 

10 

LC-05 

454413.21 

200342.91 

06/01/1997 

26.50 

1 

11 

LC-05 

454413.21 

200342.91 

06/01/1998 

29.00 

1 

12 

LC-05 

454413.21 

200342.91 

06/01/1999 

26.30 

1 

13 

LC-05 

454413.21 

200342.91 

06/01/2000 

39.80 

1 

14 

LC-05 

454413.21 

200342.91 

06/01/2001 

65.70 

1 

15 

LC-06 

455369.37 

199917.10 

06/01/1995 

88.50 

1 

16 

LC-06 

455369.37 

199917.10 

06/01/1996 

46.00 

1 

17 

LC-06 

455369.37 

199917.10 

06/01/1997 

33.40 

1 

18 

LC-06 

455369.37 

199917.10 

06/01/1998 

61.30 

1 

19 

LC-06 

455369.37 

199917.10 

06/01/1999 

72.50 

1 

20 

LC-06 

455369.37 

199917.10 

06/01/2000 

88.50 

1 

21 

LC-06 

455369.37 

199917.10 

06/01/2001 

67.30 

1 

22 

LC-108 

456129.11 

198922.98 

06/01/1995 

151.00 

1 

23 

LC-108 

456129.11 

198922.98 

06/01/1996 

15.30 

1 

24 

LC-108 

456129.11 

198922.98 

06/01/1997 

22.80 

1 

25 

LC-108 

456129.11 

198922.98 

06/01/1998 

50.50 

1 

26 

LC-108 

456129.11 

198922.98 

06/01/1999 

23.70 

1 

27 

LC-108 

456129.11 

198922.98 

06/01/2000 

14.30 

1 

28 

LC-108 

456129.11 

198922.98 

06/01/2001 

11.00 

1 

29 

LC-132 

454582.08 

200260.81 

06/01/1995 

23.00 

1 

30 

LC-132 

454582.08 

200260.81 

06/01/1996 

52.30 

1 

31 

LC-132 

454582.08 

200260.81 

06/01/1997 

59.80 

1 

32 

LC-132 

454582.08 

200260.81 

06/01/1998 

62.00 

1 

33 

LC-132 

454582.08 

200260.81 

06/01/1999 

79.00 

1 

34 

LC-132 

454582.08 

200260.81 

06/01/2000 

93.50 

1 

35 

LC-132 

454582.08 

200260.81 

06/01/2001 

102.00 

1 

36 

LC-137A 

456031.98 

198938.29 

06/01/1995 

191.43 

1 

37 

LC-137A 

456031.98 

198938.29 

06/01/1996 

93.47 

1 

38 

LC-137A 

456031.98 

198938.29 

06/01/1997 

78.20 

1 

39 

LC-137A 

456031.98 

198938.29 

06/01/1998 

132.27 

1 

40 

LC-137A 

456031.98 

198938.29 

06/01/1999 

86.70 

1 

41 

LC-137A 

456031.98 

198938.29 

06/01/2000 

122.17 

1 

42 

LC-137A 

456031.98 

198938.29 

06/01/2001 

205.33 

1 

43 

LC-144B 

455799.22 

199047.83 

06/01/1995 

195.25 

1 

44 

LC-144B 

455799.22 

199047.83 

06/01/1996 

148.40 

1 

45 

LC-144B 

455799.22 

199047.83 

06/01/1997 

136.15 

1 

46 

LC-144B 

455799.22 

199047.83 

06/01/1998 

97.00 

1 

47 

LC-16 

454905.47 

200648.62 

06/01/2001 

9.30 

1 

48 

LC-19A 

455718.37 

199063.36 

06/01/1998 

90.57 

1 

49 

LC-19A 

455718.37 

199063.36 

06/01/1999 

136.77 

1 

50 

LC-19A 

455718.37 

199063.36 

06/01/2000 

108.27 

1 
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Table  A-l.  TCE-equivalent  concentrations  used  for  Fort  Lewis  site  calibration  (continued) 


No. 

Well 

East  (m) 

North  (m) 

Date 

Concentration  (ppb) 

Aquifer 

51 

LC-19A 

455718.37 

199063.36 

06/01/2001 

163.30 

1 

52 

LC-41A 

454723.34 

199690.05 

06/01/1995 

165.00 

1 

53 

LC-41A 

454723.34 

199690.05 

06/01/1996 

170.00 

1 

54 

LC-41A 

454723.34 

199690.05 

06/01/1997 

187.50 

1 

55 

LC-41A 

454723.34 

199690.05 

06/01/1998 

160.00 

1 

56 

LC-41A 

454723.34 

199690.05 

06/01/1999 

167.50 

1 

57 

LC-41A 

454723.34 

199690.05 

06/01/2000 

162.50 

1 

58 

LC-41A 

454723.34 

199690.05 

06/01/2001 

193.30 

1 

59 

LC-44A 

455141.99 

200214.59 

06/01/1995 

45.50 

1 

60 

LC-44A 

455141.99 

200214.59 

06/01/1996 

30.00 

1 

61 

LC-44A 

455141.99 

200214.59 

06/01/1997 

16.80 

1 

62 

LC-44A 

455141.99 

200214.59 

06/01/1998 

16.00 

1 

63 

LC-44A 

455141.99 

200214.59 

06/01/1999 

24.50 

1 

64 

LC-44A 

455141.99 

200214.59 

06/01/2000 

28.30 

1 

65 

LC-44A 

455141.99 

200214.59 

06/01/2001 

30.70 

1 

66 

LC-49A 

455336.76 

199380.35 

06/01/1995 

133.00 

1 

67 

LC-49A 

455336.76 

199380.35 

06/01/1996 

155.90 

1 

68 

LC-49A 

455336.76 

199380.35 

06/01/1997 

166.15 

1 

69 

LC-49A 

455336.76 

199380.35 

06/01/1998 

173.00 

1 

70 

LC-51 

455784.81 

198661.63 

06/01/1995 

99.00 

1 

71 

LC-51 

455784.81 

198661.63 

06/01/1996 

130.00 

1 

72 

LC-51 

455784.81 

198661.63 

06/01/1997 

142.50 

1 

73 

LC-51 

455784.81 

198661.63 

06/01/1998 

162.50 

1 

74 

LC-51 

455784.81 

198661.63 

06/01/1999 

170.00 

1 

75 

LC-51 

455784.81 

198661.63 

06/01/2000 

162.50 

1 

76 

LC-51 

455784.81 

198661.63 

06/01/2001 

153.30 

1 

77 

LC-53 

455473.31 

198707.04 

06/01/1995 

115.00 

1 

78 

LC-53 

455473.31 

198707.04 

06/01/1996 

167.50 

1 

79 

LC-53 

455473.31 

198707.04 

06/01/1997 

155.00 

1 

80 

LC-53 

455473.31 

198707.04 

06/01/1998 

162.50 

1 

81 

LC-53 

455473.31 

198707.04 

06/01/1999 

195.00 

1 

82 

LC-53 

455473.31 

198707.04 

06/01/2000 

210.00 

1 

83 

LC-53 

455473.31 

198707.04 

06/01/2001 

200.00 

1 

84 

LC-66A 

454812.20 

200218.85 

06/01/1995 

76.75 

1 

85 

LC-66A 

454812.20 

200218.85 

06/01/1996 

112.00 

1 

86 

LC-66A 

454812.20 

200218.85 

06/01/1997 

115.90 

1 

87 

LC-66A 

454812.20 

200218.85 

06/01/1998 

111.75 

1 

88 

LC-66A 

454812.20 

200218.85 

06/01/1999 

119.15 

1 

89 

LC-66A 

454812.20 

200218.85 

06/01/2000 

106.65 

1 

90 

LC-66A 

454812.20 

200218.85 

06/01/2001 

91.20 

1 

91 

PA-381 

454330.00 

199657.72 

06/01/1995 

56.50 

1 

92 

PA-381 

454330.00 

199657.72 

06/01/1996 

35.80 

1 

93 

PA-381 

454330.00 

199657.72 

06/01/1997 

35.50 

1 

94 

PA-381 

454330.00 

199657.72 

06/01/1998 

39.30 

1 

95 

PA-381 

454330.00 

199657.72 

06/01/1999 

48.50 

1 

96 

PA-381 

454330.00 

199657.72 

06/01/2000 

47.80 

1 

97 

PA-381 

454330.00 

199657.72 

06/01/2001 

39.00 

1 

98 

LC-137B 

456035.56 

198940.36 

03/25/2004 

263.60 

1 

99 

LC-137B 

456035.56 

198940.36 

08/29/2005 

630.67 

1 

100 

LC-137B 

456035.56 

198940.36 

08/29/2006 

538.27 

1 
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Table  A-l.  TCE-equivalent  concentrations  used  for  Fort  Lewis  site  calibration  (continued) 


No. 

Well 

East  (m) 

North  (m) 

Date 

Concentration  (ppb) 

Aquifer 

101 

LC-137B 

456035.56 

198940.36 

08/31/2007 

326.51 

1 

102 

LC-137B 

456035.56 

198940.36 

09/02/2008 

727.02 

1 

103 

LC-137B 

456035.56 

198940.36 

08/31/2009 

269.33 

1 

104 

LC-137B 

456035.56 

198940.36 

03/02/2010 

322.48 

1 

105 

LC-160 

456044.89 

198862.04 

12/18/2006 

77.73 

1 

106 

LC-160 

456044.89 

198862.04 

12/12/2007 

42.51 

1 

107 

LC-160 

456044.89 

198862.04 

12/10/2008 

34.22 

1 

108 

LC-160 

456044.89 

198862.04 

12/07/2009 

22.11 

1 

109 

LC-160 

456044.89 

198862.04 

03/02/2010 

13.40 

1 

110 

LC-27 

456110.34 

198690.28 

09/28/2006 

28.53 

1 

111 

LC-27 

456110.34 

198690.28 

03/02/2007 

16.59 

1 

112 

LC-27 

456110.34 

198690.28 

03/19/2008 

8.70 

1 

113 

LC-27 

456110.34 

198690.28 

02/24/2009 

12.09 

1 

114 

LC-27 

456110.34 

198690.28 

03/02/2010 

11.01 

1 

115 

MT-1 

455984.84 

198929.11 

12/12/2006 

127.66 

1 

116 

MT-1 

455984.84 

198929.11 

12/12/2007 

106.31 

1 

117 

MT-1 

455984.84 

198929.11 

12/10/2008 

139.14 

1 

118 

MT-1 

455984.84 

198929.11 

12/07/2009 

41.07 

1 

119 

MT-1 

455984.84 

198929.11 

06/07/2010 

85.83 

1 

120 

MT-2 

456065.61 

198750.19 

12/12/2006 

15.75 

1 

121 

MT-2 

456065.61 

198750.19 

12/12/2007 

11.58 

1 

122 

MT-2 

456065.61 

198750.19 

12/10/2008 

13.02 

1 

123 

MT-2 

456065.61 

198750.19 

12/07/2009 

10.51 

1 

124 

MT-2 

456065.61 

198750.19 

06/07/2010 

7.51 

1 

125 

MT-3 

456103.63 

198642.71 

12/18/2006 

21.22 

1 

126 

MT-3 

456103.63 

198642.71 

12/12/2007 

13.65 

1 

127 

MT-3 

456103.63 

198642.71 

12/10/2008 

7.42 

1 

128 

MT-3 

456103.63 

198642.71 

12/07/2009 

14.17 

1 

129 

MT-3 

456103.63 

198642.71 

06/07/2010 

12.20 

1 

130 

MT-4 

456154.73 

198581.30 

12/18/2006 

26.30 

1 

131 

MT-4 

456154.73 

198581.30 

12/12/2007 

12.05 

1 

132 

MT-4 

456154.73 

198581.30 

12/10/2008 

4.68 

1 

133 

MT-4 

456154.73 

198581.30 

12/07/2009 

11.84 

1 

134 

MT-4 

456154.73 

198581.30 

06/07/2010 

4.57 

1 

135 

PW-3 

456057.00 

198897.68 

12/20/2006 

505.04 

1 

136 

PW-3 

456057.00 

198897.68 

12/18/2007 

265.24 

1 

137 

PW-3 

456057.00 

198897.68 

12/10/2008 

271.59 

1 

138 

PW-3 

456057.00 

198897.68 

12/07/2009 

176.84 

1 

139 

PW-4 

456069.12 

198840.78 

12/20/2006 

314.34 

1 

140 

PW-4 

456069.12 

198840.78 

12/18/2007 

162.52 

1 

141 

PW-4 

456069.12 

198840.78 

12/10/2008 

103.84 

1 

142 

PW-4 

456069.12 

198840.78 

12/07/2009 

54.26 

1 

143 

PW-4 

456069.12 

198840.78 

09/20/2010 

41.53 

1 

144 

PW-5 

456112.97 

198807.55 

12/20/2006 

118.60 

1 

145 

PW-5 

456112.97 

198807.55 

12/28/2007 

114.49 

1 

146 

PW-5 

456112.97 

198807.55 

12/10/2008 

57.65 

1 

147 

PW-5 

456112.97 

198807.55 

12/07/2009 

40.63 

1 

148 

PW-5 

456112.97 

198807.55 

09/20/2010 

54.72 

1 

149 

PW-6 

456143.41 

198753.97 

12/18/2006 

74.60 

1 

150 

PW-6 

456143.41 

198753.97 

12/18/2007 

86.56 

1 
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Table  A-l.  TCE-equivalent  concentrations  used  for  Fort  Lewis  site  calibration  (continued) 


No. 

Well 

East  (m) 

North  (m) 

Date 

Concentration  (ppb) 

Aquifer 

151 

PW-6 

456143.41 

198753.97 

12/10/2008 

35.99 

1 

152 

PW-6 

456143.41 

198753.97 

12/07/2009 

20.90 

1 

153 

PW-6 

456143.41 

198753.97 

09/20/2010 

18.76 

1 

154 

PW-7 

456168.18 

198700.00 

12/20/2006 

37.79 

1 

155 

PW-7 

456168.18 

198700.00 

12/18/2007 

44.40 

1 

156 

PW-7 

456168.18 

198700.00 

12/10/2008 

30.77 

1 

157 

PW-7 

456168.18 

198700.00 

09/20/2010 

16.57 

1 

158 

PW-8 

456195.19 

198649.27 

12/18/2006 

62.05 

1 

159 

PW-8 

456195.19 

198649.27 

12/14/2007 

38.89 

1 

160 

PW-8 

456195.19 

198649.27 

12/10/2008 

30.81 

1 

161 

PW-8 

456195.19 

198649.27 

12/07/2009 

25.82 

1 

162 

PW-8 

456195.19 

198649.27 

09/20/2010 

18.86 

1 

163 

LC-126 

453981.55 

200510.92 

06/01/1995 

93.00 

2 

164 

LC-126 

453981.55 

200510.92 

06/01/1996 

117.50 

2 

165 

LC-126 

453981.55 

200510.92 

06/01/1997 

115.00 

2 

166 

LC-126 

453981.55 

200510.92 

06/01/1998 

103.00 

2 

167 

LC-126 

453981.55 

200510.92 

06/01/1999 

108.80 

2 

168 

LC-126 

453981.55 

200510.92 

06/01/2000 

81.80 

2 

169 

LC-126 

453981.55 

200510.92 

06/01/2001 

82.30 

2 

170 

LC-40D 

454232.17 

200231.42 

06/01/1995 

15.50 

2 

171 

LC-40D 

454232.17 

200231.42 

06/01/1996 

16.00 

2 

172 

LC-40D 

454232.17 

200231.42 

06/01/1997 

17.30 

2 

173 

LC-40D 

454232.17 

200231.42 

06/01/1998 

15.90 

2 

174 

LC-40D 

454232.17 

200231.42 

06/01/1999 

11.80 

2 

175 

LC-40D 

454232.17 

200231.42 

06/01/2000 

4.00 

2 

176 

LC-40D 

454232.17 

200231.42 

06/01/2001 

13.30 

2 

177 

LC-66D 

454815.25 

200223.18 

06/01/1995 

42.00 

2 

178 

LC-66D 

454815.25 

200223.18 

06/01/1996 

53.00 

2 

179 

LC-66D 

454815.25 

200223.18 

06/01/1997 

53.00 

2 

180 

LC-66D 

454815.25 

200223.18 

06/01/1998 

57.80 

2 

181 

LC-66D 

454815.25 

200223.18 

06/01/1999 

26.30 

2 

182 

LC-66D 

454815.25 

200223.18 

06/01/2000 

2.90 

2 

183 

LC-66D 

454815.25 

200223.18 

06/01/2001 

32.00 

2 

184 

LC-67D 

454256.85 

199869.25 

06/01/1995 

50.50 

2 

185 

LC-67D 

454256.85 

199869.25 

06/01/1996 

52.30 

2 

186 

LC-67D 

454256.85 

199869.25 

06/01/1997 

56.50 

2 

187 

LC-67D 

454256.85 

199869.25 

06/01/1998 

61.80 

2 

188 

LC-67D 

454256.85 

199869.25 

06/01/1999 

51.80 

2 

189 

LC-67D 

454256.85 

199869.25 

06/01/2000 

51.30 

2 

190 

LC-67D 

454256.85 

199869.25 

06/01/2001 

55.70 

2 

191 

LC-69D 

454756.91 

199682.96 

06/01/2001 

150.00 

2 

192 

LC-72D 

453770.55 

200173.05 

06/01/1995 

45.50 

2 

193 

LC-72D 

453770.55 

200173.05 

06/01/1996 

47.80 

2 

194 

LC-72D 

453770.55 

200173.05 

06/01/1997 

54.00 

2 

195 

LC-72D 

453770.55 

200173.05 

06/01/1998 

57.30 

2 

196 

LC-72D 

453770.55 

200173.05 

06/01/1999 

43.80 

2 

197 

LC-72D 

453770.55 

200173.05 

06/01/2000 

13.60 

2 

198 

LC-72D 

453770.55 

200173.05 

06/01/2001 

37.00 

2 

199 

LC-73D 

453627.71 

199977.88 

06/01/1995 

25.50 

2 

200 

LC-73D 

453627.71 

199977.88 

06/01/1996 

34.30 

2 
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No. 

Table  A-l. 

Well 

TCE-equivalent  concentration  selected  in  Fort  Lewis  (continued) 

East  (m)  North  (m)  Date  Concentration  (ppb) 

Aquifer 

201 

LC-73D 

453627.71 

199977.88 

06/01/1997 

34.30 

2 

202 

LC-73D 

453627.71 

199977.88 

06/01/1998 

28.80 

2 

203 

LC-73D 

453627.71 

199977.88 

06/01/1999 

20.60 

2 

204 

LC-73D 

453627.71 

199977.88 

06/01/2000 

13.50 

2 

205 

LC-73D 

453627.71 

199977.88 

06/01/2001 

21.30 

2 

206 

LC-74D 

453425.18 

199565.97 

06/01/1995 

40.00 

2 

207 

LC-74D 

453425.18 

199565.97 

06/01/1996 

50.30 

2 

208 

LC-74D 

453425.18 

199565.97 

06/01/1997 

54.50 

2 

209 

LC-74D 

453425.18 

199565.97 

06/01/1998 

53.50 

2 

210 

LC-74D 

453425.18 

199565.97 

06/01/1999 

72.70 

2 

211 

LC-74D 

453425.18 

199565.97 

06/01/2000 

68.00 

2 

212 

LC-74D 

453425.18 

199565.97 

06/01/2001 

62.30 

2 

213 

PS-13 

453382.38 

199122.79 

06/01/2001 

0.50 

2 

Notes: 

1.  When  a  measuring  date  is  not  clear  or  an  annual  average  is  used,  the  date  set  to  the  middle  of  the  year,  e.g.,  June  1 

2.  Aquifer  1  and  2  indicate  Vashon  and  SLA,  respectively. 


Table  A-2.  Source  flux  measurement  at  EGDY  in  For  Lewis 


Source 

Date 

Source  flux  (kg/day) 

1 

11/01/2003 

0.754 

1 

04/01/2006 

0.005 

2 

11/01/2003 

0.323 

2 

04/01/2006 

0.003 

3 

04/01/2006 

0.420 

Notes: 

1.  Nov  1,  2003  is  pre-ERH  for  Areas  1  to  3  and  Apr  1,  2006  is  post-ERH  for  Areas  1  and  2  and  pre-EREI  for  Area  3. 

2.  The  model  uses  numeric  years,  e.g.,  08/04/2004  =  2004.59,  etc. 


Table  A-3.  DNAPL  mass  removed  by  TSR  at  EGDY  in  For  Lewis 


Source 

Date 

Mass  removal  (kg) 

1 

08/04/2004 

2990 

2 

08/05/2005 

1340 

3 

01/26/2007 

1120 

Note:  The  model  uses  decimal  years,  e.g.,  08/04/2004  =  2004.59,  etc. 
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